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ABSTRACT 

A  procedure  is  presented  that  allows  an  engineer  to  coordinate 
protective  devices  consisting  of  fuses,  circuit  breakers,  and  over- 
current  relays  based  upon  load  and  fault  currents.   The  resulting 
digital-computer  program  makes  many  of  the  decisions  requiring  very 
little  prior  in-depth  knowledge  of  protective  devices  and  power  systems 
for  its  utilization. 

The  coordination  procedure  consists  of  six  separate  programs  . 
that  must  run  in  sequence,  as  the  output  from  one  program  forms  the 
input  for  the  next  program.   This  reduces  in  size  of  internal  computer 
storage  required. 

The  theory  is  presented  with  an  explanation  of  each  algorithm 
followed  by  a  basic  flow  chart.   An  example  computer  run  is  presented 
at  the  end. 
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1.   INTRODUCTION 

A  prime  objective  of  electrical  engineers  working  on  the  design, 
operation,  or  maintenance  of  electrical  power  systems  is  to  provide 
reliable,  efficient  power  to  their  customers.   The  distribution  system 
has  always  been  a  weak  link  in  the  power  chain.   Often,  the  catas- 
trophic results  of  a  fault  are  easily  seen,  but  the  cause  may  be 
extremely  difficult  to  find.   A  distribution  system,  either  in  the 
utility  or  the  user's  facilities  that  is  not  adequately  protected  from 
electrical  faults  is  both  dangerous  and  costly.   An  unsatisfactorily 
protected  system  can  cause  considerable  inconvenience  and  damage  to 
equipment.   Adequate  protection  is  not  a  set  definition.   Each  load 
system  is  different  with  its  own  set  of  boundary  values. 

In  a  distribution  system,  the  flexibility  of  design  reduces 
near  the  generator  or  the  load.   Usually  the  generator  or  load  protec- 
tion is  designed  as  an  integral  part  of  the  equipment  by  the  manufacturer 
and  any  changes  should  only  be  undertaken  with  their  counsel.   Moving 
away  from  the  generator  or  load,  the  priority  shifts  from  load  pro- 
tection to  distribution  protection.   It  must  never  be  forgotten  that 
the  load,  distribution,  and  generation  systems  must  be  protected  in 
such  a  way  that  a  device  that  operates  satisfactorily  in  one  system 
cannot  have  a  detrimental  effect  on  another  system.   In  this  thesis 
the  distribution  section  is  attached  from  an  infinite-bus  power 
source  to  the  load-protection  device.   The  programs  provide  the 


proper  setting  of  a  chosen  device  in  order  to  provide  good  protec- 
tive coordination.   Proper  coordination  is  achieved  when  a  device 
closest  to  the  fault  clears  the  fault  or  overload  condition  before 
the  rest  of  the  system  is  disturbed.   In  order  to  achieve  this  end  by 
manual  means,  extensive  mathematical  calculations  and  experience  are 
required.   As  a  direct  result, proper  coordination  is  often  not  achieved. 
The  use  of  digital  computers  to  solve  the  coordination  problem  makes 
a  more  exacting  solution.   However,  computers  with  large  internal 
memory  are  required  to  solve  these  problems.   Several  approaches 
to  computer  solutions  have  been  published  [1-4]. 

Problem  Statement 

The  purpose  of  this  thesis  is  to  develop  the  procedure  and  an 
interactive  computer  program  for  a  generalized  case  of  device  coordina- 
tion.  This  coordination  should  include  fuses,  relays,  and  molded  case 
circuit  breakers.   The  phases  of  research  are  as  follows. 

1.  Develop  a  computer  procedure  to  compute  the  one-line  program 
per  unit  values. 

2.  Develop  a  computer  procedure  to  perform  a  load-flow  analysis. 

3.  Develop  a  computer  procedure  to  compute  line-to-line  ground 
and  three-phase  faults. 

4.  Develop  and  apply  procedures  for  mathematical  representation 
of  line  X/R  data,  fuse,  molded-case  circuit  breaker,  and 
relay-curve  time-current  characteristics. 


5.   Develop  a  computer  procedure  to  coordinate  the  devices 
from  phase  4  with  the  data  derived  from  phases  2  and  3. 

Although  the  coordination  programs  are  written  for  the  PDPIO 
computer,  it  requires  no  more  than  25K  bytes  of  internal  memory  and 
a  disk  unit  to  solve  any  system  of  30  buses  and  100  elements.   Almost 
any  large  system  can  be  reduced  to  its  smaller  parts  composed  of  30 
buses  or  less.   As  each  segment  is  coordinated,  it  can  be  represented 
by  a  load  and  its  protective  device.   In  this  manner  a  300-bus  system 
could  be  divided  into  smaller  segments  and  in  11  runs  would  be  com- 
pletely coordinated.   The  step-by-step  approach  used  here  allows  for 
insertion  or  extraction  of  data  before  proceeding  to  the  next  step. 
Figure  1.1  is  a  basic  flow  chart  of  the  coordination  program.   Here- 
after, it  will  be  referred  to  by  its  individual  parts. 

When  the  coordination  process  has  been  completed,  a  load  flow 
study  and  a  fault  analysis  has  been  performed  as  well.   The  major 
assumptions  that  have  been  made  are  listed  here.   Other  minor  ones 
are  included  in  the  corresponding  chapters.   These  assumptions  have 
been  made  in  an  effort  to  keep  the  solution  as  general  as  possible. 

1.  The  incoming  line  is  assumed  to  be  an  infinite  bus. 

2.  Assymetrical  currents  are  neglected.   Devices  are  coordi- 
nated on  the  basis  of  symmetrical  fault  currents. 

3.  Only  synchronous-motor  contributions  to  the  fault 
currents  are  recognized.   If  a  large  induction  motor 
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is  present,  it  may  be  added  as  a  fictitious  synchronous 
motor. 

4.  Line-charging  shunt  admittance  is  ignored. 

5.  Prefault  current  and  voltages  are  neglected  in  fault- 
current  computations. 

6.  Different  sizes  of  fuses  have  the  same  time-current  curve 
shape. 

7.  Overload  conditions  are  125%  of  full  load. 

8.  The  system  has  a  balanced  three-phase  load. 

In  this  thesis  each  phase  of  research  in  the  problem  statement 
has  a  chapter  devoted  to  it.   After  describing  the  computational 
methods  used  in  curve  smoothing, the  coordination  problem  begins  with 
data  input.   Next  the  fault  currents  are  solved  followed  by  the 
actual  coordination  program  and  device  response  graphing.   Finally, 
a  sample  distribution  system  coordination  problem  is  described. 


2.   CURVE  SMOOTHING  AND  DIGITIZING 

General 

In  this  chapter,  the  methods  for  using  a  data  digitizing  "numonics' 
device  and  the  curve  smoothing  program  of  Appendix  A  are  explained. 
The  various  device  equations  that  were  found  that  adequately  represented 
the  particular  device's  time  current  curve  are  also  described. 

A  method  was  described  in  Wagner  [1]  for  performing  this  by 
using  the  PDPIO  computer.   This  least-square  curve-fitting  prograjn  [1] 
was  expanded  to  eliminate  any  external  subroutines  that  may  not  be 
available.   The  coordination  program  makes  extensive  use  of  polynomial 
functions  whose  coefficients  were  derived  by  using  the  least-sqnares 
curve-smoothing  technique. 

Curve  Smoothing 

Electrical  Cables.   Coefficients  for  the  mathematical  modeling 
of  resistance  and  reactance  were  taken  from  data  for  mine  power 
cable  [5].   This  is  used  as  only  one  type  of  cable  or  wire.   Any  other 
type  may  be  modeled  by  the  same  procedure.   For  instance,  either  actual 
cable  or  line  data  could  be  used  directly  [6], or  the  theoretical  resis- 
tance and  reactance  for  each  size  could  be  calculated  [5-8]. 


Fuses.   Two  types  of  fuses  are  used  in  the  coordination  pro- 
gram.  These  are  current-limiting  and  solid-material  boric-acid  power 
fuses.   In  each  case, the  curve  shape  is  assumed  to  be  the  same  for 
different  sizes.   The  low-current  asymptotic  value  is  called  the 
"of set."   Each  fuse  has  a  minimum  and  a  typical  operating  time.   In 
fuse  coordination, both  values  must  be  used  in  describing  the  operating 
envelope  of  a  particular  fuse.   The  higher  typical  operating  time 
curve  is  referred  to  as  "Of set  1."   By  this  means,  any  fuse  can  be 
completely  described  by  four  equations,  where: 


n 


Log  T  .   =  (  E  a.  (Log  Current -Of set)^  -  2  (2.1) 

mm     .    i 
j=o  -J 


Log  T   =  (  Z  a, (Log  Current -Of setl)  -  2  (2.2) 

yp    k=0  ^ 


Fuses  are  selected  also  by  using  mathematical  models.   To 
select  a  fuse,  the  program  uses  the  overload-current  value  as  the 
variable  in  the  polynomial  function.   The  result  is  the  device  number. 
In  curve  smoothing,  variables  can  assume  no  value  less  or  greater  than 
those  originally  used  as  data  for  the  curve-smoothing  program,  or  large 
errors  can  result. 

Molded-Case  Circuit  Breakers.   Like  the  fuses  discussed  above, 
each  molded-case  circuit  breaker  consists  of  an  operating  envelope 
bounded  by  a  minimum  value  and  a  maximum  value.   Unlike  fuses,  however. 
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a  table  look-up  approach  is  used  in  selection  of  a  device  and  its 

setting  because  of  the  smaller  number  of  sizes  and  magnetic-trip  settings, 

Relays .   Wagner  [1]  discusses  the  representation  of  the  family 
of  relay  curves  for  each  type  with  one  polynomial  equation  comprised 
of  four  coefficients.   The  0   power  coefficient  was  left  out  as  a 
variable  for  coordination.   Unfortunately,  this  method  proves  unsatis- 
factory when  working  with  more  than  one  relay  type.   It  was  found  that 

t"  v» 
each  relay  could  be  represented  by  five  coefficients.   The  0   power 

coefficient  is  necessary  in  that  it  provides  the  time  displacement 
necessary  for  minimum  time  dial  setting.   The  curves  were  derived 
using  the  procedure  below. 

Westinghouse  Company  type  relays  were  used.   An  average  curve 
was  taken  as  time-dial  4.   This  curve  was  reduced  to  a  polynomial 
function.   The  time  difference  between  time-dial  4  and  time-dial  1/2 
taken  at  20  times  the  pickup  value  is  subtracted  from  the  0    coef- 
ficient.  Twenty  times  pickup  was  a  good  value  for  computing  the  time 
difference  as  the  time  value  changes  very  little  for  increases  in 
current  beyond  that  point.   The  resultant  curve  was  the  average  shape 
of  a  time-dial  setting  of  4  now  shifted  to  minimum  setting.   This  allows 
for  time  values  computed  during  coordination  to  be  summed  with  the  0 
coefficient.   The  curves  are  then  represented  as  follows: 


i 
Log  Time  =  (    a . (log(Current  Value)  )  + 

j=0  ^ 


Log  (Time  Setting) 


n 


This  allows  the  use  of  more  than  one  type  of  relay  in  each 
coordination  as  it  separates  time  setting  value  from  relay  type. 

Curve  Digitizing.   The  device's  time  versus  current  curves 
were  digitized  by  using  a  NUMONICS  electronic  graphics  calculator 
(EGC).   The  NUMONICS  EGC  is  designed  to  translate  graphic  information 
on  any  medium  into  digital  information.   Complete  information  on  its 
usage  may  be  found  in  -NUMONICS  Electronic  Graphics  Calculator  [9]. 
The  NUMONICS  EGC  then  provided  a  convenient  way  to  enter  X-Y  coordinate 
data  representing  any  curve  onto  disk.   Since  this  device  works  in 
length  in  inches,  the  data  has  to  be  normalized  to  its  proper  units. 

A  simple  program  DIGIT. F4  was  written  to  normalize  and  store 
onto  disk  the  NUMONICS  EGC  information.  The  program  flow  chart  is 
shown  in  Figure  2.1. 

Curve  Smoothing.   The  curve-smoothing  program  contained  in 
Wagner  [1]  was  used  as  a  basis  for  the  curve  smoothing  program  DEV.F4. 
Unnecessary  parts  were  removed  and  a  double  precision  matrix  inver- 
sion routine  was  added.   This  prevented  the  necessity  of  using 
subroutines  contained  in  the  computer's  SSP  routine  files.   The 
inversion  routine  was  converted  from  The  Pennsylvania  State  University 
Computer  Center's  matrix  inversion  routine  CMINV  which  inverts  a 
complex  square  matrix.   This  routine  in  turn  was  converted  from  the 
SSP  routine  MINV  [10] • 

Since  it  has  its  own  inversion  routine,  this  curve  smoothing 
program  can  be  easily  expanded  to  any  number  of  points  and  any 
dimension  of  fit  by  changing  the  dimension  statements.   The  flow 
chart  and  procedure  remains  the  same  as  in  [1] .   A  complete  program 
listing  may  be  found  in  Appendix  A. 
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Figure   2.1.      Continued. 
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Summary 

Polynomial  functions  describe  the  various  device-response 
curves,  line  X/R  data,  and  fuse-selection  criteria.   The  coefficients 
for  these  polynomial  functions  are  found  from  data  entered  using  a 
digital-data  device  or  manual  input.   This  data  is  reduced  by  a  curve- 
smoothing  program  to  the  desired  function. 

The  start  of  the  coordination  problem  begins  in  the  next  chapter 
with  line  and  load  data  entry. 
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3.   THE  INPUT  OF  DISTRIBUTION-SYSTEM  DATA 

General 

Usually,  the  first  step  in  performing  any  calculations  on  a 
distribution  system  is  to  reduce  all  system  values  to  per  unit.   The 
methods  used  to  do  this  in  the  program  will  be  described  in  this 
chapter.   Errors  can  be  introduced  easily  when  converting  to  the  per 
unit  system  using  manual  calculations.   Program  INPUT. F4  was  written 
to  perform  this  function  as  well  as  computing  line  resistance  and 
reactance  values  and  storing  all  system  data  on  disk  for  use  in  sub- 
sequent programs. 

Per-Unit  Reduction 

The  program  uses  the  following  formulas  to  convert  to  the  per- 
unit  system  [8].   (See  Table  3.1  for  variable  definitions  used  in  this 
thesis.  ) 


Base      Base      Base 


Z^  -,  =  Z(fi)/Z^ 
P.U.         Base 


I_     =  KVA^    /KV^ 
Base      Base    Base 
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Table  3.1.   Variable  definitions, 


Z{Q.)  element  impedance 


Z  impedance  in  per  unit 


X,         subtransient  reactance 
a 


P  real  power  bus  p 

P 


Q  reactive  power  bus  p 


E  ,         complex  voltage  bus  p 


I  complex  current  bus  p 

P 


Y  Y  bus  element  bus  p  to  q 

pq 


y  admittance  value  from  primitive  admittance  matrix 

pq,rs 


Z  Z  bus  element  bus  p  to  q 

pq 


z  impedance  value  from  primitive  impedance  matrix 

pq,pq 


z_         fault  impedance 


y  fault  admittance 
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KV„      .     2  KVA^ 
_  ,      Base  given.  ,    Base  new  . 

P.U.new    P. U. given  ^  KV„        '   ^KVA^      .    ' 

Base  new       Base  given 


Utilizing  the  above  formulas  it  is  possible  then  to  represent  all 
elements  in  the  distribution  system  on  a  per-unit  basis.   Loads  need 
to  be  included  also  in  the  one-line  representation. 

The  program  allows  loads  to  be  represented  either  by  induction 
motor  horsepower,  synchronous  motor  horsepower,  kVA,  or  current. 

Loads  are  entered  in  the  following  format:   Bus  number/horse- 
power, kVA,  current/power  factor.   The  program  assumes  100%  efficiency 
so  one  must  use: 


HP^     J  =  HP/Efficiency  (3.1) 

Entered     '  y 


If  information  on  efficiency  and  power  factors  are  unknown,  the  follow- 
ing equations  may  be  used  instead  of  3.1  [9]. 

For  induction  motors,  use: 


HP^  ^    J  =  HP  X  1.1627 
Entered 


Power  Factor  =  .86  (3.2) 

For  synchronous  motors  use : 

Unity  Power  Factor  HP  =  1.139  x  Horsepower  (3.5) 

.8  Power  Factor  UP  =  1.18  x  Horsepower  (3.6) 
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KV„      .     2  KVA„ 
,        _  ,      Base  given,  ,    Base  new  . 

^P.U.new  "   P. U. given   KV„        '     KVA„      .    ^ 

Base  new       Base  given 
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Loads  are  entered  in  the  following  format :   Bus  number/horse- 
power, kVA,  current/power  factor.   The  program  assumes  100%  efficiency 
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If  information  on  efficiency  and  power  factors  are  unknown,  the  follow- 
ing equations  may  be  used  instead  of  3.1  [9]. 

For  induction  motors,  use: 


HP^     ^  =  HP  X  1.1627 
Entered 


Power  Factor  =  .86  (3.2) 

For  synchronous  motors  use : 

Unity  Power  Factor  HP  =  1.139  x  Horsepower  (3.5) 

.8  Power  Factor  HP  =  1.18  x  Horsepower  (3.6) 
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Formulae  for  Line  Resistance  and  Reactance 

Element  data  is  entered  in  the  following  format:   line  bus/ 
load  bus/cable  length  in  feet/wire  size/transformer  P.U.  impedance/ 
size  in  MVA/line  voltage  in  KV/load  voltage  in  KV. 

Wire  sizes  that  are  entered  in  American  wire  gauge  (AWG)  are 
converted  to  MCM  wire  gauge  by  using  the  following  formula  [1] : 


MCM  =  [(.005)  [(1.229)  ^-^^  ^^'^''*]   xlOOO  (3.7) 


Once  the  wire  size  is  in  MCM,  the  line  resistance  and  reactance  are 
computed  using  the  polynomial  coefficients  that  were  found  by  the 
procedure  in  Chapter  2. 


Resistance  or  reactance  =   I  a.[log(MCM)]  (3.8) 


j=o  J 


The  program  should  not  be  limited  to  just  one  wire  type.   The  program 

does  this  by  reading  wire  coefficients  a.  above  from  disk  so  that  each 

J 

time  the  program  is  run  it  can  compute  element  data  using  different  wire 
types  until  the  entire  system  has  been  completely  entered. 

Transformers  become  very  important  in  that  they  can  become  the 
highest  source  of  resistance  and  reactance.   By  their  connections  and 
neutral  impedance,  they  determine  the  magnitude  of  zero-sequence  current 
that  will  flow  in  the  particular  element.   The  program  allows  for  this 
input  and  writes  this  information  in  the  particular  disk  files  for  use 
later  on  in  fault  analysis. 
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Transformer  inrush  and  withstand  values  must  also  be  calcu- 
lated.  These  values  must  be  used  in  coordinating  protective  devices. 
These  devices  must  pass  inrush  currents  but  open  under  conditions  of 
exceeded  transformer  withstand  currents.   The  algorithm  in  the  program 
is  based  on  average  values.   If  test  data  is  available  this  should  be 
entered  directly  by  changing  the  transformer  data  disk  file  directly 
after  the  program  is  run  and  before  proceeding  with  coordinations. 
Inrush  and  withstand  values  are  computed  as  follows  [11] : 


Inrush  =  12x (Transformer  Base  Current)  (3.9) 


Withstand  (32  Sec  =  25x  (Transformer  Base  Current)       (3.10) 


Withstand  @5  Sec  =  14 . 5x  (Transformer  Base  Current)      (3.11) 

Transformer  withstand  is  some  value  between  the  two  second  value  to  the 
five  second  value  and  depends  upon  the  transformer  impedance.   It  is 
fairly  easy,  however,  to  have  the  device  protect  the  complete  range 
of  withstand  values  and  the  program  is  written  to  accomplish  this  task. 
Synchronous  subtransient  reactance  is  computed  as  follows  [11] : 


(Per  unit  KVA)  x  (.17)  =  iX,  (3.12) 

d 
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One  of  the  problems  encountered  in  computing  transmission  line 
impedance  is  a  reasonable  calculation  of  the  zero-sequence  impedance. 
The  computer  program  assumes  that  this  value  is  two  times 
the  positive-sequence  impedance  [8]  if  nothing  is  entered  when  the 
computer  asks  for  the  zero-sequence  impedance  multiple.   If  a  different 
value  is  entered,  the  computer  multiplies  the  positive  sequence  reactive 
value  by  this  multiple  to  arrive  at  the  zero-sequence  value. 

A  basic  flow  chart  follows  in  Figure  3.1  with  a  complete  program 
listing  in  Appendix  B. 

Summary 

The  first  step  toward  solving  the  coordination  problem  is  com- 
pleted.  All  element  data  has  been  converted  to  the  per-unit 
representation.   The  next  logical  step  is  to  perform  a  load-flow 
analysis  to  determine  the  full  load  currents  in  each  element. 


Flow  Chart  Input 


Input 


Enter  base 
MVA,  kV 


^ 


Read  wire 
coefficient 


Read 

element 

data 
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Enter 
element  data 


Compute 
here  per 
unit 
values 


Enter  bus 
voltages 


Figure  3.1.   Per-unit  calculations  flow  chart, 
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Input 


Enter  loads 


Compute 
P.U.  values 
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Compute 
Xd 


Write 
"^^s  /   data  to 
disk 


Qtoi 


Figure  3.1.   Continued. 
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4.   LOAD- FLOW  ANALYSIS 

General 

In  order  to  establish  overload  conditions  on  each  element, 
the  load  currents  must  be  determined.   This  is  done  by  using  a  load 
flow  analysis.   The  Newton-Raphson  Method  using  Y-bus  [12]  was  chosen 
because  it  required  the  least  amount  of  iterations  to  obtain  a  solu- 
tion to  the  load -flow  problem. 

Load  Flow 

In  this  method  Y-bus  must  first  be  formed  from  the  primitive- 
impedance  network.   When  mutual  coupling  is  ignored,  Y-bus  can  be 
formed  easily  by  algorithm.   Any  system  can  be  completely  represented 
by  a  set  of  equations  [8,12]. 


^bus  ^  ^^bus^^^bus^  ^^_^^ 


^bus    ,  ^bus        ,         ...  r-       1  1 

E    and  I    are  column  matrices  with  one  entry  for  each  bus 

in  the  system.   Y     is  a  square  matrix  of  order  equal  to  the  >iumber 

of  buses  in  the  system.   Therefore,  once  all  bus  voltages  are  xnown  and 

Y-bus  is  formed,  matrix  multiplication  gives  the  desired  load  currents. 
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However,  there  is  no  direct  approach  to  obtain  the  bus  voltages, 
The  only  approach  to  use  is  that  of  iteration.   The  Newton-Raphson 
method  accelerates  the  iteration  process  so  that  normally  an  accurate 
answer  is  obtained  after  three  iterations. 

The  power  at  bus  p  is  (from  4.1): 


P  -jQ  =  E  I 
P   P    P  P 


P  -jQ   =  e''  Z  Y   E  (4.2) 


So 


P 


p  V    I  pq  q 


E   =  e  +jf  Y    =  G   -IB  (^-3) 

P     p    p  pq     pq    pq 


P   =   E   {e  (e  G  +f  B   )+(f  G   -e  B   )}  (^.^) 

P    ^^i         P   q  pq   q  pq     q  pq   q  pq 


Q_,  -   Z   {f  (e  G  +f  B   )-e  Cf  G   -^  R   ^^ 
P    q  =  l    P   q  pq   q  pq^   P^  q  pq  ^%^  ^ 


(A. 5) 


The  Newton-Raphson  method  uses  the  Jacobian  matrix  of  partial- 
differential  equations  as  follows: 
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n-1      -r 

J9e 
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n-1 
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3e, 


3Q 


de 
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3e 


n-1 
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9Qi 
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9Q 


n-1 
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1 

9P      , 
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1 

3Ql 
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1 

="^n-l 
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1 

9P. 


9f 


n-1 


9P 
Jf 


n-1 
n-1 


9Q, 


3f 


n-1 


3Q 


n-1 


9f 


n-1 


Ae. 


Ae 


n-1 


Af. 


Af 


n-1 


~ 

•■ 

~                 — 

AP 

'l 

h 

A 
e 

AQ 

■'3 

^4 

'^f 

*                                     ^B_ 

L_ 

_ 

^                  P« 

(4.6) 


The  elements  of  the  Jacobian  matrix  are  found  from  the  following 
equations : 


J-,  :  9P 

—-2-  =  e  G   -f  B    q^p 
9e      p  pq   p  pq 


9P 

I 

9e 


=  2e  G   +f  B   -f  B   +  I  (e  G   +f  B   ) 

p  pp  p  pp  p  pp  -^    q  pq  q  pq 

q=p 


(A. 7) 
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I 

8f 


=  e  B   +f  G    q?^p 
p  pq   p  pq  ^ 


8P 

I 

3f 


P  _ 


e  B  +Z^  G   -e  B   +  E  (f  G   -e  B   ) 

p  pp  fp  pp  p  pp   1  q  pq  q  pq 
q?^p 


(4.8) 


8e 


-  e  B   +f  G    q/p 
p  pq   p  pq 


30, 

de 


f  G  -f  G  +Z 
P  PP   P  pp 


ep  PP  q^i 

q^p 


-  L  (f  G   -e  B   ) 
q  pq   q  pq 


(4.9) 


3Q 


P  _ 


e  G  +f  B    q  7^  p 

p  pq  p  pq 


3Q 


P  _ 


=  e  G  +Z 


-e  G 


p  PP  fp  PP  p  PP  ^_i  q  pq 


n 
+  E  (e  G 
q=l 

p^q 


+f  B  ) 

q  pq 


(4.10) 
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Bus  //I  or  the  swing-bus  voltage  is  known.   The  swing  bus  is 
the  bus  that  will  provide  whatever  power  the  system  requires.   The 
program  begins  by  forming  Y-bus  by  using  the  primitive-impedance  matrix 
generated  by  INPUT. F4  of  Chapter  3.   Since  mutual  coupling  is  ignored, 
the  y-admittance  matrix  can  be  found  by  inverting  each  member  of  the 
z-impedance  matrix  found  in  disk  file  TEST. DAT.   Y    matrix  is  formed 
using  the  algorithm  below. 

Diagonal  elements  are  formed  as  the  sum  of  all  admittances 
connected  to  that  bus  or  node.   Off-diagonal  elements  are  the  negative 
of  the  admittance  connected  between  the  buses  or  nodes. 

If  more  than  10  iterations  are  performed,  the  program  stops  and 
informs  the  user  that  either  an  error  has  been  made  in  design  or  data 
entry.   Upon  successful  completion  of  the  program,  a  load-flow  analysis 
may  be  obtained  by  printing  file  DISK. DAT.   The  program  puts  the  load 
currents  in  the  proper  format  and  file  for  later  use  in  the  coordination 
program.   A  flow  chart  follows  in  Figure  4.1,  and  a  complete  listing 
is  given  in  Appendix  B. 

Summary 

This  load-flow  analysis  obtained  can  be  very  useful  in  analyzing 
the  distribution-system  design.   It  can  be  used  to  point  out  locations 
for  power  factor  correction  as  well  as  potential  problems  due  to  under- 
sized wiring. 

The  load  currents  found  here  are  stored  on  disk  to  be  used  in 
the  coordination  program.   The  next  step  before  beginning  the 
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coordination  program  is  to  solve  for  all  fault  currents.   Z-bus  must 
be  formed  to  compute  these  currents.   The  methods  used  to  do  this 
follows  in  the  next  chapter. 
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Figure  4.1.   Newton-Raphson  load-flow  chart. 
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5.   FORMATION  OF  IMPEDANCE-BUS  NETWORK 

General 

012 
The  formation  of  the  impedance  bus  network  (Z-Bus    )  by 

algorithm  is  the  first  step  toward  solving  for  fault  currents.   In 

this  chapter  this  algorithm  will  be  discussed  as  well  as  other  problems 

encountered  in  matrix  compression. 

Formation  of  Z-Bus 

Once  the  one  line  diagram  has  been  completed,  the  primitive 
impedance  matrix  that  was  found  in  INPUT. F4  is  a  partially  filled 
matrix.   However,  it  is  not  in  the  proper  format  for  the  formation  of 
Z-bus  as  it  contains  no  coupling.   The  primitive  impedance  matrix  must 
change  from  one  of  order  equal  to  the  number  of  buses  to  one  of  order 
equal  to  the  number  of  elements.   This  matrix  becomes  a  very  large 
sparse  matrix.   The  memory  required  to  store  this  matrix  becomes  too 
large  to  fit  in  a  small  computer.   A  method  was  derived  to  only  store 
the  non  zero  values  and  their  location  in  the  matrix.   Since  the  matrix 
is  triangular  only  the  top  half  is  stored  (Figure  5.1). 

Diagonal  elements  found  in  INPUT. F4  are  stored  in  "ZD."   ZD 
is  a  three-dimensional  matrix,  ZD(200,2)  where  (X,l)  is  the  element 
number  and  (X,2)  is  the  diagonal  value. 
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(A) 


ZD(200,2) 
ZD(200,1) 


Elements  X,Y,Z 


ZM( 200,1) 


ZM(200,2) 


Number  locations  matrix  (A)  =  (//  elements)   X  2 

Number  of  locations  of  alternate  method  =  8  X  (#  elements) 


Figure  5.1.   Transformation  of  square  sparse  primitive  impedance, 
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The  off  diagonal  mutual  coupling  Zm  is  a  three  dimensional 
matrix  also  where  (X,l)  is  the  location  value  A+jB.   Element  A  is 
coupled  to  element  B.   Location  (X,2)  is  the  coupling  value.   Fortun- 
ately, partitioning  makes  it  necessary  to  assemble  a  square  matrix  for 
inversion  of  only  the  coupled  elements.   Since  coupling  only  exists  in 
zero  sequence,  this  applies  only  to  the  zero-sequence  Z-bus  matrix. 

In  the  formation  of  the  Z-bus  matrix,  bus  1  is  the  reference 
bus  for  all  calculations.   Z-bus  is  formed  by  starting  with  element 
one  and  proceeding  with  each  element  in  turn.   An  element  is  either  a 
link  or  branch.   A  link  means  that  both  buses  have  been  used  earlier; 
a  branch  means  they  have  not  been  used  earlier. 

Like  Y-bus  in  Chapter  4  any  system  can  be  represented  by 


^,   012    ^,   012  ^,   012 
Ebus    =  Zbus     Ibus 


012 
Formation  of  Zbus     can  be  formed  either  by  the  incidence, 

network-matrices  method  [12]  or  by  algorithm.   Formation  by  algorithm 
is  much  simpler  and  is  used  here. 

After  forming  the  primitive  impedance  network,  it  is  a  straight- 
forward application  of  the  equations  of  Table  5.1. 


33 


Table  5.1.   Equations  for  the  formation  of  ZBUS , 
[11] 


Add  p-q     Mutual  Coupling  No  Mutual  Coupling 


Branch      Z  .  =  Z  .+y   ,  ■  (Z  .-Z  . )  Z  .  =  Z  . 

qi     pi   pq  rs   ri   si  qi    px 


i  =  1,  2 ,  .  ,  .  ,  m 
i  +   q 


1+y   ,   (Z   -Z   ) 
pq  rs   pq   rs 

Z    =  Z   +      ^— ^ Z    =  Z   +z 

qq     pq       ^pq''^^  °^^  P^   P^.P^ 


y  ,  (z  .-z  .) 

Link        Z,  .  =  Z  .-Z  .+  pq  ^"   ^"   ""         Z,  .  =  Z  .-Z  . 
li     pi   qi       y  ll     pi   qi 

pq,pq 


i  =  1,  Z,  .  .  .  ,  m    i  ?^  1 


1+y   ,^s(Z  ^-^   ,) 
^11  "  ^>l"^ql^      ^^  ^11  "  ^pl~^ql 


pq.pq 


Elimination  of  1   Node 


+  z 

pq,pq 


z.^z^. 

Z .. (modified)    =    Z.. (before    elimination)    -  — -— — 


i,j    =   1,    2,    3,    .     .     .     ,    m 
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For  example,  a  primitive  impedance  matrix  looks  like  this 
for  five  elements: 


1(1-2)    2(2-3)    3(3-4)    4(1-4)    5(2-4) 


.l+j.5 
j5 


J5 
.l+j.5 


.l+j.5  jlO  j4 
jlO  .l+j.5  j8 
j4       j8       .l+j.5 


[z  ] 


where  the  off  diagonal  values  represent  mutual  coupling  between  ele- 
ments.  Using  this  primitive  matrix  and  the  equations  in  Table  5.1 
yields  Zbus. 

ZBUS.F4  is  one  of  the  longest  programs  in  the  series.   The 
program  asks  for  the  zero-sequence  impedance  for  all  loads  in  the 
system  and  the  mutual-coupling  impedances.   If  the  zero-sequence 
impedance  of  a  load  is  not  known,  the  program  assigns  it  the  same  value 
as  the  positive  sequence  impedance. 

All  loads  that  are  added  are  links,  and  only  loads  that  are 
added  are  those  of  synchronous  motors.   The  subtransient  reactance 
is  used  instead  of  its  actual  impedance  since  the  Z-bus  is  going  to  be 
used  for  fault  calculations.   Any  link  from  a  bus  to  the  reference  bus 
will  cause  a  fault  current.   Prefault  voltages  and  currents  will  be 
ignored  as  will  all  contributions  from  loads  other  than  synchronous 
motors.   A  flow  chart  follows  (Figure  5.2)  and  a  complete  program  list- 
ing in  Appendix  B. 
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Figure  5.2.   Continued, 
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Figure  5.2.   Continued. 
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Summary 

012 
Zbus     is  formed  by  algorithm  allowing  for  mutual  coupling 

between  elements.   The  next  step  in  the  coordination  procedure  is  to 

solve  for  fault  currents.   This  procedure  follows  in  the  next  chapter. 
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6.   FAULT  CALCULATIONS 

General 

The  final  program  required  before  executing  the  coordination 
program  is  FALT.F4.   This  program  computes  all  the  possible  line-to- 
ground  and  three-phase- fault  currents  and  prepares  all  other  load 
data  for  entry  into  the  coordination  program.   The  procedure  to 
accomplish  this  will  be  discussed  in  this  chapter. 

Fault  Currents.   The  equations  describing  the  fault  currents 
need  to  be  derived  [12]. 

Three-Phase  Fault.   For  a  fault  at  bus  P: 


^P(.),-  =  v-c„VV")-\co) 


012 


(6.1) 


where  U  is  unity  matrix,  E  (0)  is  prefault  voltage,  and  Y   is  the 

P  ^ 

fault  admittance  matrix  for  a  three-phase  fault. 


012 


0  0  0 
0  10 
0   0    1 


(6.2) 


Substituting    (6.2)    into    (6.1) 
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(6.3) 


Since   Zpp' 


Zpp      and    z      =    1/y 


/  3 


^P(F) 


S-^^PP 


I„   =   0 


^^o 


(6.4) 
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Using  the  same  method  the  following  equations  are  found  [12]. 


^F  "^PP 


(6.5) 


and 


"i(F)   =  ^  3  ^hm 


^ip'^EpW 


z^+Z 
f   pp 


(1) 


i^P   E.°''=0   E  ,,  °''=0 
1  P(f) 


(6.6) 


These  equations  are  identical  to  those  found  for  a  three-phase 
grounded  fault  [12] . 

The  line-to-ground  fault  admittance  matrix  is  given  by  Y 

012 
below.   Substituting  Y      into  (6.1)  the  equations  below  result  [12 

F 


012 


^012  ^  ^ 


1  1 
1  1 
1    1 


1 


012 


/  3  E 


P(0) 


Z°-.2z(^)4-3z, 
PP    PP     F 


(6.7) 
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012 


/  3  E 


PP      PP       F 


-z  <") 

PP 

Z     +Z     +j2t^ 
PP      PP       F 


-Z 


(1) 


PP 


(6.8) 


012  ^ 
i(F)      ^i(O) 


ir^P 


0 

/I 
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/  3  E 


P(0) 


(0) 


PP 


+2Z 


(1) 


PP 


+3z. 


(0) 


iP 

Z. 
iP 


(1) 
(1) 


iP 


(6.9) 


The  value  of  z   is  entered  from  the  keyboard  when  asked  by  the  program. 
F 

All  of  the  equations  are  available  now  to  solve  for  all  fault  currents. 
All  post-fault  bus  voltages  are  solved.   Using  these  bus  voltages 
each  element  fault  current  is  determined.   These  currents  are  stored 
onto  disk  so  that  they  may  be  used  in  both  the  coordination  and  graph- 
ing programs. 

Summary 


The  programs  thus  far  have  solved  for  all  load  and  fault  cur- 
rents but  no  coordination  has  been  discussed.   Now  the  coordination 
problem  is  ready  to  be  addressed.   In  the  chapter  that  follows,  the 
protective  devices  for  the  distribution  system  will  be  coordinated. 
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Figure  6.1.   Flow  chart  for  FALT.FA, 


44 


7.   PROTECTIVE-DEVICE  COORDINATION 

General 

In  this  chapter, the  protective-device-coordination  problem 
will  be  discussed,   A  computational  method  is  put  forward  that  allows 
a  digital  computer  to  solve  the  problem. 

Device  Coordination 

One  of  the  best  methods  to  approach  coordination  is  a  step  by 
step  device  selection  and  setting  working  from  the  load  toward  the 
source.   A  load  is  the  best  place  to  start  because  there  is  not  much 
flexibility  for  device  setting  at  a  load.   The  load  flow  analysis  of 
Chapter  4  provides  the  overload  values  for  the  various  loads.   This 
enables  a  definite  overload  setting  for  each  device  protecting  a  load. 
Since  overload  currents  are  a  primary  concern  at  these  buses,  all  fault 
settings  can  be  placed  at  six  times  the  full  load  currents.   This  should 
allow  for  normal  starting  currents.   Care  must  be  taken  so  that  an 
instantaneous  tripping  element  will  not  open  under  normal  starting 
conditions.   From  the  loads  it  is  a  straightforward  process  to  place 
a  line  side  device  curve  next  to  its  load  side  curve  to  insure  that 
proper  coordination  is  achieved.   By  taking  this  approach,  each  device 
depends  only  on  one  other  device. 
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This  program  is  written  with  the  option  of  entering  two  devices 
in  each  element — one  located  at  each  bus.   From  this  approach, 
obviously  two  possibilities  can  occur.   Either  the  device  is  on  the 
line  side  of  the  element  or  on  the  load  side  of  the  element  (Figure  7.1). 
If  it  is  the  line-side  device,  its  setting  depends  only  on  the  one 
load-side  device.   It  must  be  coordinated  for  the  fault  current  at 
the  load  bus  and  not  the  line  bus  [1] .   Devices  must  be  coordinated 
for  a  minimum  possible  fault  current  which  will  be  discussed  later. 

If  the  device  is  on  the  load  side  and  it  does  not  protect  a 
load,  then  it  must  protect  a  following  bus.  This  device  must  then  be 
coordinated  with  all  line-side  devices  connected  to  the  bus  it  pro- 
tects.  This  load-side  device  must  be  coordinated  for  fault  currents 
at  the  protected  bus. 

Devices  should  be  set  to  (1)  open  under  overload  conditions, 
(2)  open  under  minimum  fault  conditions  with  proper  coordination,  and  (3) 
open  under  maximum  fault  conditions  with  proper  coordination.   Often, 
however,  all  three  conditions  cannot  be  met  simultaneously  unless  the 
proper  devices  are  chosen.   Criterion  (1)  is  not  as  important  as 
Criterion  (2)  or  (3)  when  in  the  distribution  system  not  at  a  load  bus 
since  an  overload  condition  below  fault  conditions  should  not  do  damage 
to  the  wiring.   For  this  reason  wire  sizes  must  be  checked  to  insure 
that  upon  final  coordination  they  can  withstand  the  overload  and  fault 
conditions  that  can  occur  [5,13].   At  this  point  it  may  be  necessary 
to  begin  the  entire  coordination  procedure  again  if  an  element  cannot 
withstand  the  short-circuit  conditions.   Allowable  short-circuit  cur- 
rents can  be  checked  by  using  the  final  coordination  graph  and  the  fault 
currents  obtained  in  FALT.F4. 
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Bus  1  Bus  2 

(a)   Coordinating  line  side  device  of  an  element 


A 


B 


■o- 


E 
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(b)   Coordinating  load  side  device  of  an  element  feeding  a  bus  with 
load  elements  connected  W,  Z.  Y. 


Figure  7.1.   Two  possibilities  for  coordination  element  Z.   In 

Figure  (7.1a)  coordination  device  A  according  to  curve 
device  B  up  to  fault  bus  2.   In  Figure  (7.1b)  coordina- 
tion device  B  depending  on  load  flow  current  element  Z. 
Coordinated  with  device  curves  C,  D,  E,  up  to  fault  at 
Bus  2. 
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This  coordination  routine  utilizes  fault  currents  and  load 
currents  in  achieving  protective-device  coordination.   The  load  cur- 
rents determine  the  overload  or  pickup  setting  in  each  device.   The 
coordination  problem  is  then  in  two  parts. 

Devices  must  be  coordinated  for  overload  conditions.   This 
provides  the  most  important  coordination  step  for  fuses  in  that  they 
contain  no  variable  fault-protection  settings  as  do  relays  or  molded- 
case  circuit  breakers.   Also,  another  reason  is  that  based  upon  only 
fault  currents  and  load-side  device  settings,  it  is  impossible  to 
satisfactorily  achieve  coordination. 

This  inability  to  achieve  coordination  enters  when  setting  a 
device  that  is  to  protect  a  bus  with  other  elements  attached  to  that 
bus  as  loads.   Just  laying  a  line  curve  next  to  a  load  curve  with  the 
highest  time  value  at  any  current  value  will  provide  an  erroneous 
answer.   This  device  must  be  set  by  determining  its  overload  current 
value  first  since  this  current  is  the  sum  of  all  currents  in  the  load 
elements  that  are  attached  to  that  bus.   Once  this  point  has  been 
established,  each  load  device  curve  must  be  polled  in  turn  to  insure 
that  the  device  setting  responds  to  a  fault  in  a  load  element  slower 
than  the  load  element  line  side  device.   The  device  needs  to  be 
coordinated  up  to  the  maximum  fault  current  that  the  device  will  see. 
It  is  possible  for  this  device  to  see  a  smaller  fault  current  than  that 
in  a  load  element. 

Coordinating  a  device  located  at  a  line-side  bus  is  much  simpler 
since  the  pickup  or  overload  setting  is  the  same  as  the  load-side 
device.   The  line-side  device  must  be  coordinated  up  to  the  maximum 
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fault  current  that  will  occur  at  the  bus  where  its  load-side  device 
is  located.   A  higher  fault  current  will  occur  the  closer  one  gets  to 
the  line-side  device;  however,  this  higher  current  value  will  only  cause 
the  line-side  device  protecting  the  element  to  respond  quicker.   The 
method  of  coordination  will  be  to  first  find  a  load  point  based  on  a 
load  flow,  and  second  coordinate   the  device  at  any  current  value  up 
to  some  maximum  fault  current  such  that  the  line-side  device  has  a 
longer  response  time  than  its  load-side  device. 

Since  the  program  coordinates  three  different  types  of  devices, 
nine  different  possibilities  can  occur.   Each  device  may  be  coordinated 
with  one  of  its  own  type  or  one  of  the  other  two. 

The  program  structure  lends  itself  well  to  the  use  of  subrou- 
tines. Two  of  these  are  SELECT  and  SETDEV.  SELECT  chooses  a  device 
and  sets  its  own  load  point,  and  SETDEV  coordinates  it. 

Fuse  Selection  and  Coordination 

Due  to  the  large  number  of  fuse  sizes  available,  look-up  tables 
are  not  as  efficient  as  using  polynomial  functions  to  select  devices. 
Once  a  load  current  is  known,  this  value  is  substituted  into  one  func- 
tion providing  the  proper  device  number.   This  device  number  is  next 
substituted   into  two  other  functions  providing  the  OFSET  and  OFSET  1 
values  referred  to  in  Chapter  2.   Since  there  are  two  types  of  fuses 
in  the  program,  current-limiting  and  solid-material  boric-acid  power 
fuses,  there  are  six  functions  used  in  selecting  fuses.   These  all  have 
the  following  form: 
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I    a.X^  (7.1) 

j=0  J 


X  will  be  the  logarithm  of  overload  current  or  device  number.   Table  7.1 
lists  the  various  coefficients  a.. 

After  the  fuse  is  selected,  it  must  be  coordinated.   The  program 
uses  the  values  of  OFSET  and  OFSET  1  with  the  applicable  fuse  curve  to 
solve  for  a  time  value  for  a  set-current  value.   The  program  searches 
for  the  load-side  device  to  find  the  highest  time  value  to  that  current 
value.   If  the  time  value  is  larger  for  the  size,  it  is  properly  coor- 
dinated.  If  not,  it  increases  the  fuse  size  by  one,  solves  for  OFSET 
and  OFSET  1,  and  tries  again  continuing  until  it  has  found  the  proper 
fuse  size. 

Relay  Selection  and  Coordination 

Relays  provide  the  most  difficult  coordination  problem.   First 
a  turns  ratio  for  the  current  transformer  (CT)  must  be  chosen.   The 
computer  selects  the  proper  tap  setting  to  prevent  saturation.   It 
does  this  by  always  selecting  the  highest  turns  possible  while 
selecting  the  tap  setting  on  the  relay  as  well  [12].   The  computation 
uses  a  600/5  CT  [13]  with  tap  settings  and  CO  relays  [14]  with  tap 
settings.   This  combination  of  tap  setting  for  the  CT  and  the  relay 
combined  with  the  load  current  establish  the  pickupi  point: 

Pickup  = Log  [(Turns  Ratio  CT  x  Turns  Ratio  Relay)]     (7.2) 


Table  7.1.   Fuse-selection  coefficients, 
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Current  to  Device 


Offset 


Offset  1 


Current 

Limiting  Fuses 

^0 

-  11.15221 

^1 

23.865 

^2 

-  16.487 

^3 

2.9288 

^ 

3.19672 

S 

-  1.48877 

a^ 

.1776431 

.9532 

.046433 

.117885 

.029352 

.00337 

.0001847 

.0000039 


1.1619 

■  .110148 
.10626 

■  .02249 
.002452 

■  .0001365 
.00000304 


Solid  Material  Boric  Acid 

a  88.1643 

a  -215.2235 

a  188.225 

a  -  69.13403 

a,  7.4695 
4 

a  1.65757 

a,  -    .347786 


1.34359 

1.41734 

.151393 

.1585889 

.0239735 

-  .035007 

.00046627 

.007180 

.0004321 

.0008126 

.0000182 

-  .0000377 

.0000078 

0 
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Since  relay  curves  are  in  units  of  multiples  of  pickup  this 
pickup  value  will  be  used  often.   Only  the  relay  pickup  has  been  set, 
its  time  dial  must  be  coordinated.   Table  7.2  lists  the  relay-curve 
coefficients. 

Due  to  the  inverse  time  shape  of  a  relay  curve,  many  points  on 
the  curve  must  be  checked  with  its  load-side  device  to  insure  that 
proper  coordination  exists.   The  manner  in  which  the  family  of  relay 
curves  are  represented  simplifies  the  problem  somewhat.   A  relay 
operating  time  of  O.A  seconds  is  assumed  throughout.   At  each  current 
value,  the  load-side-device  operating  time  is  summed  with  the  relay 
operating  time.   If  this  time  is  greater  than  the  line-device  time,  the 
difference  is  the  relay  time-delay  setting.   Thus,  when  all  currents 
up  to  the  fault  currents  are  checked,  all  points  on  the  line-device 
curve  are  at  least  0.4-seconds  higher  than  all  points  on  its  load-side 
device. 

Molded  Case  Circuit  Breaker 

The  molded-case  circuit  breaker  (MCCB)  with  its  different 
characteristics  provides  a  different  set  of  problems.   The  MCCB  curve 
cannot  be  represented  by  one  equation.   This  device  has  one  set  of 
curves  for  the  thermal  element  plus  another  for  its  magnetic  element. 
The  MCCB  used  in  the  program  is  a  General  Electric  type,  K-1200. 
Although  there  is  some  difference  in  the  curve  shape  for  the  current 
ratings,  300-1200  A,  they  are  similar  enough  to  be  represented  as  one 
type.   Each  MCCB  has  one  curve  corresponding  to  minimum  total  clearing 
time  and  a  maximum  total  clearing  time. 
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The  magnetic  portion  of  the  time-current  curve  cannot  be 
represented  adequately  by  a  polynomial  function,  but  is  represented 
by: 

If  Current  ^  Magnetic  Setting 

Minimum  Time  =  0,  Maximum  Time  =  .025  seconds. 

As  will  be  seen  later,  .025-seconds  response  time  is  quite  fast  and 
usually  does  not  present  a  problem  during  later  coordination. 

The  thermal  elements  of  molded-case  circuit  breakers  are 
susceptible  to  temperature  change  so  the  operating  temperature  must 
be  entered.   The  temperature  merely  shifts  the  curve  to  the  left  or 
right  so  it  does  not  affect  the  curve  shape.   This  offset  is  added 
or  subtracted  from  the  multiple  of  current  rating  when  using  the  thermal 
curve  polynomial  functions. 

The  problem  of  coordinating  the  thermal  elements  is  the  same  as 
that  for  fuses.   The  magnetic  element  coordination  is  much  simpler  in 
that  the  magnetic  element  is  set  separately.   Its  load-side  device 
operating  time  is  checked  at  a  current  value.   If  the  operating  time 
is  greater,  the  magnetic  element  setting  is  increased.   The  end  of  the 
curve  for  the  magnetic  element  is  determined  by  the  maximum  fault  cur- 
rent at  the  bus.   The  total  decoupling  of  the  overload  from  the  fault 
device  of  a  MCCB  makes  coordination  much  simpler  than  the  fuse  or  relay. 
The  more  or  less  proportional- time  shape  of  the  curve  will  ensure  that 
it  will  be  properly  coordinated  once  its  proper  current  rating  has  been 
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established.   Table  7.3  contains  the  polynomial  coefficients  used  in 
the  program. 

A  basic  flow  chart  (Figures  7.2  and  7.3)  follows  with  a  complete 
program  listing  in  Appendix  B. 

Summary 

In  this  chapter,  the  coordination  problem  was  solved  to  allow 
fuses,  molded-case  circuit  breakers,  and  relays  to  function  so  as  to 
clear  a  fault  in  the  area  of  the  fault  before  interrupting  the  rest  of 
the  distribution  system.   There  is  enough  flexibility  in  the  program 
to  allow  device  settings  to  be  changed  and  checked  using  the  program 
PL'OTD.FA.   An  example  problem  would  be  beneficial  at  this  point  to 
show  how  to  run  the  programs.   In  the  next  chapter,  a  sample  distribu- 
tion system  will  be  coordinated  using  the  coordination  programs. 
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Table  7.3.   Molded-case  circuit  breaker. 


log(time)  =   Z  a  [log(current)  -  (ofset  +  .1122)]'^  -  2 
n=0  " 


6 
log(time)  =   E  a . [log(current)  -  (ofset  +  .2315)^  -2 


j=o  J 


ofset  =  log(ckt  breaker  size) 


Minimum  Time : 


a  =    5.76908 

a  =  -10.88261 

a  =   23.65418 

a  =  -34.7922 

a,  =   33.7747 

a  =  -21.14536 

a^  =    6.43177J 
b 


Maximum  Time ; 


a  =   5.8596721 

a  =  -  4.3943349 

a  =  -   .50690283 

a  =  -10.402739 

a,  =   -17.0649561 

a  =  -13.011543 

a^  =  -  4.0107239 
6 


COORD. FA 
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Figure  7.2.   Coordination  flow  chart, 
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Figure  7.2.   Continued 
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Figure  7.3.   Device-graphing  flow  chart. 
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8.   OPERATING  THE  PROGRAMS 

General 

In  this  chapter  the  capabilities,  limitations,  and  use  of 
the  programs  will  be  described.   The  example  system  of  Figure  8.1 
will  be  coordinated  using  these  programs. 

One-Line  Diagram 

INPUT. F4.   All  data  is  entered  via  the  keyboard  in  the  format 
detailed  in  Chapter  3.   Data  entry  is  a  straightforward  process 
except  when  entering  transformer  data.   The  program  is  written  for  a 
"Y-Y"  configuration.   If  one  or  both  sides  are  "A"  configured,  then  a 
large  number  must  be  entered  for  the  transformer  neutral  impedance 
"(ZN)"  on  the  bus  connected  to  the  A  side.   After  the  program  is 
completedj  the  zero  sequence  value  for  the  transformer  element  con- 
taining the  transformer  may  be  assigned  an  arbitrary  high  value  if  the 
transformer  wiring  is  such  that  there  can  be  no  ground  current  through 
the  element.   This  is  done  by  changing  its  value  in  disk  file  TEST. DAT 
(see  Appendix  C  for  the  file  format) . 

When  entering  the  loads  the  program  asks  which  type  of  load 
to  be  entered.   The  various  loads  are  described  as: 
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Figure  8.1.   One-line  diagram  of  sample  system. 
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HP   -  Synchronous  motor 

HPl  -  Assynchronous  motor 

VA   -  In  KVA 

CUR  -  In  amperes 

S    -  No  more  loads 


Load  Flow 

INPUT. F4 .   There  are  no  entries  from  keyboard  for  this  program. 
Printing  DSK.DAT  will  provide  a  load-flow  analysis  in  per-unit  values. 

Solution  of  Fault  Currents 

012 
ZBUS . F4 .   In  addition  to  forming  Zbus    ,  this  program  allows 

the  input  for  mutual  coupling  and  for  zero-sequence  impedance  for  loads. 
Mutual  coupling  tends  to  increase  fault  currents  in  the  mutually  coupled 
elements  and  small  neutral-to-ground  impedances  can  create  higher  line- 
to-ground  fault  currents  than  three-phase  faults. 

FALT . FA .   Besides  the  solution  of  fault  currents  described  in 
Chapter  6,  FALT.F4  also  allows  the  input  of  fault  impedances.   This  can 
provide  an  interesting  study  of  the  system  because  the  user  can  go 
back  and  place  different  fault  impedances  into  the  fault  analysis  to 
see  if  the  protective  devices  will  operate  as  desired. 
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Device  Coordination  and  Plotting 

COORD. F4 .   The  program  forms  a  two-dimensional  complex  matrix 
called  ADEV  that  contains  all  the  information  required  to  plot  all  the 
devices  in  the  system  plus  all  actual  settings  that  are  to  be  made  on 
the  devices.   An  example  row  3  in  the  matrix  is  3,(3,4),  (-4,2), 
(1.26,203),  (1.03,1.26),  (3,667,0),  (.67,0),  (20,0)  in  which: 

1.  For  3,  3  is  the  row  location  of  the  matrix  and  corresponds 
to  the  element  number. 

2.  For  (3,4),  it  is  the  first  column  location  and  is  the  bus 
numbers  of  the  element. 

3.  For  (4, 2), -4  is  a  relay  type  #4;  2  is  a  current  limiting 
fuse  type  //2. 

4.  For  (1.26,203),  1.26  is  the  logarithm  of  the  pickup  point; 
203  contains  the  tap  information.   Two  is  the  tap  of  the 
CT,  and  03  is  the  tap  of  the  relay. 

5.  For  (1.03,1.26),  1.03  is  the  logarithm  of  the  current  for 
the  minimum  fuse  open;  1.26  is  the  typical  fuse  open  point. 

6.  For  (3.667,0),  3.667  is  the  magnetic  setting  of  the  relay; 
0  since  the  fuse  had  no  magnetic  setting. 

7.  For  (.67,0),  .67  is  the  time-dial  setting  for  the  relay. 
The  actual  time  dial  setting  must  be  read  from  the  manu- 
facturer's family  of  relay  curves  using  the  formula  below 
at  the  maximum  fault  current. 
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time-dial-setting  found  =  (found-time  curve)  - 
(time-dial  1/2  setting) 

8.   For  (20,0),  20  is  the  turns  ratio  of  the  current  trans- 
former. 

Device  numbers  the  programs  asks  for  correspond  to  the 
devices  below: 

1  -  Current  limiting  fuse 

2  -  Boric  fuse  high  voltage  only 

3  -  Molded-case  circuit  breaker 

4  -  Overcurrent  relay 

In  device  4  the  program  will  ask  later  what  particular  type 
of  relay  to  use  from  Table  8.1. 

Upon  completion  of  coordination,  a  device  listing  is  typed 
on  the  keyboard,  and  a  printout  of  the  same  is  available  by  typing: 
PRINT*. LPT. 

When  performing  the  plot,  having  a  device  listing  as  well  as  a 
fault  current  listing  is  helpful. 

PL0TD.F4.   This  program  plots  the  devices  that  were  coordinated 
above.   By  using  a  graph,  it  is  possible  to  visually  check  coordination. 
The  program  arrangement  allows  the  user  to  change  any  device  setting 
or  size  in  ADEV.DAT  to  tailor  the  coordination  to  his  particular  require- 
ments (see  Appendix  C  for  the  file  format) .   The  program  uses  the  lowest 
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Table  8.1.   Overcurrent-relay  types 


Program  Type  //  Relay  Type 


CO-2 


CO-3 


CO-4 


CO-5 


CO-6 


CO-7 


CO-11 


Relay  type  numbers  are  listed  as  negative  numbers  by  the  program. 
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voltage  in  the  system  as  a  base  for  plotting  currents.   The  user  can 
select  any  plot  offset.   This  has  the  effect  of  plotting  only  the 
area  of  concern.   In  this  way  a  very  detailed  graph  can  be  obtained. 

Tables  8.1-8.4  provide  a  listing  of  program  device  type  numbers 
to  the  actual  device  type  and  size. 

Table  8.5  details  what  is  contained  in  the  various  disk  loca- 
tions created  by  the  programs. 

To  coordinate  the  devices  in  Figure  8.1,  the  procedure  begins 
by  entering  element  data  using  INPUT. F4.   The  program  asks  if  data  is 
to  be  read  from  disk  or  from  the  keyboard.   The  first  element  must  be 
entered  using  the  keyboard. 

The  proper  procedure  is  to  group  the  elements  to  be  entered 
by  wire  or  cable  type  as  the  program  only  handles  one  type  of  wire  or 
cable  at  a  time.   When  completed  with  one  type, answer  zero  for  all 
other  questions  and  the  computer  writes  the  element  data  that  it  has 
computed  on  disk  in  the  file  RELA.DAT  and  stops.   The  new  wire  type 
polynomial  coefficients  are  copied  into  RES. DAT  and  the  program  is 
ready  to  be  run  for  the  new  wire  type. 

The  cable  coefficients  used  in  the  test  problem  are  for  mine 
power  cable  [5].   The  coefficients  are  listed  in  Table  8.6. 

After  selecting  keyboard, the  base  power  and  volts  are  entered. 
The  computer  then  informs  the  user  to  begin  element  data  entry.   The 
data  format  (Chapter  3)  is  listed  to  aid  the  user.   For  example,  the 
element  from  bus  one  to  two  of  Figure  8.1  is  entered  as  follows: 
1/2/3000/6.   For  the  element  containing  the  1  MVA  transformer  the  follow- 
ing is  input:   2/3/35/6/ . 001 . 005/1/15/4 . 1.   The  bus  base  voltages  are 


rH 
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Table  8.2.   Current-limiting  fuses. 


Program  Type  //  Limiting  Fuses 


1  5E 

2  7E 

3  lOE 

4  15E 

5  20E 

6  25E 

7  30E 

8  40E 

9  50E 

10  65E 

11  •    80E 

12  lOOE 

13  125E 

14  150E 

15  200E 
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Table  8.3.   Boric  fuses 


Program  Type  #  Size 


16  15E 

17  20E 

18  25E 

19  30E 

20  40E 

21  50E 

22  65E 

23  80E 

24  lOOE 

25  125E 

26  150E 

27  175E 

28  200E 

29  250E 

30  300E 

31  400E 

32  2-250E 

33  2-300E 

34  2-400E 
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Table  8. A.   Molded-case  circuit  breakers, 


Program  //  Size 


34  300  amps 

35  350  amps 

36  400  amps 

37  450  amps 

38  500  amps 

39  600  amps 

40  700  amps 

41  800  amps 

42  1000  amps 

43  1200  amps 


Table  8.5.   Disc  files. 
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DIGIT. F4:  Curve-smoothing  program 

DEV.F4:  Numonics  input  program 

INPUT. F4  Element  data-input  program 

LDFL0.F4:  Load-flow  program 

012 

ZBUS.F4:  Formation  of  ZBUS    program 

FALT.F4:  Solution  of  fault-current   program 

COORD. F4:  Coordinate  device   program 

PL0TD.F4:  Plot  device  curves 


RES. DAT: 
TEST . DAT 
RELA.DAT 
XFOR.DAT 
LOCUR.DAT 
DEV.DAT 
DSK.XFR 
DSK.DAT 
ZBUS . DAT : 
ZBUS. FOR: 
MUT . CO : 
FALT . DAT 
EFLT . DAT 
LCUR.DAT 
FOUR. DAT 
ADEV.DAT 


Contains  X/R  line  coefficients 

Contains  element  and  load  data  in  per  unit 

Contains  element  and  bus  voltage  data 

Contains  transformer  inrush  and  withstand  values 

Contains  load  currents 

Contains  relay  coefficients 

Contains  load-flow  currents  in  P.U. 

Contains  load-flow  analysis  data 

012 

Contains  ZBUS 

012 

Contains  ZBUS     before  loads  are  added 

Contains  mutual  coupling  data 

Contains  per-unit  fault  currents 

Contains  per-unit  fault  voltages 

Contains  real  values  of  load-flow  currents 

Contains  real  values  of  fault  currents 

Contains  coordinated-device  listing 
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Table   8.6.      Mine   power    cable   polynomial   coefficients. 


R-resistance 


X-reactance 


a      0.758690526734m-01 


0.6731452381A07621D-01 


a     -  .1108155100578006l>f02 


-  .1A08237729509665D-01 


a^  0.6290760668063740D+01 


-  .1721580927352306D-02 


a     -  .1627214281565102IH01 


6845126542245339D-03 


a,      0.1604746751660286l>f00 
4 


0.0 
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input  next  with  the  voltage  in  volts.   Upon  completion  of  entry  of  all 
bus  base  voltages,  the  computer  types  the  computed  element  value  in 
per  unit.   The  load  data  is  entered  next.   The  program  types  the  data 
entry  format  for  the  user.   The  50  hp  motor  is  entered  as:   8/50/,85. 
The  only  other  data  to  be  entered  are  mutual  coupling  and  zero  sequence 
load  impedance  in  ZBUS.F4  and  fault  impedance  in  FALT.F4.   The  sample 
problem  uses  no  coupling  and  a  fault  impedance  of  zero.   COORD. F4  is 
run  next.   The  device  type  desired  is  entered  when  asked  for  by  the 
program.    A  sample  data  input  is  given  in  Appendix  D. 


1 ' ^v- ; 
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The  following  is  a  list  of  load  locations  and 
currents . 

Bus  number   Load  current 

4      84.0396680 
4      84.0396  6  80 

7  293.6613400 

8  48.9435570 

The  following  is  the  current  flow  in  the  system 
under  full  load  conditions  as  found  from  the 
load  flov7  analysis. 


Bus 


to 


Bus 


Q 

Current  at  P 

Current  at  Q 

2 

34.222  6  440 

34 

.2226440 

3  ' 

34  .2226630 

125 

.2048700 

4 

88.0372570 

88 

.0372570 

5 

39.  1722710 

39 

.  1722710 

6 

39.  1721210 

365 

.0129500 

7 

313.4835700 

313 

.4835700 

8 

51.6573760 

51 

.6573760 

The  following  is  the  computed  fault  currents  for  a 
3  phase  fault  and  a  line  to  ground  fault. 


Bs  P 

Element 

3  Phase 

Ln  to  Gnd 

2 

1 

9380.7892000 

7053.929  8  000 

2 

2 

2  61.6855700 

13  3.7685400 

2 

3 

6  77.8190200 

595.4  521200 

2 

4 

280.0413400 

187.2686100 

2 

5 

2  8  0.3767000 

136.2585700 

2 

6 

2  2  26.35  4  7  000 

14  96.4606000 

2 

7 

3  8  3.7966800 

270.9381000 

3 

1 

124  7 . 294  1000 

833.362070  0 

3 

2 

12  4  7.2997000 

8  2  6.2558200 

3 

3 

84  6  .65  28  100 

6  9  9.4013200 
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Fit  Bs  P 

K 1  e  m  e  n  t 

3  Phase 

Ln  to  Gnd 

3 

4 

349.82658  00 

206.0  435  6  00 

3 

5 

349.79  9  1000 

204.2971400 

3 

6 

2  7  8  0.8159000 

1648.7  6  48000 

3 

7 

479.3  6  59400 

302.5682200 

4 

1 

1017.3596000 

469.6438400 

4 

2 

1017 .3546000 

466.3491300 

4 

3 

4003.0938000 

2210.2807000 

4 

4 

285.3238G00 

121  .9024600 

4 

5 

285.3332800 

121 .3072000 

4 

6 

2268.  1567000 

972.7699300 

4 

7 

391 .0159800 

174.0103800 

5 

1 

881 .8973500 

308.3295700 

5 

2 

8  81.889990  0 

3  0  7.2671700 

5 

3 

598.60  9  47  0  0 

260. 1660700 

5 

4 

3803.2506000 

13  8  4.1725000 

5 

5 

353. 1902700 

108.597070  0 

5 

6 

2  8  07.5742000 

901 .2339600 

5 

7 

483.9  912400 

18  3.8482000 

6 

1 

791 .2624200 

148.92  8  7400 

6 

2 

7  9  1.2  6  45000 

148.4773900 

6 

3 

537 .0803000 

125.4932500 

6 

4 

3416.9009000 

668.5896000 

6 

5 

3416.8793  00  0 

569.4760500 

6 

6 

2839.6040000 

517 .3090900 

6 

7 

489 .5244300 

124.5511200 

7 

1 

350. 4617  (5  00 

10  3.1131400 

7 

2 

350.4624900 

10  3.1008700 

7 

3 

237.8819000 

70.4041400 

7 

4 

1513.3939000 

4  4  5.7008900 

7 

5 

15  13.3895000 

443.9770900 

7 

6 

14273.8  6  30  0  00 

4235.3970000 

7 

7 

216.8161500 

64.2673200 

8 

1 

156.6019600 

23.5717640 

8 

2 

156.6021300 

23.5401980 

8 

3 

106.3009100 

16.4766590 

8 

4 

676.2499600 

102.2285800 

8 

.5 

676.2485300 

100.  1874100 

8 

6 

561.99519  0  0 

8  4.0780670 

8 

7 

6769.2755000 

1059.3405000 
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Using  the  device  listing  of  ADEV.DAT  and  the  plotted  curves, 
it  is  possible  to  obtain  all  device  settings. 

Device  number  one  is  a  molded-case  circuit  breaker.   Its  size 
is  300  A  and  its  magnetic  trip  is  on  minimum  setting. 

Device  number  10  is  a  current-limiting  fuse  number  10. 

Device  number  12  is  a  solid-material  high-voltage  boric  fuse 
number  24. 

The  relay  data  follows: 


Device 

Type 

Time 

Number 

Number 

Ct 

Tap 

Dial 

Magnetic 

3 

1 

3 

(120) 

7 

1/2 



4 

1 

3 

8 

8 

6 

5 

3 

3 

14 

5 

6 

6 

4 

3 

3 

4 

6 

7 

5 

3 

4 

6 

6 

8 

5 

3 

5 

6 

6 

9 

6 

3 

6 

10 

5.2 

2 

1 

3 

13 

5 

5.64 

11 

1 

3 

4 

10 

6 

13 

7 

3 

7 

2 

6 

14 

2 

3 

8 

1 

11 

6 

The  magnetic  trip  of  device  number  9  had  to  be  reduced  after 
coordination  to  allow  the  device  12  to  set  properly  to  provide  pro- 
tection for  the  1  MVA  transformer.   The  graphs  that  follow  show  the 
devices  plotted  to  the  lowest  voltage  in  the  system,  440  V.   All  the 
graphs  were  computer  generated.   The  first  graph  shows  how  a  one-line 
computer  diagram  would  be  placed  to  aid  in  understanding  the  graphs. 
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Lgure  8.2.   Example  coordination  plot 
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gure   8.2.      Continued, 
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gure   8.2.       Continued, 
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Summary 

The  sample  distribution  system  of  Figure  8.1  was  coordinated 
using  the  coordination  programs.   Using  the  device-plotting  program 
PL0TD.F4  this  coordination  was  checked  to  see  if  the  computer  generated 
coordination  was  satisfactory  to  the  user. 

During  the  research  for  this  thesis,  many  observations  were  made 
and  many  features  could  be  added  to  these  programs  in  the  future.   This 
will  be  discussed  in  the  chapter  that  follows. 
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9.   CONCLUSIONS 

In  this  thesis  the  author  set  out  to  represent  a  one-line 
distribution  system  in  a  per-unit  system  and  to  solve  for  three- 
phase  and  line-to-ground  faults  and  to  perform  a  load-flow  analysis. 
The  various  protective  devices-response  curves  needed  to  be  represented 
by  polynomial  functions.   These  requirements  were  fulfilled  using  the 
methods  detailed  in  Chapters  2-7. 

The  last  task  was  to  develop  a  computer  procedure  to  coor- 
dinate a  distribution  system's  protective  devices.   Fuses,  molded- 
case  circuit  breakers,  and  relays  were  coordinated  and  their  coordi- 
nation graphs  were  plotted.   There  are  many  extensions  for  further 
work  that  could  be  done.   These  extensions  will  be  discussed  later. 

While  working  with  the  various  programs,  several  observations 
have  proved  useful  when  performing  coordination.   Coordination  should 
not  be  an  afterthought  in  designing  a  distribution  system.   It  is 
easy  to  design  a  system  that  is  impossible  to  coordinate  correctly. 
A  radial  system  is  easier  to  coordinate.   Moving  from  a  load  toward 
the  source,  the  time-current  curves  shift  to  higher  current-time 
responses  so  that  when  finally  arriving  at  a  transformer,  providing 
proper  transformer  protection  and  coordination  simultaneously  may  be 
impossible . 

Coordination  of  the  thermal-tripping  elements  of  molded  case 
circuit  breakers  is  very  difficult  due  to  the  difference  between  the 
minimum  and  maximum  tripping  time  curves.   Coordination  of  these 
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devices  is  best  done  manually  with  some  overlap  of  the  thermal  time 
curves.   The  magnetic  elements  could  then  be  adjusted  to  provide 
proper  fault  coordination.   Those  manually  adjusted  devices  can  be 
added  to  ADEV.DAT  and  be  checked  using  PL0T.F4. 

In  systems  with  different  voltages,  the  transformers  should  be 
of  different  sizes  so  that  the  smaller  sizes  are  on  the  load  side. 
This  will  allow  proper  coordination  and  protection.   Relays  are 
preferable  to  other  types  of  protection  on  the  high-voltage  parts 
of  the  system  in  that  their  response  can  be  controlled. 

Recommendations  for  Further  Work 

1.  Using  an  analog  computer  and  the  programs  presented  here, 
an  entire  distribution  and  generation  system  could  be 
looked  at  in  detail.   A  transient-stability  study  would 
show  the  effect  of  various  device  settings. 

2.  The  INPUT. FA  program  could  be  enlarged  to  automatically 
account  for  different  transformer  configurations. 

3.  Additional  devices  could  be  added  easily  to  the  program. 
The  program  allows  up  to  twenty  different  relay  types.   The 
additional  coefficient  data  is  merely  entered  in  disk 
DEV.DAT.   Additional  molded-case  circuit  breakers  would  be 
slightly  more  difficult  in  that  it  would  require  some  program 
modification. 
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4.  The  plotting  program  could  be  enlarged  to  include  drawing 
a  one-line  diagram  in  the  space  provided  and  adding  other 
data  directly  on  the  graph. 

5.  The  program  could  be  expanded  for  the  case  of  having  one 
device  in  an  element. 
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C         T1TLE:DIGIT.F4 ;CREATES  FILES  FROM  DATA 
USING  NUMONICS  ARM 

1  TYPE  5 

5         FORMAT( '  INPUT  DEVICE  NR..') 

ACCEPT  8,IDEV 
8         FORMAT(I) 

OPEN(UNIT=IDEV) 

TYPE  10 
10        FORMAT( '  INPUT  X-AXIS  DIMENS ION , Y- AXI S 
DIMENSION'//'  .  .  ,  '  ) 

ACCEPT  15,XAX,YAX 
15        F0RMAT(2F) 

TYPE  20 
20        FORMAT('  INPUT  LENGTH  PER  LOG  1 0 ; X , Y .  .  '  ) 
30       ACCEPT  35,XL,YL 
35        F0RMAT(2F) 

XAX=100*XAX/XL 

YAX=100*YAX/YL 

CALL  OPEN (30, 'DEVICE', 'NUM' ,' MODE ',' IMAGE ' , 

l'ACCESS','SEQ IN', 'BUFFERS', 4) 

CALL  NUMRES(30,4) 

CALL  LOCK( 1,1) 
100       READ(30,END=200)  IDATA 

IX=( I  DATA/ "10000 00) .AND. "37 777 

IY=IDATA. AND. "37777 

IFLAGS=( IDATA. AND. "7 40000)/ "4 0000 

X=IX/XAX 

Y=IY/YAX 

WRITE ( IDEV,150)IFLAGS,X,Y 
150      F0RMAT(I,2F) 

GO  TO  100 
200       CALL  NUMCLO 

CALL  UNLOCK( 1,1) 

CLOSE(UNIT=IDEV) 

TYPE  210 

2  10       FORMATC'  CONT I NUE ?..,'$ ) 

ACCEPT  220, IC 
220       FORMATC I) 

IF( 1C.EQ.O)GO  TO  1 

STOP 

END 


C         CURVE  FITTING  PR  OCR  All  :  DEV  .  FA 

**************************  >'cAA;!cA*AA*>V***;k*A*yt  ****:<*** 

INTEGER  L(7 ) ,MM( 7 ) 

REAL  X( 100) ,YS( 100) 

REAL*  8  YM(  100),G(7,7),I](7),A(0/6),GG(7,7),  ERROR  (  100) 

TYPE  100 
100       FORMATC IX, ' CHOOSE  DIMENSION  OF  FIT   '$) 

ACCEPT  110, ND 
110       FORMAT(I) 

TYPE  120 
120       FORMATC  IX,  '  CHOOSE  NO.  OF  POINTS     '$) 

ACCEPT  110,M 

TYPE  130 
130       FORMATC /' CHOOSE  INPUT  FILES'/'   ••.,     '$) 

ACCEPT  110,IUNIT1 

DO  lAO   I=1,M 

READCIUNIT1,135)J,XC1),YSCI) 
135       F0RMATCI,2F) 

YSC I)=ALOG10CYSC I) ) 
140       CONTINUE 

DO  10  K=1,M 

DO  10  1=1 ,ND 

UCI)=UCI)+YSCK)*XCK)**CI-1) 

DO  10  J=l ,ND 

GCl,J)  =  GCl,J)  +  XCK)^'-^Cl  +  J-2) 
10        CONTINUE 

CALL  DMINVCG, 7 ,ND) 

DO  20  I=0,ND-1 

DO  20  J=1,ND 

ACI)  =  AC  I)  +  GCI  +  1  ,  J)'''UC  J) 
20        CONTINUE 

DO  30  1  =  1, M 

DO  30  J=0,ND-1 

YMCI)=YMCI)+AC J)*XCI)**J 
30        CONTINUE 

DO  40  1=1 ,M 
40        SQUERR=SQUERR+CYSC I)-YMC I)  )**2 

TYPE  160,SQUERR 
160       FORMATC'  THE  INTERGRAL-SQU ARED  ERROR  IS:     ',E) 

0PENCUNIT=4  0, FILE=' A.DAT' ) 

WRITEC40, 170)A 

CLOSECUNIT=40) 
170       FORMATCD) 

TYPE  99,IUNIT1 
99        FORMATC'  1UNIT=',I) 

DO  5  0  1  =  1  ,M 
50        ERRORCl)=DABSCCYMCl)-YSCl))/YSCl)) 

0PENCUN1T=4  1,FILE='YIN0UT.DAT') 

WRITEC41  ,  1H0)CXC 1)  ,YSC I)  ,  YMC  I)  .ERRORCI)  ,  1  =  1  ,  M) 
180       FORMATC 2F, 2D) 


c 

MATRIX 
C 


10 
15 


20 


25 


30 
35 

38 


40 
45 

50 
55 


60 
6  2 
65 


70 
75 


CLOSE(  lJ^JIT  =  4  1  ) 
STOP 
END 

TITLE :UMINV; SUBROUTINE  TO  INVERT  A  REAL  SQUARE 
DOUBLE 
PRECISION 

SUBROUTINE  DMINV( A, MDIM, N) 
REAL*8  A(MD1M,N) 
INTEGER  L( 100) ,M( 100) 

CONVERTED  FROM  SSP  ROUTNE  MINV 


DOUBLE  PRECISION  BIGA,HOLD 

DO  8  0  K=l ,N 

L(K)=K 

M(K)=K 

BIGA=A(K, K) 

DO  20  J=K,N 

DO  20  I=K,N 

IF(DABS(BIGA)-DABS(A(I,J))) 

BIGA=A(I,J) 

L(K)=I 

M(K)=J 

CONTINUE 

J=L(K) 

IF(J-:0  3  5,35,25 

DO  30  1  =  1  ,  N 

H0LD=-A(K, I) 

A(K,I)=A( J,I) 

A( J,I)=HOLD 

I=M(K) 

IF(I-IC)  4  5,4  5,38 

DO  40  J=1,N 

H0LD=-A( J,K) 

A( J,K)=A( J,I) 

A( J, I)=HOLD 

DO  55  1=1  ,N 

IF(I-K.)  50,55,50 

A( I, K)=A( I,K)/(-BIGA) 

CONTINUE 

DO  6  5  1  =  1  ,  N 

HOLD=A( I , K) 

DO  65  J=1,N 

IF(I-K)  60,65,60 

IF(J-K)  62,65,62 

A(I,J)=H0LD*A(K,J)+A(I,J) 

CONTINUE 

DO  7  5  J=l  ,N 

IF(J-K)  70,75,70 

A( K, J)=A( K, J)/BIGA 

CONTINUE 


15.20,20 
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A(K,K)=(1.0000000)/BICA 
80        COMTINUE 

K=  N 
100       K=K-1 

IF(K)  150,15  0,105 
10  5       I  =  L ( K ) 

IF(I-K)  120,120,108 
108      DO  110  J=l  ,N 

H0LD=A( J,K) 

A( J,K)=-A( J,I) 
110       A(J,I)=H0LD 
120       J=M(K) 

IF(J-K)  100,100,125 
125       DO  130  1  =  1, N 

H0LD=A(K, I) 

A(K,I)=-A( J,I) 
130       A(J,I)=H0LD 

GO  TO  100 
150       RETURN 

END 
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APPENDIX  B 


COORDINATION  AND  PLOTTING  PROGRAMS 
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5'.- *A*A**/^A*******A****A***5k***-;c  *******  A  ******  Vc  ********** 

C        INPUT.  IT  INPUTS  DATA  FPvOM  A  TERM  I NAL ,  PLACES  IN 
C         THE  PROPER  FORMAT  FOR  OTHER  PROGRAMS 
******************************************************* 

C 

c 

DIMENSION 
L(30,30,2),PU(30,30,1),RES(30,3  0,2),E(30),XL(3  0,30,2) 

DIMENSION 
S(300),ZER(30,3O,2),RC(7),XC(7),ZG(30),XD(3  0) 

COMPLEX  T1,TPU,PU,S,S1,PU1,PZ,CUR,ZER,ZG 

REAL*8  A(7),AA(7) 

REAL  L 
********************************************************** 

C         IF  THE  PROGRAM  HAS  ALREADY  BEEN  RUN  ONCE  ALL  DATA 
C         NEED  NOT  BE  ENTERED  AGAIN 
*********************************************************** 

TYPE  3 
3         FORMATC  '  DO  YOU  VJANT  TO  READ  DATA  FROM  A  FILE  OR 
FROM  KYBRD?'/ 

2'  <CR>FOR  INPUT  FROM  TTY , 1  FOR  READ...') 

ACCEPT  11,1 

IF(I.NE.0)GO  TO  AGO 

*************************  v'c  *************************  ********* 

C         INPUT  OF  DATA  ON  LINES 

vV*************  ***********************************  ******* 

TYPE  5 

5  FORMATC 2X, ' ENTER  THE  PROBLEM  BASE  IN  MVA') 
ACCEPT  lO.ZBASE 

TYPE  6 

6  FORMATC 2X, ' ENTER  THE  PROBLEM  BASE  IN  KV') 
ACCEPT  IG.EBASE 

10  FORMATC F) 

11  FORMATC I) 
GO  TO  14 

********************************************************* 

C         PROGRAM  RETURNS  HERE  AFTER  READING  DATA  AT  AGO 
******************************************************** 

12  TYPE  13 

13  FORMATC'  OLD  DATA  IS  ENTERED , RE ADY  FOR  NEW  LINE 
DATA.  .  '  ) 

ACCEPT  11,1 

14  TYPE  15 

15  FORMATC /, 2X, ' ENTER  ELEMENT  DATA  AS 
FOLLOWS  :BUSC ) TO C ) C LENGTH) 

1CSP)SIZE'  ,/,2X,'  C TRANSFORMER 
IMPEDENCEC  P.U. ) ) C TRANSFORMER  SIZE 

2IN  MVA)  SP',/,  2X, ' CTRANSFORMER  HIGH  VOLTAGE  IN  KV) 
SP',/,' 

3C TRANSFORMER  LOW  VOLTAGE  IN  KV ]...,'/ ) 
******************************************************** 
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C 
C 

PROGRAM: 


20 


21 


22 
24 


26 
27 


31 
32 
33 
34 


THESE  TWO  FILES  CONTAIN  COEFFICIENT  DATA  FO  40  BELOW 
THEY  ARE  CREATED  USING  THE  CURVE  FITTING 
CURV . F4 

*************  5'c**A**Aic*yc5!c**5!c;k- ******  ■>:t*****5^*A*5lc***A 

OPEN(UNIT=40,FILE='RES.DAT') 

OPEN(UNIT=3  5, FILE='XFOR.DAT' ) 

READ(40,20)(A(I) , 1  =  1  ,7  ) 

READ(40,20)(AA(I),I=1,7) 

FORMAT(D) 

DO  21  1=1  ,  7 

RC(I)=A(I) 

XC( I)=AA( I) 

CONTINUE 

CLOSE(UNIT=40) 

CL0SE(UNIT  =  4  1) 

ACCEPT  24  ,1  ,  J,XLL,  WIR,I  ]  ,  ][.]   ,TVH,TVL 

F0RMAT(2I, 7F) 

IF(I.EQ.O)  GO  TO  110 

TYPE  26 

JJ=1 

FORMAT('    PARALLEL  LINE     1  TO  3 

ACCEPT  27, JJ 

FORMAT( I) 

IF( JJ.EQ.O) JJ=1 

L(  I  ,  J, JJ)=XLL 

IF( WIR.EQ. 140)G0  TO  31 

IF( WIR.EQ. 130)G0  TO  32 

IF( WIR.EQ. 120)G0  TO  33 

IF( WIR.EQ. 1 10)G0  TO  34 

SA=WIR 

SB=WIR 

IF( SA.LT 

GO  TO  40 

SA=-3 

GO    TO    35 

SA  =  -2 

GO    TO    33 

SA  =  -1 

GO    TO    35 

SA  =  0 

GO    TO    35 


,$) 


37)GO    TO     35 


******************  ******yc  **************************** 
C  CON' 


******** 
35 


VERSI 
******** 


ON  FROM  AWG  TO  MCM 
************************************* 


SC  = 
SD  = 
SE  = 
SB  = 
*********** 

C         F  (J  K 
*********** 


9*''<(36-SA) 
SC 

SE 


1  .  122 
.005* 
SD**2 
1  0  0  0  * 
********************************************* 


MULAE  FOR  LINE  RESISTANCE  AND  REACTArJCE 
******************************************** 
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40 

RESC  I  ,.I,  JJ.)  =  RC(  1  )  +  RC(2)'>  (  ALOGIOC  SB)  )  +  RC(  3)''(  AL0G1()(  SB 
)**2)+ 

lRC(4)*(ALOniO(SB)**3)  +  RC(  5)  "(  ALOG]  n(  SB)'''''^4  )  +  RC(6)'<(  A 
LOCIO(SB) 

2'''"5)  +  RC(7  )>^(  ALOGin(  SB)--'^'U)) 

XL(I,  J,  JJ)=XC(  1  )  +  XC(  2)'>(  ALOGIOC  S  B  )  )  +  XC  (  3  ) ''^  (  ALOG  1  0  (  S  B  ) 
**2)  + 

1  XC(  A  )  ''^  (  ALOGl  0  (  S  B)  *  ^3  )  +  XC(  5  )  "  (  ALOGl  0  (  SB)  •-^''■-A  )  +  RC(  6  )  '•  (  A 
LOGIO(SB) 

2  *vc  5  )  +  xC  (  7  )  '•  (  ALOG  1  0  (  S  B  )  '"^  '"^  6  ) 

TYPE  5  0 , RE  S ( 1 , J , J J )  , S  A ,  S  B , XL ( I , J , J  J ) 
50        F0RMAT(2X, ' LIME  RES.  ',3F,/,'LINE  REACT.. ',F) 

IF(CABS(T1 ) .EQ.0)GO  TO  85 

PU(I,J,JJ)  =  Tl*(ZBASE/TSl)*((TVH/EBASE)*^-2) 

TYPE  55  ,  I,  J 
5  5        FORMAT( '  INPUT  TRANSFORMER  ZN  IN  P.U.  REFER  TO  BUS 
SIDE, ( BUS 

1'  ,13,'  )SP(BUS  '  ,  13,  '  )  .  .  ,  '  ) 

ACCEPT  60, SI , PUl 
60        F0RMAT(4F) 

TV  =  TV[1'^  1E3 

ZG(I)=S1*(ZBASE/TS1)*((EBASE/TVH)**2) 

ZG(J)  =  PUl*(ZBASE/TSl)^((EBASE/TVL)''-"2) 

CURB=TS1*1E6/TV 

TINRSH=12*CURB 

TWST1=25*CURB 

TUST2=1A.4*CURB 

WRITE(35, 70)1, J,TINRSH,TWST1 ,TUST2 ,TS1 
70       FORMAT( 21, 4F) 
85        IF( IMAX.GT.I)  GO  TO  90 

IMAX=I 
90        IF( JMAX.GT. J)  GO  TO  100 

JMAX=J 
100       ZER( I, J, JJ)=CMPLX(0 .0,XL(I , J, JJ) ) 

TYPE  102 ,ZER( I , J, JJ) , SSI 
102       F0RMAT(3F) 

GO  TO  22 
110       0PEN(UNIT=31,FILE='RELA.DAT') 

0PEN(UHIT=32,FILE='TEST.DAT') 

0PEN( UNI T=3 3 , FILE=' LOCUR.DAT' ) 

IF( INAX.LT. JMAX)  GO  TO  115 

IB=IMAX 

GO  TO  118 
115       IB=JMAX 

C         BUS  VOLTAGES  ARE  INPUT 

>V  *********  :*.-*****************  >V  ***********  V>5^- *************** 
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1  18 
1  19 

120 
121 
122 

125 

126 

128 


DO 
TY 
FO 
AC 
IF 
E( 
CO 
TY 
FO 
TY 
FO 
WR 
FO 
El 
l.'R 
FO 
DO 
DO 
DO 
IF 
NN 
PZ 
PZ 
PZ 
Tl 
ZE 


TY 
130       FO 

11 
IMPEDENCE' 


1 

PE 
RM 
CE 
(E 
I) 
NT 
PE 
RH 
PE 
RM- 
IT 
RM 
=  E 
IT 
RM 
1 
1 
1 
(L 
=  N 

=  ( 
Z  = 
=  P 
=  P 

R( 

PE 
RH 
X, 
,2 


20  1  = 
119, 
AT(  ' 
FT  10 
2  .EQ. 
=  E2 
INUE 
122 
AT(/, 
125, 
■AT(  2X 
E(31, 
AT(  2F 

(l)/( 
E  (  3  2  , 
AT(I, 
35  1  = 
33  J  = 
33  JJ 

(I.J, 

N+1 
CMPLX 
(ZBAS 
Z*PZZ 
Z  +  PU( 
I, J, J 
130, 
AT(  2X 
'  TRAN 
F) 


1  ,  IB 

I 

INPUT  BUS' , 13, 'VOLTAGE. 

,E2 

0)00  TO  121 


..,',$) 


2X,'DATA  IS  ACCEPTED  AND  A  LIST  FOLLOWS') 

ZBASE 

,' PROBLEM  BASE  IN  MVA..,',F) 

126)ZBASE,EBASE, IB 

,1) 

EBASE*1E3) 

128) , IB, El 

2F) 

1  ,30 

1  ,30 

=  1,2 

JJ).EQ.0)GO  TO  133 

(RES(I,J,J.J),XL(I,J,JJ)))'^(.001*L(I,J,JJ)) 
E*lE6)/(E(I)'^--^2) 

I , J, JJ) 

J)=ZER(I,J,JJ)*PZZ*(.001*L(I,J,JJ))+T1 

N  r^  L  (  I  ,  J  ,  J  J  )  ,  T  1 

, ' ELEMENT   '  ,I3,5X,  'LENGTH'  ,F, /  , 

SFORMERPER  UNIT  IMPEDENCE  PLUS  LINE 


WRITE(31,131)I,J,JJ,L(I,J,JJ),RES(I,J,JJ),PU(I,J,JJ) 
,J,JJ) 


XL(I 
131 


1 


,  ZG(I) 


132 
133 
135 


136 


137 
138 

139 

lAO 
1A5 


1,ZGU;  ,ZG(.J) 
F0RMAT(3I , 9F) 
WRITE(32, 132)NN 
F0RMAT(4I , 8F) 
CONTINUE 
CONTINUE 
VJRITE(  31,136) 
36) 
0  0 
IB 


I,J,JJ,T1,ZER(I,J,JJ),ZG(I),ZG(J) 


1  =  1, 

137)1 
2F) 


WRITE ( 3  2 
FORMATC ' 
DO  138  I 
WRITE(3  1 
FORMATC I 
CONTINUE 
WRITE(31. 
FORMAT('  0  0.0 
TYPE  14  0,13 
FORMATC lOX, 
TYPE  150 


0 


0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0') 

ECD 


139) 


) 


0.0' 
Iv^UMBER  OF  BUSSES  IS, 


1) 
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*****A**i!c5k***v'cAjic*5t*****A**AA***5kA***A****  **************** 

150 
MANNER; 


FORMAT(//'  NOW  ENTER  THE  LOAD  DATA  IN  THE  FOLLOVJING 


160 

? . .VA,H 

162 

164 


165 


168 


170 
175 


178 


180 
185 

187 


1 ' / ' (BUS   ) (LOAD 
2'/'  (BUS  )  (LOAD 
3'(BUS  )(LOAD  IN 
DO  300  N=1,IB 

TYPE  160 
FORMAT( ' 
P,OR  AMPS?  ? '  ) 

ACCEPT  16 

FORMAT( A4 

TYPE  164, 

FORMAT( 2X 

IF(MD.EQ. 

IF(MD.EQ. 

IF(MD.EQ. 

IF(MD.EQ. 

MD=' VA' 

ACCEPT  16 

FORMAT( I , 

SZTZ=AS1/ 

SZR=SZTZ* 

PCOS=ACOS 

SZI=SZTZ* 

S( J)=CMPL 

AI=AS1*1E 

WRITE(33, 

FORMATd  , 

ST  =  MD 

GO  TO  300 

ACCEPT  17 

FORMAT( I , 

SZTZ=PRH* 

SZR=SZTZ* 

PCOS=ACOS 

SZI=SZTZ* 

S(I)=CMPL 

IF(MD.EQ. 

AI=SZTZ/( 

WRITE(33 , 

FORMAT( I , 

ST  =  MD 

GO  TO  300 

ACCEPT  18 

FORMAT( I , 

WRITE(33, 

FORMAT( I , 

CURR=CURT 

PCOS=ACOS 

CURI=CURT 


IN  KVA)(PWR  FACTOR) .. .OR. 
IN  HP)(PWR  FACTOR) ... .OR 
AMPS) (PWR  FACTOR) ' ) 


VniAT  METHOD  WILL  YOU  USE 


2  ,MD 
) 

MD 

,A4,' 
'  HPl' 
'HP'  ) 
'  AMPS 
'  S'  ) 


..'/ 


ENTER. . . ' ) 
)  GO  TO  170 

GO  TO  170 
')  GO  TO  180 
GO  TO  305 


5, 

2F) 

(ZB 

PFA 

(PF 

SIN 

X(S 

3/( 

168 

F) 


J, ASl , PFAC 

ASE*1E3) 

C 

AC) 

(PCOS) 

ZR, SZI) 

E( J)*SQRT(3 .0)) 

)J,AI 


5,1, PRH, PFAC 

3F) 

746 

PFAC 

(PFA 

SIN( 

X(SZR,SZI)/(ZBASE*1E6) 

'HP' 

SQRT 

178) 

F) 


,C) 
PCOS) 


)XD(I)=. 17/CABS(S(I)) 

(3.0)''>-E(I)) 

I,  AI 


5,1, CURT, PFAC 

3F) 

187)1 , CURT 

F) 

*PFAC 

(PFAC) 

*SIN( PCOS) 
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CUR=CMPLX(CUKR,CURI) 

CUR=CUK*(SQRT(3.0)*EBASE/ZBASE) 

S(I)  =  3'''EBASE*CONJG(CUR)/SQRT(3.0) 

ST  =  MD 
300      CONTINUE 
305       DO  320  N=1,IB 

NN  =  N 

TYPE  3  10,N, S(N) ,MD 
310       F0RMAT(  '  BUS  '  , I  3 , 2F , A4 , 2F ) 

WRITE(31 , 31 2)NN, S(NN) 
312       F0RMAT(I,2F) 

WRITE(32 , 315)N, S(N) ,XD(NN) 
315       F0RMAT(I,3F) 
3  20       CONTINUE 

WRITE(3  1 ,330) 

WRITE ( 32,331) 

3  30      FORMATC'   0  0.0  0.0') 

331  F0RHAT('   0  0.0  0.0') 
WRITE(33 , 332) 

332  FORMAT( '  0  0.0') 
WRITE(35, 335) 

3  35      FORMATC'  0  0  0.0  0.0  0.0  0.0  ') 

GO  TO  900 

4  00       0PEN(UNIT=3  6, FILE=' RELA.DAT' ) 

READ(36,411)  Z BASE , EB ASE , I B 

IMAX=IB 
All       F0RMAT(2F,I) 

4  50      READ(3  6 ,45  6)I,J,JJ,XL3,RES1,PU1,XL2,T1,S1 
4  56      F0RMAT(3I, 9F) 

IF(I.EQ.O)GO  TO  4  70 

L( I, J, JJ)=XL3 

RES( I,J, JJ)=RES1 

PU( I, J, 1)=PU1 

XL( I, J, JJ)=XL2 

ZER(I,J,JJ)=CMPLX(0.0,XL2) 

ZG( I)=T1 

ZG( J)=S1 

IF(I.NE.0)GO  TO  450 
470      READ(3G,472)I,E(I) 
472      F0RMAT(1,2F) 

IF( I.NE.O)  GO  TO  470 

457  READ(36,458) I, SI 

458  FORMATC I, 2F) 

sci)=si 

IFCI.NE.O)  GO  TO  457 
REWIND  3  6 
JJ  =  0 

CL0SECUIUT  =  36) 
GO  TO  12 
9  00       CLOSEC UNIT=3  1  ) 
CLOSECUNIT=32) 
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CLOSE(UNIT=35) 
STOP 

Ei:u 

:fc  *  A  *  >V  *  *  *  ;k  Vc  A  A  v%  A  *  *  *  ii  *  ^'c  :fc  A  >V  A  ^-  ;'c  A  A  vk  :<  ;k  X  jSc  >k  A  ik  3^  *  A  *  *  ;!c  -^  A  Vc  A  5c  ;^  *  ^s  :fc  >V  5'c  ;Ic  :i  >^  :fc  ■ 

C         LDFLO.FA:A  LOAD  FLOW  PROGRAM 

AAAAAAAAAA*AAAA*A*A*AAAAAAAAA*>VA*AA*AAAAAAAAAAAs'c*AAAAAAAA 

DIMENSION 
E(33),Y1(33,33),PSCH(33),QSC}I(33),PCAL(33), 

1QCAL(33),PDIFF(33),0DIFF(33),IC(33),S(33),JAC(68,6  8) 

2PQVEC(33),EFDV(6A),QQVEC(33), 

3YS(33,33),Y(33,33),CRFLO(33,33),U(68),V(68), 

4EDIFF(33) 

INTEGER  P,Q,M,U,V,PP 

COMPLEX  E,IC,Y,Z,YS,Y1,PFL0,PL0S,YSH,EDIFF, 

IS, CRFLO, PFL01,ZG1,ZNI,ZNJ 

EQUIVALENCE  (CRFLO, Yl) 

REAL  JAC 

DIMENSION  EIDV(64) 

OPEN(UNIT=30,DEVICE='DSK' ,FIL£=' TEST. DAT' ,DISPOSE='SA 
VE'  ) 

OPEN(UNIT=3  3, DEVICE='DSK' , FILE=' DSK.XFR' ) 

CREAD  VALUES  OF  IMPEDANCES   YSHUNT  ADMITTANCES 


C         888888*****IF  SHUNT  ADMITTANCES  ARENOT  IGNORED 

REMOVE 

*A*'j'}^£  Q    XN  THE  LINE  YS(0  p)*  *'><******** '^  ************** 

READ(3n,2)  IB, El 
2         FORMAT(I,F) 

E( 1 )=CMPLX(E1 ,0.0) 

DO  10  1  =  1  ,  1000 

READ(30, 5)NN,P,Q,JJ,Z,ZG1,ZNI,ZNJ 

5  FORMATC AI , 8F) 
IF(P.EQ.O)  GO  TO  15 
YS(P,Q)=YSH 

C         YS(Q,P)=YSH 

IF(CABS(Y1(P,Q)).NE.0)Z=1.0/(Y1(P,Q)+1/Z) 

Yl (P,Q)=1 .O/Z 

Y1(Q,P)=Y1(P,Q) 

PRINT  6,P,Q,Y1(P,Q) ,YS(P,Q) 

6  FOKMAT( 21 , AF) 
10        CONTINUE 

**********************************  A AAAAAAAAAAA A *A*AAA***AA 

CCOMPILE  TERMS  FOR  YBUS 

(^A*A*AA*AAAAAA***AA*A****AAAAAAAAAAA*AAA*AAA***AAAAAA* 

1 5        DO  5  0  P=l , IB 
Y( P, P)=0 
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DO     4  0     0=1  ,  Hi 
IF(P.EQ.O)     GO     TO     4  0 
Y(P,P)=y(P,P)+Yl(P,Q)+YS(P,0) 
Y(P,Q)=-Y1(P,0) 
40  CONTINUE 

50  CONTINUE 

C  PRINT  OUT  VALUES  OF  YBUS 

PRINT  5  1 

51  FORMAT  ( '0' , lOX, ' BUS  IMPEDANCE  MATRIX , YBU S ' ) 
DO  58  1=1 , IB 

PRINT  54, (Y(I, J) , J=l , IB) 
54        FORMATC '0' , lOOF) 
58        CONTINUE 

Qii  -k  -k  -k  ^  :k  'k  i(  it  i(  -k  -k  i<  -k  :i't  y\  i<  'J:  i<  i<  -k  -k  :k  i<  -k  it  'k  -k  -k  -k  'k  -k  i(  i<  :k  -k  ic  -k  :k  :k  •/<  :k  :k  -k  ■/(  -k  -k  -k  ^  -k 

C  READ  VALUES  OF  POWER  AT  BUSSES 

Qkkkkkkk-kkk-kkkk-kkk-k-kk-kk-kik-k-kkit-k-k-k'kk-k-k-kkkkk-kkkkkkkkkk 

DO  60  P=l  ,  IB 
READ(30,59)  I,PS,QS 
IF( I.EQ.0)GO  TO  62 
IF(I.EQ. 1 )G0  TO  60 
PSCH(I)=-PS 
QSCH( I)=-QS 
FORMATC I, 2F) 
CONTINUE 

Q-k  k  k  k  -k  k  k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  -k  k  k 

C  ASSUME  VALUES  OF  VOLTAGE  AT  BUSSES 

Qkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

62        DO  70  P=2,IB 

E(P)  =  (1  .0,0.0) 

CONTINUE 

K=0 

DO  80  1=1 , IB 

TYPE  75, I ,E(I) 
FORMATC I, 2F) 
CONTINUE 
PRINT  12  3,K 
FORMAT  C '0'  ,85X,  '  ITERATION  =  ',12) 

Qk  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

C  CALCULATES  REAL   REACTIVE  BUSS  POWERS 

Qkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

DO  150  P=2,IB 

PSUM=0 

QSUM=0 

DO  140  Q=l , IB  '       . 


59 

6  0 


70 


75 
00 
122 
123 


P1  =  REALCLCP)  )^-CREALCHCQ)  )"^REALCYCP,  Q)  )-AIMAGCECQ)  )*AI 
MAGCY 


1 CP, Q) ) )  +  AIMAGCEC  P) )'U  AIMAGCECQ) )*REALC YC  P, Q)  )  +  REALCE 


99 


(Q)) 


2AIMAG(Y(P,Q))) 


Q1=AIMAG(E(P) )*(REAL(E(Q) ) *RE AL ( Y ( P , Q ) )- AIM AG( E ( Q ) ) * A 
IMAG( Y 

1(P,Q)))-REAL(E(P))=UAIMAG(E(Q)  )  *  RE  AL(  Y  (  P  ,  Q  )  )  +  RE  AL(  E  ( 
Q))*AI 


140 


145 
150 

C  PRINT 
155 


2MAG(Y 

PSUM=P 

QSUM=Q 

CO NT IN 

PCAL(P 

QCAL(P 

PRINT 

FORMAT 

CONTIN 

A  A  *  -/c  *  A 

OUT  VA 

A  A  >V  >V  A  A 

PRINT 
FORMAT 
DO  158 
PRINT 
FORMAT 
CONTIN 
M=0 

DO  160 
EP=0.0 
PDIFF( 
QDIFF( 
TDIFF= 
IF(TDI 
M=l 
160       CONTIN 

C***AA*AAA**A* 

C  PRINT  OUT  VA 

PRINT 
165       FORMAT 


157 
158 


(P,Q))) 

SUM+Pl 

SUM+Ql 

UE 

)=PSUM 

)=QSUM 

145, PCAL(P) ,OCAL(P) 

(2F) 

UE 

LUES  OF  CALCULATED  BUSS  POWERS 

155, K 

('0','REAL    AND    REACTIVE     BUSS     POV/ERS    K=  '  ,  1 2  ) 

1=2  ,  IB 
157 ,PCAL(I) ,nCAL(I) 

( '0' ,2F) 
UE 

P  =  2  ,  IB 
1 

P)=PSCH(P)-PCAL(P) 
P)=QSCH(P)-QCAL(P) 
ABS(PDIFF(P))+ABS(QDIFF(P)) 
FF.LT.EP)  GO  TO  160 

UE 

LUES  OF  DIFFERENTIAL  BUSS  POWERS 

165, K 

( '0' ,' DIFFERENTIAL  PQ    K=',3X,I2) 


DO  158  1  =  2, IB 

PRINT  167,PDIFF(I),QDIFF(I) 

FORMAT  ( '0' , 2F) 

CONTINUE 

IF(M.EQ.O)  GO  TO  400 

DO  170  1  =  2, IB 

S( I )  =  CMPLX(PCAL( I )  ,QCAL( I )  ) 

CONTINUE 

C   CALCULATE  BUSS  CURRENTS  I(P)  EXCEPT  AT  SLACK  BUS 


167 
168 


170 


100 


Qi<  -M  i<  -Jc  -k  -k  -k  -k  -k  -k  -k  -k  :k  :k  •>'<  -k  -k  -k  -k  :k  -k  -f^  -k  it  i<  ii  :l!  ■J^  -k  ^  -k  -k  -k  i<  -k  -k  if  -k  -k  :k  -k  it  -k  -k 

DO  175  1=2, IB 

IC(I)=CONJG(S( 1) )/COHJG(E( I) ) 

PRINT  17  3,IC(I) 
173       FORMAT(F) 
17  5       CONTINUE 

Q-ki(kk-kk-kii:kki(-kkkk-k-k-kk-kkkk-k-kk-k-kk-ky';k-kkk-kkk-kk-k'k 

C   CALCULATE  JACOBIAN  USING  PARTITIOr;S  '  J  1  ,  J  2  ,  J  3  ,  J4  ' 

Qk  -k  -k  it  it  k  -k  i<  k  -k  k  k  k  k  -k  -k  k  -k  k  k  -k  k  k  k  k  k  ic  k  k  i<  k  k  k  -k  k  k  k  i<  it  i<  k  k 

C       JACOBIAN     'Jl' 

Qkkkitititisiti<i<kkiikk 

IBDJ=2'^IB-2 

IBJ=IB-1 

DO  190  M=l , IBJ 

DO  180  N=1,IBJ 

P  =  M+1 

Q  =  N+1 

JAC(M,N)  =  (REAL(E(P))'>REAL(  Y(P,  Q  )  )  )  + (  AIM  AG  (  E  (  P  )  )-^AIMAG 
(Y(P,Q))) 

IF(P.NE.Q)  GO  TO  177 

JAC(M,N)=JAC(M,N)+REAL(IC(P) ) 

177  CONTINUE 

C  PRINT  178 , JAC(M,N)  ,E(P)  , Y(P, Q) 

178  F0RMAT(5F) 
180  CONTINUE 
190  CONTINUE 

ititickkkkkkkkkkkkkkkkkkkkkkkkiskkkkkitkkkkkkkitiikkkitkkiikkkkkkkkii 

C       JACOBIAN     'J2' 

kkkkitkkkkkkkitkkkkkkki<itkkkkki<kkki<kkkki<kki<i!kkkkkkkkkkkkkkkki< 

DO  210  M=l , IBJ 
DO  200  N=IB,IBDJ 
P  =  M+1 
Q=H-IB+2 


JAC(M,N)=-REAL(E(P))'^AIMAG(Y(P,Q)  )  +  AIMAG(E(P)  )*REAL(Y 

(P,Q)) 

IF(P.NE.Q)  GO  TO  19  5 

J  A  C  (  M  ,  N  )  =  J  A  C  (  M  ,  10  +  A I M  A  G  (  I  C  (  P  )  ) 
195       CONTINUE 
200       CONTINUE 
210       CONTINUE 

kkkkkkkkkkkitkkkkkkkkkkkkkkkickkkkkkkkkitkitkkitkkkkkkkkkkkkkk 

C       JACOBIAN  'J3' 

kkkkkkkkkkkkkkkititkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkitkkkkkkk 

DO  2  30  M=IB,IBnj 
DO  220  N=l ,IBJ 
P=M-IB+2 
Q  =  N+1 

JAC(M,N)=-REAL(E(  P)  )*AIMAG(Y(P,Q)  )  +  AIMAG(E(P)  )--M<EAL(Y   . 


101 


(P,Q)) 

IF(P.I,'E.Q)  GO  TO  215 

JAC(M, N)=JAC(M, N)-AIMAG( IC( P)  ) 
215       CONTINUE 
220       CONTINUE 
2  30       CONTINUE 

C   JACOBIAN  'J4' 

*AAAAAA*AAAAA***A**A**A**A*AAAA**AAA*A*A**A***AA**AA*A 


DO  250  M=IB,IBDJ 
DO  240  N=IB,IBDJ 
P=M-IB+2 
Q=N-IB+2 


JAC(M,N)=-AIMAG(E(P))*AI]1AG(Y(P,Q))-REAL(E(P))'''REAL(Y 
(P.Q)) 


2  35 
240 
250 

C 

3  33 
334 

*  *  A  A  A 

*  A  A  A  A 

C   CO 


AAA 
AAA 

MPU 


A  A  A  A  A  A  A  A 

c 

344 
345 

AAAAAAAA 

C   COMPU 

AAAAAAAA 


IF(  P 
JAC( 

CONT 
CONT 
CONT 
DO  3 
PR  IN 
FORM 
CONT 

A  A  A  A 
A  A  A  A 

TE  I 
CALL 

A  A  A  A 

P  R  I  i^J 
FORM 
CONT 

A  A  A  A 


.  Tn'  E  . 
M  ,  N  ) 
INUE 
INUE 
INUE 
3  4  M 
T  33 
AT(  1 
INUE 

A  A  A  A 
A  A  A  A 

NVER 
MIN 

A  A  A  A 


Q)  GO  TO  235 

=  JAC(M, N)  +  REAL( IC(P)  ) 


=1 , IBDJ 

3, ( JAC(M,N) ,N=1 , IBDJ) 

8F) 


AAAAAAAAAAAAAAAAAAAAAAA 
AAAAAAAAAAAAAAAAAAAAAAA 


A  A  A  A  A 
A  A  A  A  A 

SE   OF 

V(  JAC 

A  A  A  A  A 

4  ,  ( JAC(M,N)  ,N=1  , IBDJ) 
3F) 


A  A  A  A  A 
A  A  A  A  A 


A  A  A  A  A  A 
A  A  A  A  A  A 


A  A  A  A  A 
A  A  A  A  A 


JACOBIAN 
,  6  8, IBDJ) 

A AAAAAAAA A AAAAAAAAAAAAAAAAAAAAAAAAAA 


TE  D 

AAAAAAAA 


T  34 
AT(  1 

INUE 

AAAAAAAAAAAAAAA  A  >':  AAAAAAAAAAAAAAAAAAAAAAAAAA 

IFER 


IN  VOLTAGE  USING  JACOBIAN  IN 

A  A  A  A  A 


ENCE 

AAAAAAAAAAAAAAAAAAAAAAAAAAAA 


VERSE 

A  A  A  A  A  A 


AAA 

'JAC 
A  A 


346 


349 


DO  3 
M=I- 
PQVE 
QQVE 
CONT 
J  =  IB 
DO  3 
E  F  D  V 
EFDQ 
EIDD 
EIQV 
DO  3 
E  F  D  V 
CONT 
DO  3 


=  P  D I  F 
=  O  D I  F 


1=1  ,  J 


46  1=2 , IB 

1 

C(M) 

C(M) 

INUE 

-1 

5  1  M 

D  =  0 

D  =  0 

=  0 

=  0 

49  N 

D  =  EF 

IN'UE 

50 


F(I) 
F(I) 


1=1  ,  J 
DVD+J 


IB, 


AC(M,N)'^PQVEC(N) 
IBDJ 
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NQ=N-J 
EFDQD=EFDQD+JAC(M,N)*OQVEC(NQ) 

350  CONTINUE 
EFDV(M)=EFDQD+EFDVD 

351  CONTINUE 

DO  360  M=IB,IBDJ 

EIDD=0 

EIQV=0 

DO  3  58  N=1,J 

EIDD=EIDD+JAC(M,N)*PQVEC(N) 

358  CONTINUE 

DO  359  N=IB,IBDJ 

NQ=N-J 

EIQV=EIQV+JAC(M,K)*QQVEC(NQ) 

359  CONTINUE 
EIDV(M)=EIQV+EIDD 

360  CONTINUE 

C   COMPUTE  NEW  BUS  VOLTAGES  USING  'EFDV'  EFDIF FERENTI AL 
VECTOR 


368 
370 


390 
395 

CONVE 
A  00 
500 

A  *  A  *  * 

c 
c 
c 

•k  :k  ■k  -k  •/< 


503 


DO 
M  = 
N  = 
IB 
IF 
ED 
PR 
FO 
E( 
CO 
IF 
K  = 
GO 
TY 
FO 
RGE'  ) 
CO 
CL 

■k  k  k  k  i< 

CO 
CO 
CO 

*  A  *  *  A 

OP 

s( 

FO 
DO 
DO 
IF 
IF 


*AAAAAAA*>V*AAAA*>'cAAAAAycAAA**A*AAAAAA*AAAAA 

P=2 , IB 


A  A  A  A 

370 
P-1 
11+ J 
Z  =  2* 
(N.G 
IFF( 
INT 
RMAT 
P)  =  E 
NTIN 
(K.E 
K+1 

TO 
PE  3' 
RMAT( '  AN  ERROR  IS  MADE  IN  DATA,  LOAD  FLOW  DID  NOT 


(IB-1) 

T.IBZ)  GO  TO  370 

P)=CMPLX(EFDV(M) ,EIDV(N) ) 

368,P,EDIFF(P) 

(I,2F) 

(P)+EDIFF(P) 

UE 

Q. 10)G0  TO  390 

122 
95 


NTINUE 
OSE(UNIT=3  0) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

MPUTE  POWER  FLOWS  BUS  P  TO  Q 

MPUTE  CURRENT  FLOWS  AND  POWER  LOSSES 

MPUTE  SLACK  BUS  POWER 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

EN(UNIT=3  1,FILE='DSK') 
1)  =  0 
RMAT( I) 

520  P=1,IB 

515  Q=l ,IB 
(P.EQ.Q)  GO  TO  515 
( Y(P,Q) .EQ. (0.0,0.0) )  GO  TO  515 
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PP  =  P 

C  IF  YOU  ARE  NOT  NECLF.CT  SHUNT  ADMITT  REMOVE  THE  NEXT 

C 

*  A  A  *  A  A  *  Ss  A  Vc  *  V:  ;k  A  A  *  A  *  :fe  *  :fc  Vc  *  *  A  *  A  A  A  A  A  :fc  A  A-  *  :St  -^  :(.'  ^c  -.fc  A  X  ;Ic  A  :1c  ik  >!:  *  Vt  *  A  7k  *  iic  A  -A  * 

PFLO  =  CONJG(E(P)  )'>(  (E(P)-E(Q)  )-^Yl  (  P,  Q)  ) 
C  1  +  (C0NJG(E(P))AE(P)''=YS(P,Q)) 

PFL01  =  C0NJG(E(Q))*(  (E(Q)-E(P)  )-''Yl  (P,Q)  ) 
C  1  +  (C0NJG(E(Q))*E(Q)='=YS(P,Q)  ) 

CRFLO(P,Q)=PFLO/CONJG(E(P)) 

PL0S=PFL0+PFL01 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

TAKE  THE  CONJUGATE  OF  THE  POWER  FLOVJS  FOR  FORMAT 


A  A  A  A  A 

c 

P-JQ 

A  A  A  A  A 


50  5 

509 
13  U  S  '  , 

510 
515 
52  0 

540 


A  A  A  A  A 

PL 
PF 
IF 
S( 
V^R 
WR 
FO 
'T 
IF 
FO 
CO 
CO 
WR 
FO 
DO 


AAAAAAAAAAAAAvVAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


12 


WRITE( 33, 5 
IC 
FO 


543 

545 
5  50 

555 


570 
5  80 
f>00 


CO 
WR 
FO 
CL 
DO 
DO 
S( 
PF 
IF 
PF 
WR 
FO 
CO 
CO 
CL 


OS  =  C 
LO  =  C 
(P.N 
1)  =  P 
ITE( 
ITE( 
RHAT 

oM 

LOW 
RMAT 
NT  IN 
NT  IN 
ITE( 
RMAT 
550 


43)M 
RFLO 
RMAT 

NTIN 
ITE( 
RMAT 
OSE( 
600 
580 
2)  =  E 
AD  =  C 
(PFA 
AC  =  - 
ITE( 
RMAT 
NTIN 
NTIN 

0  s  !•:  ( 


ONJG 
ONJG 
E.  1) 
FLO+ 
31,5 
31,5 
(6X, 
2, 'I 
IS', 
(5X. 
UE 
UE 

31,5 

(6X, 

M=l 

DO 

IF 

,N, 

(  M  ,  N 

(21, 

CO 

UE 

33,  5 
('  0 
UNIT 
M=l 
N=l 
(M)--v 
ARS( 
D.  EQ 
(  REA 
3  1,5 
(4X, 
UE 
UE 

u  r;  I T 


( PLOS) 
( PFLO) 

GO  TO  505 
5(1) 

09  )P, Q, PFLO, CRFLO( P, Q)  , PLOS 
10)PP,E(PP) 
'POWER  FLOW 
S'  ,  2F, 2X, ' CURRENT 
2F,  ' POWER  LOSS  IS'  , 2F) 
' BUS' , 12 , 'VOLTAGE  IS',2F) 


4  0)  S ( 1 ) 

'POWER  AT  SLACK  BUS  1  IS',2F) 

,IB 

5  4  5     N  = 1  , I  3 
(CABS(CRFLO(M,N)).GT. 0.0 00001) 


) 
2F) 

NTINUE 


55) 

=  33 
,IB 
,IB 
CON 
S(2 
.0) 
L(S 
70) 
'  PO 


0  0.0    0.0') 

) 


JG(CRFL0(M, N) ) 

)) 

GO    TO    580 
(2)))/(PFAD) 

M  ,  N  ,  P  F  A  C 

WER  FACTOR  BUS' , 13 , ' TO' , 13 , ' IS' , F) 


-31) 


104 


10 
15 


20 


25 


3  0 
35 

38 


40 
45 

50 
55 


6  0 
62 
65 


70 
75 

8  0 

100 


STOP 

END 

TITLE  :MIKV  ;  SUBROUTINE  TO  I?^'VERT 

SUBROUTINE  MINV( A,MDIM , N) 

DIMENSION  A(MDIM,N) 

INTEGER  L(  100)  ,?-l(100) 

CONVERTED  FROM  SSP  ROUTNE  MINV 


DO  8  0  K=l ,N 

L(K)=K 

M(K)=K 

BIGA=A(K,K) 

DO  20  J=K,N 

DO  2  0  I=K,N 

IF(ABS(BIGA)-ABS(A(I,J)))  15,2  0,20 

BIGA=A( I, J) 

L(K)=I 

M(K)=J 

CONTINUE 

J=L(K) 

IF (J-K)  3  5,35,25 

DO  3  0  I=1,N 

HOLD=-A( K, I) 

A(K,I)=A( J,I) 

A(  J , I)  =  HOLD 

I=M(K) 

I F ( I - K )  4  5,45,38 

DO  4  0  J=1,N 

H0LD=-A( J, K) 

A( J,K)=A( J,I) 

A( J, I)=H0LD 

DO  55  1  =  1  ,N 

IF(I-K)  50,55,50 

A(I,K)=A(I,K)/(-BIGA) 

CONTINUE 

DO  6  5  1  =  1  ,N 

HOLD=A( I ,K) 

DO  65  J=l  ,  N 

I F ( I - K )  6  0,65,60 

IF(J-;0  6  2,65,62 

A( I, J)-HOLD*A( K,J)+A( I , J) 

CONTINUE 

DO  7  5  J=1,N 

I  F  (  J  -  IC  )  7  0,75,70 

A(K, J)=A(K, J)/BIGA 

CONTINUE 

A(K, K)  =  (  1  .0000000)/BIGA 

C0NTl;>JUE 

K=N 

K  =  K-1 


A  REAL  SQUARE  MATRIX 


105 


I F ( K )  15  0,150,105 
105       I=L(K) 

I  F ( I - K )  12  0,120,108 
108       DO  110  J=l ,N 

1I0LD  =  A(  J,K) 

A( J,K)=-A( J,I) 
110       A(J,I)=H0LD 
120       J=M(K) 

IF( J-IO  100, 100, 125 
125      DO  130  1=1, N 

H0LD=A(K, I) 

A(K,I)=-A( J,I) 
130       A(J,I)=UOLD 

GO  TO  100 
150       RETURN 

END 

C         TITLE:ZBUS. FA ; ACCEPTS  DATA  FROM  TEST. DAT  AND  DSK.XFR 

FROM 

C  LOAD  FLOW  PROGRAM  LDFL0.F4  AND  FORMS  ZBUS  BY 

C  ALGORITHM  THIS  WILL  BE  USED  TO  COMPUTE  ALL 

FAULT 

C  CURRENTS. LOAD  FLOW  COMPUTED  MAXIMUM  LOAD 

CURRENTS. 

C         READ  DATA  FROM  TEST. DAT 

DIMENSION 
ZD(30O,3),ZM(30O,2),RID(3  0,3  0,3),ZG(3  0),ZN(30), 

1ID(200),RI(30),ZB(3  0,30),ZB0(30,30),ZMU(5,5),IDD(5), 
ZL(30) , 

2ZBR(30),ZG0(30),ZL0(30),EB(30) 

COMPLEX  ZD,ZM,ZG,T1,ZN,ZN1,ZNI,ZNJ,ZG1,Z0, 

1ZB,ZB0,ZMU,ZBB,ZL,ZLL,ZBR,ZG0,ZL0,ZLL0,EB 

NNB  =  1 

EQUIVALENCE(ZBR, ZLO,EB) 

0PEN(UNIT=3  1,FILE='TEST.DAT') 
C         OPEN(UNIT=32,FILE=' DSK.XFR') 

READ(31 , 10)IB, EBl 
10        FORMAT(I,F) 

20        READ(31,30)NB,I,J,JJ,T1,ZG1,ZNI,ZNJ 
3  0        F0RMAT(4I , 8F) 

I  F  (  N  B  .  E  f ) .  0  )  G  0  T  0  A  0 

DN  =  NB 

NN  =  NB 

ZD(NN, 1 )=CMPLX(DN, 0.0) 

Z  D  (  N  N  ,  2  )  =  T  1 

Z  D  (  N  N  ,  3  )  =  Z  G  1  +  3  *  Z  N  (  I  ) 

Z  !;  (I  )  =  Z  N  I 
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40 
C 

43 
45 

C 

46 

48 
49 
5  0 


ZN( J)=ZNJ 

RID( I, J, JJ)=NN 

IF(NB.NE.0)GO  TO  20 

DO  4  5  1=1 , IB 

READ(32, 43)T1 

F0RMAT( 2F) 

CONTINUE 

DO  48  1=1 ,IB 

READ(32,46)J,T1 

F0RHAT(I,2F) 

EB( J)=T1 

CONTINUE 

READ(3  1 , 50)I,T1 ,XD1 

FORMATC I, 3F) 

IF(I.EQ.O)  GO  TO  60 

IF( CABS(T1 ) .EQ.O)  GO  TO  49 
C         T1=C0NJG(T1) 
C         ZG(I)=3*(CABS(EB( I) )A*2)/T1 

IF(XD1.NE.0)ZG(I )=CMPLX(0.0, XDl ) 

ZGO( I)=ZG(I) 

IF( I . NE.O)  GO  TO  4  9 

A*A*****^***5!cyc*>S:*****A5k***i!c**:S:*AA*A*;S:*A;!cAA5V  *********  A 

C         ENTER  THE  ZERO  SEQ  IMP  BUS    TO  REF( GROUND) 

**********  vV  *  *  *  **********************************>'<******* 

6  0        TYPE  6  2 

62        FORMATC'  ENTER  THE  SYNC   ZERO  GEQUEI.'CE  IMPEDENCE  FOR 

THE  LOADS 

lAT  EACH'/'  BUS, 3*ZN, MUST  BE  ADDED, IMPED  IN  PU..,') 

DO  65  1=1 , IB 

ACCEPT  64,1I,T1 

IF(II.EQ.0)GO  TO  66 

ZGO( II)=T1 

format( i, 2f) 
coi:tinue 

DO  68  1=1 , IB 

IF(CABS(ZGO( I)).En.O)ZGO(I)=ZG(I) 

CONTINUE 
******************************************************** 

C         ENTER  THE  MUTUALLY  COUPLED  ELEMENTS 
******************************************************** 

CL0SE(Ur;iT  =  3  1  ) 

TYPE  70 
70        FORMATC'  ENTER  THE  MUTUALLY  COUPLED  ELEMENTS  IN  PER 
UNIT'/' 

1  IN  THE  FOLLOWING  FORM AT , TE S T . D AT  HAS  THE  ELEMENT 
S..,'/' 

2  CELEMENT  )SP  TOCELEMENT  )SPC0  SEQ  VALUE)  ') 
75        ACCEPT  80, EL, ELL, ZO 

80        FORMATC 4 F) 

IFCEL. EO.O.O)  GO  TO  100 
11=11+1 


6  4 
65 
6  6 

6  8 
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ZM(1I,  1)  =  CMPLX(EL,ELL) 
ZH(II, 2)=Z0 
GO  TO  7  5 


c 

NEGATIVt: 

c 

******** 
******** 

100 


ALL  ELEMENT  DATA  PREPARED  PERFORM  POSITIVE  AND 


SEQUENCE  COMPUTATIONS  TO  FORM  ZBUS 


************************************************ 
************************************************ 

DO  500  IEL=1,NN 
BRA  =  0 
ID( IEL)=1 
DO  120  1=1 , 30 

DO  120  J=l , 30 

DO  120  JJ  =  1  ,  3 

IF(RID(I, J, JJ) .NE. lEL)  GO 


TO 


120 
110 
115 


110 
1  15 


115 


IF(RI(I) .EQ.0)GO  TO 

IF( RI( J)  .En.0)GO  TO 

BRA=1.0 

IPP=MINO( I , J) 

IQQ=MAXO(I , J) 

GO  TO  130 

IF(RI( J ) .EQ.0)GO  TO 

RI(  I)  =  l  .0 
IPP  =  J 
IQQ  =  I 
GO  TO  130 
RI( J)=l  .0 
R1(I)  =  1  .0 
IPP=MIN0( I, J) 
IQQ=MAX0(I , J) 
GO  TO  130 


CONTINUE 


120 
***************************************************** 

C         BRANCH  OR  LINK  HAS  BEEN  DETERMINED  BRA=0: BRANCH 

*************************************•,■:******************* 


130 

lyi 


186 
190 
195 


DO  195  1L=1,IB 

ZBB=CMPLX(0 .0,0.0) 

IF(BRA.NE.O)  GO  TO  136 

IF( IQQ.EQ. IL)  GO  TO  190 

ZBR(IL)=ZB(IPP,IL) 

GO  TO  190 

ZL(IL)  =  ZB(IPP,IL)-ZB(ir)(^,IL) 

CONTINUE 
CONTINUE 

ZBB=CMPLX(0.0, 0.0) 
IF(BRA.NE.O)  GO  TO  2  50 
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2  00 
240 


250 


DO  2  00  1L=1,IB 

ZB( IQQ, IL)=ZBK( IL) 

ZB(IL,IQQ)=ZBR(IL) 
CONTINUE 

ZB(I0Q,IQQ)=ZB(IPP,IQQ)+ZD(IEL,2) 
ZBB=CMPLX(0 .0,0.0) 
GO  TO  300 

ZLL=ZL(IPP)-ZL(IQQ)+ZD(IEL,2) 
ZBB=CMPLX(0. 0,0.0) 
DO  290  IP=1,IB 
DO  290  IQ=1,IB 

IF( IP.GT. IO)GO  TO  290 


ZB(IP,IQ)=ZB(IP,IQ)-(ZL(1P)«Z 
L(IQ) 

1/ZLL) 


IF(IP.EQ.IQ)  GO  TO  290 
ZB(IQ,IP)=ZB(IP,IQ) 


350 


380 
400 


290       CONTINUE 
300  DO  350  1=1,5 

DO  350  J  =  l  ,  5 

ZMU( I , J)=CMPLX(0.0,0.0) 
IDD( I)=0 
CONTINUE 
DO  400  1=1, IB 
DO  400  J=1,IB 

PR.INT  3<S0,I,J,ZB(I,J) 
FORMATC 21, 2F) 
CONTINUE 
NNB=1 
500       CONTINUE 

C  ZBUS  POS/NEG  IS  FORMED 

DO  50  5  1=1 , NN 

ID(I)=0 
505       CONTIN'UE 

DO  5  10  1=1 , IB 

RI( I)  =  0  .0 
5  10       CONTINUE 

DO  1000  IEL=1,NN 
BRA  =  0 
ID( IEL)=1 
NNB=1 
DO  620  1=1 , 30 

DO  620  J=l , 30 

DO  6  2  0  JJ=1 , 3 

IF(  RI1)(  1  ,  ,1,  JJ)  .NE.  lEL)  GO  TO 


620 
610 


IF(  RI(  I)  .  E').0)GO  TO 
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615 


610 
615 


615 


1F( Rl( J) .EQ.O)  GO  TO 

BHA=1 .0 

IPP=MIN0(I, J) 

IQQ=MAX0(1 , J) 

GO  TO  625 

IF( RI( J) .EQ.O)  GO  TO 

RI(1)  =  1  .0 
IPP  =  J 
IQQ=I 
GO  TO  625 
R1(I  )  =  1  .0 
RI(J)  =  1  .0 
IPP=I 
IQQ  =  J 
GO  TO  625 
CONTINUE 


620 
***** 

C 
********************************************************* 


SEARCH  FOR  MUTUAL  COUPLING  IN  ZERO  SEQUENCE 


625 


1DD( 1 )=IEL 

DO  650  1=1 , NN 


630 


DM=IEL 

IF(REAL(ZM(  I  ,  1 ))  .EQ.DM)MUT=1 

IF( A1MAG( ZM( 1 , 1) ) . EQ . DM)HUT=1 

IF(MUT.NE.1)G0  TO  650 

MUT  =  0 

II1=IF1X(REAL(ZM( 1 , 1 ) ) ) 

II2=IFIX( A1MAG(ZM( 1,1))) 

IF( ID(II1  )  .EQ.0)GO  TO  650 

IF(ID(II2).EQ.0)G0  TO  650 

IF( III  .EQ. IEL)GO  TO  630 

113=111 

II  1  =  112 

112=113 

I  n  D  (  1  )  =  I  1  1 

NNB=NNB+1 

1DD(NNB)=1I2 

DM=1I2 

DO  64  5  J1=1,NN 


1F( REAL(ZM( Jl , 1) ) .EQ.DM)HUT=1 

1F( A1HAG(ZM( J  1  ,  1 ) )  .EQ. DM)MUT  = 
1 


1F( MUT. NE. 1 )  GO  TO  645 
MUT  =  0 

I111=1F1X( KEAL(ZM( Jl , 1 ) ) ) 
1 1  I  2  =  1  F  1  X  (  A  I  M  AG  (  Z  M(  J  1  ,  1  )  )  ) 
1F(I1KIII1).EQ-^)G0  TO  6  45 
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643 


645 
650 


IF(ID(III2) .EQ.0)GO  TO  645 

IF(  nil  .EQ.DM)  GO  TO  643 

1113=1111 

1111=1112 

1112=1111 

1F( 1112 .EQ. Ill )  GO  TO  645 

NNB=NNB+1 

IDD(NNB)=I112 
CONTINUE 
CONTINUE 
DO  660  1=1 ,KNB 

ZMU(I,1)=ZU(1DD(1),3) 
DO  660   J=l , NNB 

DD1=1DD(1) 

DD2=IDD( J) 

DO  66  0   KK=1 ,NN 


IF(ZM(KK, 1) .EQ.CMPLX( 
DDl ,DD2) ) 

1MUT=1 


1F(ZM( KK, 1 ) . EQ.CMPLX(DD2,DD1 ) 

)  MUT=1 

IF(MUT.NE. 1 )G0  TO  66  0 
MUT  =  0 

ZMU( I , J)  =  ZM(KK,  2) 
ZMU( J, 1 )  =  ZM(  KK,  2) 

660  CONTINUE 

C         ZMU  IS  FORMED 

DO  6  6  5  1  =  1  ,  N !>! B 

DO  665  J=l , NNB 

PRINT  663, 1, J, ZMU( 1 , J) 
663  F0RMAT(21, 2F) 

665  CONTINUE 

AA:!cA>'cA*AAA5lcyr5lc*AAA*A*A****AA*:!c**V!cA-A*5^*5Sr*A>V****5'cA**A:!c55:)!cA*A* 

CALL  CMX1NV(ZMU, 5, NNB) 

DO  668  1=1, NNB 

PRINT  667,(ZMU(I,J),J=1,NNB) 

667  FORMAT(IOF) 

668  CONTINUE 

DO  69  5  IL=1 , IB 

IF(NNB.EQ.l)  GO  TO  681 

ZBB=CMPLX(0 .0,0.0) 

DO  6  80   J4=l , NNB 
J=J4+1 
I=1DD( J) 

IF(l.EQ.O)  GO  TO  680 
DO  67  0  I  1  =  1  ,  30 


Ill 


DO  670  12=1,30 
DO  r,  7  0 


13=1  ,  3 

I  F  (  R  I 
0(11,12,13 

1)  .IJE.  I)  GO  TO  670 


IP  =  I1 

IQ=I2 

6  70 
680 
681 


686 
6  90 
695 


700 


720 


ZBB=(ZMU(1,J)/ZMU(1,1))*(ZB0(IP,IL)-ZB0(IQ,IL))+ZBB 

CONTINUE 
CONTINUE 

IF(BRA.NE.O)  GO  TO  686 
IF(IQQ.EQ.IL)  GO  TO  690 
ZBR(IL)=ZBO(IPP,IL)+ZBB 
GO  TO  690 

ZL(IL)=ZBB+ZB0(IPP, IL)-ZBO(IQQ,IL) 
CONTINUE 
CONTINUE 

ZBB=CMPLX(0 .0, 0.0) 
IF(BRA.NE.O)  GO  TO  750 
DO  7  00  IL=1,IB 

ZBOdQQ,  IL)  =  ZBR(  IL) 
ZBO(IL, IQQ)  =  ZBR(  IL) 
CONTINUE 

IF(NNB.EQ.l)  GO  TO  7  40 
DO  730  J=2,NNB 

I=IDD( J) 

DO  720   11=1 , 30 

DO  72  0  12=1,30 

DO  720  13=1,3 
IF(RID(I 1  , 12,  13)  .NE. I)  GO  TO  720 

IP  =  I1 
IQ  =  I2 
CONTINUE 


ZBB=(ZMU(  1  ,  J)'==(ZBO(IP,  lQO)-ZB0(lg,  IQQ 
))/ 

1ZMU( 1  ,  1  ))  +  ZBB 
7  30  CONTINUE 

740  ZBO(IQO,IQQ)=ZBO(IPP,IQQ)+ZBB+1.0/ZMU(1,1) 

ZBB=CMPLX(0 .0,0.0) 
GO  TO  800 
7  50  IF(NNB.EQ.l)  GO  TO  78  1 

DO  7  80  J=2,NNB 

I=IDD( J) 

DO  770   11=1,30 

DO  770   12=1,3  0 

DO  7  70   13=1,3 


112 


770 


112  ,  13)  .NE. I)GO  TO  770 


C0NTi:ajE 


1 F  (  R  1  D  (  1  1  , 

1P=I  1 
1 Q  =  I  2 


Z  B  B= ( Z  MU (I,J)*(ZL(IP)-ZL(IQ))/ZMU(1,1 

)  )  +  ZBB 

780  CONTINUE 

78  1  ZLL  =  ZL(IPP)-ZL(IQQ)  +  ZBB+( 1  . 0 / ZMU ( 1  ,  1  )  ) 

ZBB=CMPLX(0. 0,0.0) 

DO  790  IP=1  ,  IB 

DO  790  IQ=1,IB 


IF( IP.GT. IQ)G0  TO  790 


ZB0(IP,1Q)=ZB0(1P,1Q)-(ZL(IP) 
*ZL( IQ) 

1 /ZLL) 


IF(IP.EQ.IQ)  GO  TO  790 
ZBOdQ,  IP)  =  ZBO(  IP,  IQ) 


7  90 
800 


CONTINUE 


50 


830 
900 


DO  850  1=1 , 5 
DO  850  J=l , 5 

ZMU( I, J)=CMPLX(0 .0,0.0) 

IDD(I)=0 
CONTINUE 
DO  900  1=1 , IB 
DO  900  J=l , IB 

PRINT  880, I, J,ZBO(I , J) 

F0RMAT(2I , 2F) 
CONTINUE 
NNB  =  1 
1000      CONTINUE 

A  A  A  A  A  A 

C  ZBUS  ZERO  IS  FORMED 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

OPE N(UNIT=33, FILE ='ZBUS. FOR') 

DO  1010  1=1 , IB 

DO  10  10  J=l , IB 

URITE(3  3 ,1005)I,J,ZB0(I,J),ZB(I,J),ZB(I,J) 

FORMATC  21 , 6F) 

PRINT  1007 

FORMATC 3X, ' BUS   TO  BUS  ',3X,'ZER0  SEQ  ',5X,'P0S 


100  5 

1007 
SEQ'  ) 

1008 
1010 

A  A  A  A  A  A  * 

c 


PRINT  1008, I , J, ZBO(I , J) ,ZB( I , J) 

FORMATC  216 , 4F) 

CONTINUE 

CLOSECUNIT=32) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

NOW  TO  ADD  ALL  LOADS.  ALL  LOADS  WILL  RE  LINKS  WITH 
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NO 

C         MUTUAL  COUPLING. 

DO  1200  1  =  1, IB 

IF(  CABP.(ZG(I)  )  .EQ.0)GO  TO  1200 
DO  1 100  IL=1 , IB 

ZL( IL)=-ZB(I, IL) 
ZL0(IL)=-ZB0(I, IL) 
1100  CONTINUE 

ZLL=-ZL(I)+ZG(I ) 
ZLL0=-ZL0(I)+ZG0(1) 
DO  1150  J=l , IB 

DO  1150  K=l ,IB 

IF( J .GT.K)G0  TO  1150 


ZB(J,K)  =  ZB(J,K)-ZL(J)'«ZL(K)/Z 
LL 


ZB(K, J)  =  ZB(  J,K) 


ZBO(  J,K)  =  ZBO(  J,K)-ZLO(  J)'^ZLO( 
K)/ZLLO 


ZBO(K, J)=ZBO( J, K) 


1  150 
1  200 


1250 
1300 


1310 

1315 
1  320 


CONTI 
OPEN( 
DO  13 

WRITE 
FORMA 
CONTI 
CLOSE 
CLOSE 
OPEN( 
DO  14 


C 

NUE 

UNIT 

00  I 

D 

(34, 

T(  21 

NUE 

(UNI 

(UNI 

UNIT 

00  I 

E 

I 

E 

I 

I 

D 

D 

D 

I 

I 

G 

I 

I 

G 

I 

I 

C 


ONTINUE 


=34, FILE='ZBUS.DAT' ) 

=  1  ,  IB 

0  1300  J=l , IB 

1250)1, J, ZB(I, J) ,ZBO(I , J) 

.4F) 


T  =  34 
T  =  33 
=  36, 
=  1,N 
L1=R 
F(EL 
L2  =  A 
EL1  = 
EL2  = 
0  13 
0  13 
0  13 
F(RI 
F(RI 
0  TO 
BII  = 
BIJ  = 
0  TO 
BJI  = 
BJJ  = 
ONTI 


) 
) 

FI 

N 

EA 

1  . 

IM 

IF 

IF 

20 

20 

20 

D( 

D( 

1 
II 
Jl 

1 
II 
J  1 
NU 


LE='HUT. CO' ) 


L(ZM 

EQ.O 

AG(Z 

IX(E 

IX(E 

11  = 

Jl  = 

Kl  = 

II,. J 

I  1  ,  J 

320 


320 


(1,1 
.0)G 
M(I, 
LI) 
L2) 
1  ,  IB 
1  ,  IB 
1  ,  IB 
1,K1 
1  ,  Kl 


)) 

0  TO 

D) 


1322 


)  .EQ.ELl )G0  TO 
).EQ.EL2)G0  TO 


13  10 
1315 


J 
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1322 


IF(EL1 . EQ.0)GO  TO  1400 


\;RITP:(36,1350)EL1,EL2,1BII,IBIJ,IBJ1,IBJJ,ZM( 

1,2) 

1350 

1400 


1450 


10 
15 


20 


25 


30 
35 

38 


40 
45 

50 
55 


CO  NT  I 

WRITE 

FORMA 

CLOSE 

STOP 

END 

TITLE 

SUBRO 

COM  PL 

INTEG 


FORMAT( 2F,4I3 , 2F) 
NUE 
(36,1450) 

T( '   0.0  0.0  0  0  0  0  0.0  0.0') 
(Ut'IT  =  36) 


:CMXINV; SUBROUTINE  FOR  INVERT  COMPLEX  MATRIX 
UTINE  CMXINVC A, MDIM,N) 
EX  A( MDIM,N) , BIGA, HOLD 
ER  L( 100) ,M( 100) 


CONVERTED  FROM  SSP  ROUTNE  MINV 


K=l  ,N 


25 


6  0 
62 


DO  00 

L(K)=K 

M(K)=K 

BIGA=A(K,K) 

DO  2  0  J=K,N 

DO  2  0  I=K,N 

IF( CABS( BIGA) 

BIGA=A( I, J) 

L(K)=I 

M(K)=J 

CONTINUE 

J=L(K) 

IF(J-K)  3  5,35 

DO  3  0  1=1,  N 

HOLD=-A(K, I) 

A(K,I)=A(J,I) 

A( J, I)=HOLD 

I=M(K) 

IF(I-K)  45,45 

DO  4  0  J=l  ,  N 

HOLD=-A( J,K) 

A( J,K)=A( J,I) 

A( J, I)=HOLD 

DO  5  5  1  =  1  ,N 

IF(I-IC)  50,55. 

A(  I ,  K)  =  A(  I , K)/(-BIGA) 

CONTINUE 

DO  6  5  1  =  1  ,N 

HOLD=A( I,K) 

DO  6  5  J=1,N 

I F ( I -  lO  6  0,65,60 

IF(J-K)  6  2,65,62 

A(I,  J)  =  H0LD''^A(  K,  J)  +  A(  I  ,  .1) 


CABS( A( I , J) ) )  15,20,20 


38 


50 


115 


65 

70 
75 

80 

100 
105 
108 


110 
120 

125 


CONTINUE 

^0  75  J=1,N 

IFCJ-K)  70,75,70 

A(K,J)=A(K,J)/BIGA 
CONTINUE 


CONTIN^E'-'"''°°"°'''-0'"'0000)/BIr;, 


1  30 
150 


K=N 
K  =  K-1 

IF(K)  150,150,105 
I=L(K) 

IF(I-K)  120,120,108 

DO  110  J=1,N 

HOLD=A( J,K) 

A(J,K)=-A(J,I) 

A( J, I )=HOLD 

J=M(K) 

I^^J-K)  100,100,125 
DO  130  1  =  1, N 
HOLD=A(K,I) 
A(K,I)=-A(j,i) 

A(J,I)=HOLD 
GO  TO  100 
RETURN 


TlTLE:FALT.F4;TlKEr"""""**'"""'""  — 


**  *i<^ 


END 

c 

DIMENSION 

REAnr-^/   ^N  .    ^^^  ^BUS.DAT') 
KLAD(34,5)  IB, El 

^         F0RMAT(I,F) 

IBB=IB*IB 

DO  7  1  =  1,  IB 

1^B(I)  =  CMPLX(E1,0.0) 
CONTINUE  ^ 

READ(34,10),NI,I  J  K 
FORMATCAI)      '-^'^.^^ 

IF(NI.E().0)GO  TO  15 
IF(K.EQ.2)l=-i 

^^'<^<-EQ.3)j  =  -j 
R1  =  I 

HJ  =  J 
lDELE(NI)=CMPLx(Ki,Rj) 


7 
9 
10 
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IF(NI . NE.0)CO  TO  9 
15        DO'30  11  =  1, IliB 

READ(3  5, 25) 1, J, ZBl ,ZB01 
25        F0RMAT(2I, 4F) 

ZB( I, J)=ZB1 

ZBO(I, J)=ZB01 
30        CONTINUE 

CLOSE(UNIT=35) 

OPEN(UNIT=3  7 , FILE='MUT. CO' ) 

OPEN(UNIT=38,FILE='FALT.DAT') 

OPEN(UNIT=39,FILE='EFLT.DAT') 

0PEN(UHIT=4  0, FILE='RELA.DAT' ) 

OPEN( UNIT=A 1 , FILE='FCUR.DAT' ) 

OPEN(UNIT  =  4  2  ,  FILE=' DSK.XFK'  ) 

OPEN(UNIT=4  3, FI LE= ' LCUR . DAT ' ) 

READ(40,32)  RMVAB,RKVBA, IRR 

32  F0RMAT(2F,I) 
RKVBA=RKVBA*1E3 

FIPU  =  RMVAB*  lE6/(  SQRT(  3  .0) --RKVBA) 

33  READ(40, 34) IRR 

34  FORMAT(I) 

IF( IRR.NE.O)GO  TO  33 

36  READ(40 , 37) J,  E 

37  FORMAT(I,F) 
IF(J.EQ.0)GO  TO  38 
E I  B  (  J  )  =  E 

GO  TO  3  6 

38  TYPE  35 

35  FORMAT('  ENTER  THE  FAULT 
IMPEDENC£-REAL(SP)IMAGINARY. . , ' ) 

ACCEPT  40, ZF 
40        F0RMAT(2F) 

DO  500  1FLT=1 , IB 

IF(  CABS(ZB(  IFLT,  IFLT))  .EQ.OGO  TO  500 
DO  50  1  =  1, IB 

IF( I.EQ. IFLT)GO  TO  45 

EFALT(  I,  2  ,  1  )  =  S0RT(3  .0)''=(EB(  I)-(  (ZB(  I  , 
IFLT)*EB 

I(IFLT) )/(ZF  +  ZB(IFLT,  IFLT)  ))  ) 

EFALT(  I,  1  ,  2)  =  -(SQRT(  3  .0)'-EB(  IFLT)*ZBO 
(I, IFLT)/ 

1(ZB0(1FLT,IFLT)+3*ZF+(2*ZB(IFLT,IFLT)))) 

EFALT(I,2,2)=SQRT(3.0)*EB(I)-((SQRT(3 
.0)*EB 

l(IFLT)>^ZB(I,IFLT))/(ZB0(IFLT,IFLT)  +  (2^>-(ZB(IFLT,IFLT) 
)+3*ZF))) 
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EFALTd,  3,  2)  =  EFALT(  1,  2,  2)-(SQRT(3.0)* 
EB(I)) 

GO  TO  50 
45 

EFALT(IFLT,2,1)  =  SQRT(3.0)''^(ZF*EB(1FLT 
)/(ZF+ 

1ZB(  IFLT, IFLT)  )  ) 

EFALTCIFLT,  1,2)=-SQRT(3.0)--^EB(IFLT)'^-Z 
BO (IFLT, 

1IFLT)/(ZB0(IFLT,IFLT)  +  (2^'ZB(IFLT,IFLT))  +  3*ZF) 

EFALTCIFLT, 2,2)=SQRT(3.0)*EB(IFLT)*(Z 
BO(IFLT, 

1IFLT)+ZB(IFLT,IFLT)+3*ZF)/(ZB0(IFLT,IFLT)+(2*ZB(IFLT 
,IFLT))+3* 

2ZF) 

EFALTCIFLT, 3,2)=-SQRTC3.0)*EBCIFLT)*Z 
BCIFLT.IFL 

1T)/CZB0CIFLT,IFLT)+C2*ZBCIFLT,IFLT))+3"ZF) 
50  CONTINUE 

REVJIND  3  4 
REWIND  37 

WRITEC39, 53) C IFLT, I, EFALTCI,  2,  1 ) ,EFALT(I,  1  ,  2) 

lEFALTC I, 2, 2), EFALTCI, 3, 2), 1=1, IB) 
53  FORMATC2I3,8F) 

READC34 , 55)11, E 
55  FORMATCl,F) 

DO  90  K=l , IBB 

READC3  4, 60)NN, I, J, KK, ZRl , ZNl , ZNI , ZNJ 
60  F0RMATC4I, 8F) 

IFCNN.EQ.0)GO  TO  91 

FCURCNN,2,1)=(EFALTCI,2,1)-EFALTCJ,2,1))/ZR1 

FCURCNN,1,2)=CEFALTCI,1,2)-EFALTCJ,I,2))/CZN1 
+  3* 

ICZNI+ZNJ) ) 

FCURCNti,2,2)  =  CEFALTCl,2,2)-EFALTCJ,2,2))/ZRl 
FCURCNN,3,2)=CEFALTCI,3,2)-EFALTCJ,3,2))/ZR1 
no  8  5  1 2  =  1, IBB 


RE ADC  37,6  5)RN1,RN2,IBII,IBIJ, 
IBJI, IBJJ 

1  ,ZN2 
f)5 


F0RMATC2F, 41 , 2F) 
NN1=IFIXCRN1 ) 
N  N  2  =  I  F  I  X  C  R :«'  2  ) 
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70 


74 

05 
86 
90 
91 


IF(NN1  .Eq. NN)GO  TO  70 
IF(NN2.Eg. NN)  GO  TO  7  4 
IF(NN1 .EQ.O)GO  TO  86 
GO  TO  8  5 
FCUR(KN,  1 , 2)  =  ( (EFALT( IBJl,  1,2 )-EFALT( IBJJ,  1 ,  2) )/ 
1(ZN2) )+FCUR(NN, 1 , 2) 

GO  TO  85 
FCUR(NN, 1,2)=((EFALT(IBII,1, 2 )- EF ALT( I B IJ , 1,2))/ 
1(ZN2) )+FCUR(NN, 1,2) 

CONTINUE 
REWIND  3  7 
CONTINUE 
REWIND  3  4 
DO  300  NN=1,IBB 

IF( CABS(FCUR(NN, 2, I ) ) .EQ.0)G0  TO  300 


WRITE(38,250)IFLT,NN,FCUR(NN,2,1),FCU 
R(NN, 1,2) 

1  , FCUR(NN, 2, 2)  ,  FCUR(NN, 3, 2) 
250  F0R!1AT(  213  ,8F) 

300  CONTINUE 

DO  400  NN=1 , IBB 

IF( CABS(FCUR( NN, 2, 1 ) ) .EQ. 0)G0  TO  400 

FCUR1=(RKVBA/EIB(IF1X(ABS(REAL(IDELE( 
NN) 

1)  )  )  )  )'''FCUR(NN,  2,1) 


FCUR2=(RKVBA/EIB(IFIX(ABS(REAL(IDELE( 
NN) 

1) ) ) ))*(FCUR(NN, 1 , 2)+FCUR(NN, 2, 2)+FCUR(NN, 3, 2) ) 

FCUR1=FCUR1*FIPU 
FCUR2=FCUR2*FIPU 
IRI=( IFIX(REAL( IDELE(NN) ) ) ) 
IRJ=( IFIX( AIMAG( IDELE(NN) ) ) ) 
CUR1=CABS(FCUR1 ) 
CUR2=CABS(FCUR2) 
IF( CURl .LT. 0.001 )  GO  TO  400 
WRITE(4  1  ,  350)IFLT,  N  \' ,  CURl  ,  CUR2 
350  F0RMAT(2I, 2F) 

400  CONTINUE 

500       CONTINUE 

DO  600  1=1, IBB 
READ(42, 52  5) IL, JL, ALCUR 
IF( IL.EQ.0)GO  TO  610 
IF( REAL(ALCUR) .LE.0.0)GO  TO  600 
BLCUR=CABS( ALCUU) 

CLCUR  =  FIPU*BLCUR-irKVBA/i:iB(JL) 
BLCUR  =  FIPU"BLCUR*R1CVBA/EIB(  IL) 
WRITE(43,525)1L,JL,BLCUR,CLCUR 
5  25       FORMATC 21 , 2F) 
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6  00 
610 

60  2 


C  0 1 1 
l/RI 
W  R I 
FOR 
CLO 
CLO 
CLO 
CLO 
CLO 
CLO 
STO 
END 


TI 
TE 
TE 
tlA 
SE 
SE 
SE 
SE 
SE 
SE 
P 


NUE 
(M 
(43 
T(  ' 
(UN 
(UN 
(UN 
(UN 
(UN 
(UN 


,602) 
,602) 

0  0  0 
IT=36) 
IT=37) 
IT=38) 
IT=3  9) 
IT=A2) 
IT=43) 


0  0') 
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***5k**>k:fcA***AAA**yc  A**  ***>iA*5k*AA)!c**5'{  **********  >c  ****** 

COORD.  FA  ;USi:S  DATA  FROM  1  NPUT  .  F4  ,  F  ALT  .  FA  TO 
OTECTION  DFVICFS  FOR  THE  WHOLE  DISTRIBUTION 

*****  A  *****  Vc  A  ***********************  * 

NIT=30,FILE='RELA.DAT' ) 
NIT=31,FILE='L0CUR.DAT') 


c 

TITLE: 

COORDINATE- 

C 

THE  PR 

SYSTEM. 

AAA*********** 

OPEN(U 

OPEN(U 

OPEN(U 

OPEN(U 

OPEN(U 

OPEN(U 

OPEN(U 

OPEN(U 

DIMENS 

1 , AIDE 

) 

2 ,DEV( 

COMPLE 

*A**A****AAA*A 

C 

READ  I 

A* A AAA******** 

READ(3 

5 

FORMAT 

READ( 3 

10 

FORMAT 

1  1 

READ(3 

12 

FORMAT 

I  F  (  I  .  N 

13 

READ(3 

14 

FORMAT 

I  F  (  I  .  E 

EB(  I)  = 

GO  TO 

15 

DO  20 

READ(3 

16 

FORMAT 

I  F  (  J  .  E 

CURLO( 

20 

CO NT IN 

21 

DO  24 

READ(3 

IF(K.E 

22 

FORMAT 

24 

CONTIN 

A  A  A  A  A  A  * 

*  **  A  A  *  * 

C 

NOW  TO 

A  A  A  A  A  A  * 

*  A  *  A  A  A  A 

23 

TYPE  2 

LCUR.DAT' ) 

FCUR.DAT' ) 

TEST. DAT' ) 

XFOR.DAT' ) 

DEV .DAT' ) 

ADEV.DAT' ) 
ION  EB(30),CUREL(200),CURLO(30),Z(10  0) 
VL(30,A),AIDEV(100,7),XFR(2  00,2),IID(3),RID(100 


NIT=32,FILE= 
NIT=33, FILE= 
NIT=3A , FILE= 
NIT=3  5, FILE= 
NIT=36, FILE= 
NIT=37 ,FILE= 


20,0/A) 

X  Z,AIDEV,C1,C2,Z1,Z2,RID,CUREL,CU1,CURL0,XFR 

AAAAAA****AA*AA***AAA*************A 

N  LOAD  AND  BUS  VOLTAGE  DATA 

A***A***A*AAA*A**AAA**A************AA 

0,5)R 

(F) 

A  ,  10)IB,E1 

(I,F) 

0,12)1 

(I) 

E.O)GO    TO     11 

0,  1A)I,E1 

(I.F) 

Q.0)GO  TO  15 

El 

13 

1=1 , IB 

1  ,  16) J,CU2 

(I,F) 

Q.0)GO  TO  2  1 

J)=CMPLX(CU2*1 .25,0.0) 

UE 

1=1 ,20 

6, 22)K, (DEV( I, J) , J=0, A) 

Q.0)GO  TO  23 

(I,5F) 

UE 

**************  A* *********A*AAAA**** 

SET  ALL  LOAD  DEVICES 

*  A  *******  A  A  *  A  *  A  *******  A  *  *  A  *  *  *  A  *  A  A  A  A  *  A 
5 
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25        FORMATC '  ALL  LOAD  CURRLNTS  ARE  READY, USE  DEVICE 
NUMBER  LIST,  '  ) 

DO  100  1=1 , IB 

IF(CABS( CURLO( I) ) .EO.O)GO  TO  100 

ID=I 
30  TYPE    35, ID 

35        FORMATC'  SELECT  THE  DEVICE  FOR  THE  LOAD  AT 
BUS'  , 13,  '  .  .  ,  ' $) 

ACCEPT  40,IDEV 

AIDEVLCID, 1)=FL0AT(IDEV) 
40        FORMAT(I) 

CURLI  =  REAL(CURLO( I)  ) 

CALL 
SELECT(CURLI,ER(I),A1DEVL(I,1),AIDEVL(I,2),AIDEVL(I,3) 

1  ,  AIDEVLd  ,  A)  ) 
100       CONTINUE 

C         NOW  TO  BEGIN  COORDINATING  ELEMENTS  IN  THE 
DISTRIBUTION 

■ki^-k'k'k-k-k-k-k-k:k-kk'k-kk-k:k-k-ki(i:-kk-k-k-k'k'k-kit^k->'<:k-k-ki(i<'ki<-k-k-kiiki<:k 

110       READ( 34 , 120)NN, I , J, K,R,X 
120       FORMATC 41 , 2F) 

IFCNN. EQ.0)GO  TO  130 

NE  =  NN 

CI  =  I 

CJ  =  J 

RIDCNN)  =  CMPLXCCI,  C.T) 

zCnn)=cmplxCR, X) 

GO    TO     110 
130       DO  140  IE=1,NE 

READC32, 135)1, J, GUI 

135  FORMATC 21 , 2F) 
IFCl.EQ.O)GO  TO  140 
IA=0 

CALL  SERCHCNE,RID,IE,I,J,IA,IID) 

Z1=0.0*Z1 

Z2=0.0*Z2 

DO  1 36  IEL=1 , lA 

Zl  =  l /CZ2  +  ZC 1  IDC IKL) ) ) 

Z2=1/Z1 

136  CONTINUE 

DO  138   IEL=1,IA 

CURELCIIDCIA))  =  CZ2)'>C1/ZCIIDCIEL)))'^CU1 
138  CONTINUE 

IFC  lA.GT.  DTYPE  101,  J 

FORMATC'  CAUTION  THERE  IS  MORE  THAN  ONE 


101 
SOURCE 

140 


lOF  POWER  TO  THE  BUS' 14) 
CONTINUE 


122 


D0150   1=1,  I^E 
READ(35 

l/«5  FORMATC 

IF(  IJ.  E 
CALL  SE 
IA=1 

147 

IA=0 

150  CONTINUE 
*********************** 

C         NOW  TO  DE TERM IN 

CONNECTED 

*********************** 

1 51  DO  190   1=1 ,NE 

IR=IKIX 
IM=IFIX 
IF( AIDE 
IF( AIDE 
GO  TO  1 

155  II=IFIX 

JJ=IFIX 
GO  TO  1 

156  II=IFIX 

JJ=IFIX 

157  CII=II 

CJJ=JJ 
AIDEV(  I 
AIDEV(I 
AIDEV(  I 
IF( AIDE 

lAIDEVLC 11,4)) 
IF( A1DEVL( 11,1) 


,  14  5)IJ,  J,T1NRSH,TIJST1  ,TWST2  ,TS1 

21, 4F) 

Q.0)GO 

RCH(NE, 


TO     151 

RID, I, IJ, J, lA, IID) 


XFR(IID(IA),1)=CMPLX(TS1,T1NRSH) 
XFR( I1D( lA) , 2)=CMPLX(TWST1 ,TWST2) 


********************** 

E  WHAT  ELEMENTS  HAVE  LOAD  BUSSES 

***************************** 


( REAL( R 

(AIMAG( 

VL(1R, 1 

VL(IM,  1 

90 

(REAL(R 

( A1MAG( 

5  7 

( AIMAG( 

(REAL(R 


,  1)  =  CMP 
,  2)  =  CMP 
,  4 )  =  CMP 
VL(II, 1 


ID(I))) 

RID(l))) 

) .NE.O)  GO  TO  155 

) .NE.0)GO  TO  156 

1D(I))) 
RID(I))) 

RID(I))) 

in(i))) 


LX(CJJ, CIl) 

LX(0.0,  AIDEVLdI,  1  )) 

LX(A1DEVL(II,2),AIDEVL(1I,3)) 

)  .GE. 35)AIDEV( I , 5)  =  CMPLX(0 .0, 


.LT.0)AIDEV(I,7)=CMPLX(0.0,AIDEVL(II,A 


)) 

159 

158 
BUS' 


TYPE  158,JJ,II 
FORMATC'  SELECT  DEVICE  FOR 
,  13  ,'  PROTECTING  ELEMENT 

IFEEDING'//'  LOAD  AT  BUS ', 14 ,'..,'$ ) 
ACCEPT  160,IDEV 
160       FORMATC I) 

AIDEV1=FL0ATC IDEV) 
DURR=AIMAGC  AIDEVC I , 4 ) ) 
CURR=REALCAIDEVCI,4)) 

IFC AI M AG C AIDEVC I,2)).GT.0)CURR=10*aCDURR) 
CURR=REALC  CURELC  I ) ) 
CALL  SELECTCCURR,EBCJJ)  ,AI1)EV1,AIDEV2,AIDEV3,CT1) 
CALL 
S  E  T  D  E  V  C  C  U  R  E  L  ,  1  ,  J  J  ,  I  I  ,  A I  D  E  V  1  ,  A I  I)  F.  V  2  ,  A I  D  E  V  3  , 

lAlDEV.XFRC I ,  1 ) ,XFRC I , 2)  , RE, EHC  J  J)  ,EBC II)  , RID, CTl  ,DEV 
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) 

190 
200 


205 

210 

ELEMENT 

230 


EV) 

250 


73,0.0) 

EV) 
300 


400 
4  10 


CT1=0 .0 

1F(RE.NE.O)CO  TO  159 
CONTINUE 
IEL  =  0 

NEL=NE-IEL 

IF(REAL(AinEV(NEL, 1 ) ) .NE.0)GO  TO  300 
1NDEX=1 

I=IFIX(REAL(RID(NEL) ) ) 
J=IFIX(AIMAG(RID(NEL))) 

IF( AIMAG(A1DEV(NEL, 2) ) .NE.0)GO  TO  250 
TYPE  210,J,NEL 
RE  =  0 

FORMAT( '  ENTER  THE  DEVICE    FOR  THE  BUS', 13, 'OF 
,IA) 

ACCEPT  230,IDEV 
FORMATC I) 

IF( IDEV.EQ. 100)GO  TO  400 
AIDEV1=FL0AT( IDEV) 
CURR=AIMAG( CUREL( NEL) ) 

CALL  SELECT(CURR,EB( J) , AIDEVl , AIDEV2, AIDEV3, CTl ) 
IF(AIIJEV1.GT.34)AIDEV(NEL,5)  =  CMPLX(0.0,.414  973) 
CALL  SETDEV(CUREL,NEL,I,J,AIDEV1,AIDEV2,A1DEV3, 
1AIDEV,XFR(NEL,1),XFR(NEL,2),RE,EB(I),EB(J),RID,CT1,D 


205 


IF( RE.NE.0)GO  TO 
TYPE  210, I, NEL 
ACCEPT  2  3  0, IDEV 
AIDEV1=FL0AT( IDEV) 
RE  =  0 

CURR=REAL(CUREL(NEL)  ) 
CALL  SELECTC CURR,EB( I) 


AIDEV1,AIDEV2,AIDEV3,CT1) 


IF(AIDEV1.GT.34)AIDEV(NEL,5)=AIDEV(NEL,5)+CMPLX(.4149 

CALL  SETDEV(CUREL,NEL,I,J,AIDEV1,AIDEV2,AIDEV3, 
1AIDEV,XFR(NEL,1),XFR(NEL,2),RE,EB(I),EB(J),RID,CT1,D 


TO  250 


TO 
TO 


200 
400 


I  F  (  R  E  .  N  E  .  0  )  G  0 

IEL=IEL+1 

IF( lEL. LT.NE)  GO 

IF( INDEX. EQ.0)GO 

IEL=0 

INDEX=0 

GO  TO  2  00 

TYPE  410 

PRINT  4  10 

FORMATC'  DEVICE  LISTING') 

DO  440  1=1 ,NE 

TYPE  420,  I,  ( AIDEV( 1 , J)  ,  J=l  ,  7  ) 
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420 
440 

442 
0.0 

2000 


PRINT  4  20,  I  ,  ( AIDEV( I, J )  , J=l  ,  7) 

WRITE(37,420)I,(AIDEV(I,,J),J=1,7) 

FORMAT ( I,  14F) 

CONTINUE 

WRITE( 37 , 442) 

F0RMAT( '  0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 


1  0. 
CLOS 
CLOS 
CLOS 
CLOS 
CLOS 
CLOS 
CLOS 
STOP 
END 

•k  -k  -k  -k 
k  :k  -k  k 


0  0.0  0.0' ) 
E( UNIT=30) 
E(UNIT=3  1 ) 
E(UNIT=32) 
E(UNIT=3  3) 
E(UNIT=3  4) 
E(UNIT=36) 
E(UNIT=3  7) 


kkkkkk-kkk-kk-k-kkkkkkkkkkkkkkk-kk-kkkkkkkit-k 
■kkkkkkkkkk-kk-kk-kkk-k'J<k'k-k-ki<kkkkk-kkkkkkkkk 


C 

TWO 

c 

^kk-kkkk-kkkkk 

SUB 
DIM 
COM 
IA  = 
DO 


SUBROUTINE  SERCH; FINDS  THE  ELEMENT  CORRESPONDING  TO 


BUSSES. 


kk-kk'k-kkkkkkkkk-kkkkkkkkkkkkkkkk-kk-kkkk-M'kk 

rout: 

ENSI 
PLEX 
0 
30 


INE  SERCH(NE,RID,NN,I,J,IA,IID) 
ON  RID( 100) , IID( 3) 
R I D  ,  C  1  ,  C  2 


20 

3  0 
31 

:k-k-kkkkkk 
Vc  ******* 


IA  = 

IID 

CON 

RET 

END 
**  * 

*  *  * 


IEL=1 , 100 
CJ  =  J 
CI  =  I 

C1=CMPLX(CI , CJ) 
C2=CMPLX(CJ, CI) 
IF(RID(IEL)  .EQ.Cl )G0 
IF(RID( lEL) .EQ. C2)G0 
GO  TO  30 

I  EL 


TO 
TO 


20 
20 


IA+1 
(lA): 

T I  K  U  1 
URN 


*************************************** 
************************************* 


SUBROUTINE  REFALT;READS  FCUR.DAT  TO  OUTPUT  FAULT 


CURRENTS 

C 

************************************************* 


AT  THE  BUSSES 


10 
15 


SUBROUTINE  REFALT(IFALT,NEL,FALT,CUREL) 
READ(3  3 , 15) IFF, NN,TFAL, SEAL 
F0RMAT( 21 , 2F) 


125 


20 


IF(  IE 
1 F  (  I  E 

I  F  (  N  N 
FALTl 
FALT2 
FALT  = 
FLOCU 
I  F  (  F  A 
FALT  = 
HEWIN 
RETUR 
END 


F  . 
F. 

=  T 
=  S 

AM 

R  = 

LT 

AM 

D 

N 


EQ.O 

;j E  .  1 

E.NE 

FAL 

FAL 

INK 

3.0* 

.LE. 

AX1( 

33 


) G  0  TO  2  0 
FALT)GO  TO  10 
L)GO  TO  10 


FALTl , FALT2) 

CUREL 

FL0CUR)FALT=AMAX1( FALTl , FALT2) 

FALTl , FALT2) 


ic  ■fc  -it  -k  -k  -k  ■k  i;  -k  i(  -k  -k  i<  -k  i<  -^  -^  -k  -k  it  -k  :k  i<  -k  it  •k  -k  -k  -f;  y'c  :k  -k  -k  ■><  i<  i<  •/<  -k  -k  i<  -k  i(  :k  iz  :>'<  -i!  'J(  -k  -k  ■!<  :k 

c 

ELEMENT 
C 

DOES 
C 


SUBROUTINE  SETUEV;  SETS  THE  DEVICE  INSTANTANEOUS 

DEPENDING  ON  THE  FAULT  CURRENT  AT  TflAT  BUS 


SETDEV(C 
AIDEV(  10 


30 
SETREL(C 


'k  k  k  ii  i;  k 

SUBROU 
UREL, I 
1  ,XFR1 
DIMENS 
0,7)  ,C 
COMPLE 
IDEV  1  = 
I  F  (  I  D  E 
UREL, I 
lAIDEV 
IF( IDE 
IF( IDE 
IF( IDE 
1),I,  J 


TIN 
EL, 
,XF 
ION 
URE 
X  A 
IFI 
VI  . 
EL, 
2, A 
VI  . 
VI  . 
VI  . 
,IE 


OORDINATE    WITH    OTflER     DEVICES 

ki<-k-k'k-kk'k-k'kkki<'ki<kki!i<i<i<:ki:i<^-kki<k-ki<:l<-k 

E 

I, J, AIDEVl , AIDEV2, AIDEV3, AIDEV 

R2,RE,EB1,EB2,RID,CT,DEV) 

L(10  0),DEV(20,0/4),RID(10  0) 

II)EV,XFR1,XFR2,CUREL,RID 

X(AIDEVl) 

LT.O)CALL 

I, J, AIDEV, AIDEVl, 

IDEV3,XFR1 , XFR2,RE, EBl ,EB2 , Rin, CT, DEV) 

LT.0)GO  TO  100 

EQ.99)GO    TO     100 

LT.35.AND.IDEV1.GT.0)CALL    SETFUS(CUREL(IEL 

L,AIDEV,IDEV1,AIDEV2,AIDEV3,XFR1,XFR2,RE,EB1 


.EB2, 

S  E  T  C  K  B  (  C 
100 

•ki::kkkkk-k 
kic^kk-k^kk 

C 

■k-kkkkk'kk 

SETFUSC  C 
A I DEV (  10 


2  D  E  V  ,  R I  D  ) 

IF(  I  DEVI .GE. 3  5 .AND. I  DEVI .LT. 99) CALL 

UREL(  I  EL)  , 

1I,J,IEL,AIDEV,IDEV1,AIDEV2,AIDEV3,CT,EB1,EB2) 

RETURN 
END 

k  :k  k  k  k  k  k  k  k  k  k  k  ■k  k  -k  'k  -k  k  k  k  •!<  -k  ^  k  k  k  k  -k  k  k  k  k  -k  k  k  -k  -k 
k'kkkkkkk'i!'k'kkki<i(kk'kkkkkk'k:k'k'kitk:ki<-kkkkkkkkk 

SUBROUTINE  S ETF U S ; C OORD I  NATE S  AND  SELECTS  A  FUSE 

k  k  k  k  k  k  k  k  k  k  k  k  "k  k  -k  k  k  -k  k  k  k  k  k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

S  U  B  R  0  U  T  I  N'  E 

Ur>EL,  I  ,  J  ,  lEL,  AIDEV  ,  I  DEVI  ,  AIDEV2  ,  AIDEV3 

1  , XF  R 1  , X  F  R  2 , R  E , E  B  1  , E  n  2 , DE V , R I D ) 

D  I  M  E  K'  S  I  0  N 

0  ,  6  )  ,  CU  R  E  L (  1 0  0  )  ,  D  E  V ( 2  0 , 0 / 4  )  ,  R I D ( 1 0  0 ) 
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20 


25 


35 


40 


A3 
45 


COMPLK 
JJ=IFI 
IF( JJ. 
IF( JJ. 
I  F  (  J  J  . 
I  F  (  J  J  . 
I  F  (  J  J  . 
I  F  (  J  J  . 
IF( JJ. 
IF( IDE 
IF( IDE 
I  F  (  J  J  . 
CI  =  I 
CJ  =  J 
AIDEV( 
AIDEV( 
AIDEV( 
GO  TO 
CI  =  I 
CJ  =  J 
AIDEV( 
AIDEV( 
AIDEV( 
GO  TO 
IF(CAB 
AX=ALO 
IF(REA 
IF( OFS 
I F  (  0  F  S 
GO  TO 
I F  (  0  F  S 
I  F  (  0  F  S 
AIDEV( 
AIDEV( 
AR=AL0 
CALL  F 
TVAL=1 
1F(TVA 
GO  TO 
TYPE  4 
FORMAT 


X  AIDEV, 
X(AIMAG( 
EQ.O)CAL 
NE.O)CAL 
0)CUR 
0)CUR 
0.  AND 
0)C0O 
0  .  A  M  D 
L  E  .  1  5 
GT.  1  5 
0)  GO 


EQ 
NE 
EQ 
NE 
NE 
VI 
VI 
EQ 


XFRl  , 
AIDEV 
L  XFD 
L  XFD 
R=AIM 
R=REA 
.  GOOF 
F  =  (  E  B 
.  GOOF 
)  GAL 
)  GAL 
TO  2 


XFR2, 
(I  EL, 
EV(  1  , 
EV(2, 
AG(  CU 
L  (  C  U  R 
.LT.C 
2  /  E  B  1 
•  LT.C 
L  SFU 
L  SBO 
0 


CUKLL, RID 

2))) 

J, ILL, AIDEV, GOOF) 

J ,  I  EL, AIDEV,  GOOF) 

R  E  L  (  I  E  L  )  ) 

EL( lEL) ) 

URR)GOOF=GURR 

)''=COOF 

URR)GOOF=GURR 

SE(G00F,IDEV1,0FSET,0FSET1) 

RFU(GOOF,IDEV1,OFSET,OFSET1) 


lEL,  1  )  =  CMPLX(CI ,  GJ 
lEL, 2)=AIDEV( lEL, 2 
lEL, 3)=GMPLX(0FSET 
25 


) 

)+GMPLX(FL0AT( IDEVl ) ,0.0) 

, OFSETl ) 


IE 

IE 

IE 

50 

S( 

Gl 

L( 

ET 

ET 

40 

ET 

ET 

IE 

IE 

Gl 

US 

0* 

L. 

50 

5 

(' 


L,  1  )  =  GMPLX(GI,  GJ 
L, 2)=CMPLX(0 .0, F 
L, 4)=GMPLX(0FSET 

XFRl )  .Eq.0)GO  TO 
0( ( AIMAG(XFR1 ) )/ 
AIDEVdEL,  2)  )  .LE 
.LE.  AX)GtJR=10*-'^A 
.LE.AX)CALL     SBOR 


) 

LOAT( IDEVl  ) ) 
, OFSETl ) 

5  0 
3  .0) 

.  15)G0  TO  35 
X 
FU(GL'R,  I  DEV.OFSET,  OFSETl) 


.X 


E(GUR,1DEV,0FSET, OFSETl) 

(  IDE  V),  AIM  AG  (AIDEVdEL,  2))) 

.OFSETl ) 


.LE. AX)CUR=10**A 

.LE.AX)GALL  SFUS 

L, 2)=CMPLX(FL0AT 

L, 3)=CMPLX(0FSET 

0(REAL(XFR2) ) 

E(2  ,  RE  AL(  AIDEVdEL,  2  )  )  ,  AR  ,  OF  SET  1  ,  TV  AL  1  ) 

*TVAL1 

GT. 2 .0)G0  TO  43 


UNABLE  TO  PUT  FUSE  IN  THIS  LOCATION 


SUBSTITUE 

1  ANOTHER  DEVICE'  ) 

RE=1  .0 

GO  TO 
50        CALL  R 

PAUSE 

IF( JJ. 
SFALTl ( AIDEV, G 


300 

EFALT( J ,  I  EL, FALT, CUREL( I  EL) ) 

55 

EQ.O)CALL 

U  R  E  L  ,  R  I  D  ,  1 1:  L  ,  I  ,  J  ,  T  V  A  L  1  ,  D  E  V  , 
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55 


200 

250 

300 

•f;  ■k  "k  »'c  -k 

*  *  *  *  A 

c 

CUKRE 

*  *  *  A  * 


ICVAL 

IF(  JJ 

IFALT 

PAUSE 

ADEV  = 

CALL 

IF(TV 

ADEV  = 

NEGDE 

IF(  AD 

IF(  AD 

GO  TO 

IF(  JJ 

AIDEV 

AIDEV 

GO  TO 

AIDEV 

AIDEV 

RETUR 


) 

,c 

FL 
FU 
AL 
AD 
V  = 
EV 
EV 

5 
.  N 
(I 
(I 

3 
(I 
(1 
li 


E.O) 

VAL) 

0 

OAT( 

SE(  1 

.GT. 

EV+1 

-IFI 

.LE. 

.GT. 

5 

E.O) 

EL,  2 

EL,  4 

00 

EL,  2 

EL,  3 


CALL  SFALT2(AIDEV,IEL,I,J,TVAL1 ,DEV, 


IDEVl ) 

,ADEV,CVAL,OFSET,TVAL) 
TVAL1)G0  TO  200 

X( ADEV) 

15)CALL  SFUSE(C00F,NEGDEV,0FSET,0FSET1) 

15)CALL  SB0RFU(C0OF,NEGDEV,OFSET,0FSETl) 

GO  TO  250 
)=CMPLX(0.0, ADEV) 
)=CMPLX(0FSET,0FSET1) 

)=CMPLX( ADEV, AIM AG (AIDEV ( lEL, 2) ) ) 
)=CMPLX(OFSET,OFSETl ) 


vt  A  * 
•k  k  -k 


E  N  D 

•k  k  -k  k  -k 
k  k  k  k  k 

SUBROUTINE  FU S E ; COM PUTE S  TIME  VALUE  FOR  A  GIVEN  LOAD 


k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 
k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  i<  ■k 


NT 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

SUBROUTINE  FUSE(ID,ADEV,CVAL,OFSET,TVAL) 
TVAL=0.0 


IDEV=IFIX( ADEV) 
1F(I15EV.GT.  15)G0 
X=CVAL-OFSET 
IF( ID.NE. 1)G0  TO 


TO  50 


AO 


(X**3)) 

40 
*3))  + 


50 


70 


TVAL=A.048  0+(-36.4234*X)+(223.7  52  6*(X**2))+(-751.279* 

l+(13A7.01*(X**4))+(-1222.67*(X**5))+(4A0.813*(X**6)) 

TVAL=TVAL-2 .0 

GO  TO  100 

X=CVAL-0FSET+. 1 

T  V  A  L=  6  .  4  2  6  +  (  -  3  7  .  9  2  2  •■''  >  >T  r  1  2  1  .  9  0  9  *  (  X  '<  *  2  )  )  +  (  -  2  1  3  .  4  2  0  *  (  X  * 

1(201. 771*(X**4))  +  (-9 6. 9836 *(X**5))  +  (13.6046*(X**6)) 

TVAL=TVAL-2 .0 

IF( X.GT. 1 .4)TVAL=-1 . 95 

GO  TO  100 

X=CVAL-0FSET 

IF( ID.KE. 1 )G0  TO  70 

TVAL=3 .87698+(-4 . 26  27  68*X)+( 1 . 219340*(X**2)) 

TVAL=TVAL-2 .0 

GO  TO  100 

X=CVAL-0FSET+. 1 
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TVAL  =  A  .6  1A8  4  0+(-6.18322A'''X)  +  (3.2318  1*(X*=«^2))  +  (-.6  0  571 
36*(X**3)) 

TVAL=TVAL-2 .0 


ino 


RETURN 
END 


A  *  A  >'c  A  :!c  A  *  A  ifc  *  ;!:  ;fc  :1c  *  *  A  v^  vk  -^  >^  *  Vc  A  >k  *  *  :%  j'c  A  *  A  A  >V  :fc  X  *  *  v'c  ^fC  rf  -.'c  :k  -Jc  *  A  v'c  y.- 
A*A**5'cAA**AAAA*AAAAAAAAAAAAAA*/C5'cAA*AAAAAAxAA*AAAA 


CURRENT 
C 

A  A  A  A  A  A  A 


10 


SUBROUTINE  RELAY ; COMPUT E S  TIME  VALUE  FOR  A  GIVEN 
VALUE  FOR  A  SPECIFIC  RELAY 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  RELAY  (  NDEV  ,  CURR  ,  TV  AL  ,  TDS  ET  ,  DEV  ,  AIIST  ) 

DIMENSION  DEV(20,0/4) 

TVAL=0.0 

X=ALOG10( CURR) 

IF(X.GT.  1  .2)X=1  . 2 

DO  10  1=0,4 

TVAL=TVAL+DEV(NDEV, I)" (X**I) 

CONTINUE 

IF(NDEV.EQ.7)TVAL=TVAL-2 .0 

TVAL=TVAL+TDSET 

IF(CURR.GT.  AIIST.  AND.  AMST.NE.0)TVAL  =  -2  .0 

RETURN 

END 

AAA  vV  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  Vc  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  -A  A  A  A  A  A  A 

C         SUBROUTINE  CRT  Bl'Jl ;  COMPUT  E  S  TIME  VALUE  FOP.  A  GIVEN 

CURRENT 

C  VALUE 

A  A  A  A  A  A  A  A  A  A  A  A  A  Vt  A  A  A  A  A  A  A  A  A  A  A  A  vV  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 

SUBROUTINE  CKTBKR ( I D , MDEV , D VAL , OF  SET , TV AL , AMVAL) 

CVAL=DVAL-OFSET 

IF(ID.NE.1)G0  TO  50 

CVAL=CVAL+. 1 122 

IF (AMVAL. EQ.0.0)AMVAL=10 

IF( CVAL.GT. AMVAL) GO  TO  9  0 

TVAL=5  .76908  +  (-10.88261*(CVAL))  +  (23.65419*(CVAL'^='<2))  + 
(-34.7922 

I*(CVAL**3))  +  (  3  3.7  7  47  2'''(CVAL**4  )  )  +  (-2  1  .  1  4  5  3  6  *  (  CV  AL*  "  5 
))+(6.4 

231778'KCVAL**6)  ) 

TVAL=TVAL-2 .0 

GO  TO  100 
50        CVAL=CVAL+.23  15 

AMVAL=AMVAL+.07  5 

IF( CVAL.GT. AMVAL)GO  TO  95 

TVAL=5  .85967  2  2  +  (-4  .  394  3  3  5*CVAL)  +  (-  .  50r>90284*(  CVAL*''2) 
)+( 10.40 

12739*(CVAL*'^3)  )  +  (-17  .064956*(CVAL**4)  )  +  (  1  3.01  1  543*(C 
VAL**5)  ) 
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90 

95 

100 


2  +  (- 
TVAL 
GO  T 
TVAL 
GO  T 
TVAL 
RETU 
END 
A  A  *  *  *  *  * 

*  *  *  *  A  *  -/i 

SUBR 


4 .010724 

=  TVAL-2  . 

0  100 

=  -2.0 

0  100 

=AL0G10( .025) 

RN 


^CCVAL'^^'-G)  ) 
0 


>V  *  A  *  A 
A  A  sV  *  Vc 

c 


■k  *  -k-k 
'k  -k  -k  -k 

OUTI 


SUBROUTI 
SETREL(CUREL, lEL 
IXFRl , XF 
DIMENSIO 
AIDEV(  100, 7)  , IID 
1  , DEV(20 
COMPLEX 
COOF=0.0 
CURR=0.0 

cnuR=o.o 

TYPE  5 
F0RMAT( ' 


■k  ■k  -k  -k 
*  A  A  * 

NE  S 

A  *  A  * 

NE 

,I,J 

R2  ,  R 

N 

D(30 

,0/A 

AIDE 


A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 
AAAAAAAAAAAAAAAAAAAAAAAAAAA 

ETREL;C00RDINATES  AND  SELECTS  A  RELAY 

AAAAAAAAAAAAAAAAAAAAAAAAAAAA 

,AIDEV,AIDEV1,AIDEV2,AIDEV3, 
E.EBl ,EB2,RID, CT, DEV) 

),RID(100),IID(3),CUREL(100) 

),FT(30,2) 

V,CUREL,XFR1,XFR2,RID 


5 

YES  0 


DOES  THIS  RELAY  HAVE  AN  INST.  ELEMENT, IF 


') 


ACCE 
FORM 
JJ=I 
MDEV 
IF  (J 
IF(  J 
IF(  J 


PT  7 , IL 

AT(I) 

FIX(AIMA 

=IFIX(-A 

J.EQ.O)A 

J.EQ.O)A 

J.NE.O)A 


G( AIDEV( lEL, 2) ) ) 

IDEVl ) 

IDEV(IEL, 1)=CMPLX(FL0AT(I) ,FLOAT(J)) 

IDEVdEL,  2)  =  CMPLX(0  .0,  AIDEVl  ) 

IDEV(IEL,2)=CMPLX(AIDEV1,0.0)+AIDEV(IEL,2 


IF( 
IF( 
IF( 
IF( 
IF( 
IF( 
CDU 
IF( 
IF( 
IF( 

SELECTC  GOOF 
lAI 
IF( 

SELECT(C00F 
1  AI 
GUR 


JJ.EQ.O)CALL  XFDEV(1,J,IEL,AIDEV,C00F) 

JJ.NE.0)CALL  XFDEV(2 , J, lEL, A I DEV, GOOF) 

JJ.EQ.O)CURR=AIMAG(CUREL(IEL)) 

JJ.NE.O)CURR=REAL(CUREL(IEL)) 

JJ.EQ.0.AND.COOF.LT.CURR)GOOF=CURR 

JJ.NE.O)COOF=(EB2/EB1)*COOF 

R=ALOG10(GOOF) 

JJ.NE.0.AND.GOOF.LT.GURR)COOF=CURR 

CDUR.GT.AIDEV2)A1DEV1=A.0 

CDUR.GT.  AIDEV2.  AND.  JJ.  E().0)GALL 

, EB2, AIDEVl , 

DEV2, AIDEV3, CT) 

CDUR.GT. AIDE V 2 . AND . J J . NE . 0 ) G ALL 

, EBl  , AIDEVl  , 

DEV2, AIDEV3,CT) 

L=AIMAG( GUREL( lEL) ) 
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)))) 


I F  (  J  J 
I F  (  J  J 
1F(  JJ 
IF(  JJ 
I F  (  J  J 
I F  (  1  L 
CALL 
IF(  JJ 
FALTl 
FALTl 
IN  =  0 
DO  10 


8 

10 

CONTI 

DO     30 

.NE.O)  AIDEVdEL,  3)  =  CMPLX(  AII)EV2,  AIDEV3) 

.EQ.O)AIDEV(IEL,  4)  =  CMPLX(  AIDEV2,  AI.DEV3) 

.NE. 0)CURL  =  REAL( CUREL(  lEL) ) 

.NE.O)AIDEV( lEL, 7)  =  CMPLX(CT,0.0)  +  AinEV(IEL,  7) 

.EQ.0)AIDEV(IEL,7)=CMFLX(0.0,CT) 

.NE.O) GO  TO  100 

REFALT(J,IEL, FALTl,  CURL) 

.NE.0)G0  TO  32 

=FALT1*(EB1/EB2) 

=  FALTl/( 10A*(REAL(  AIDEV(  lEL,  4) ) ) ) 

JI  =  1  ,  30 

CALL  SERCH(1,RID,IEL,J,JI,IA,IID) 

IF( IA.EQ.0)G0  TO  10 

DO  8  JA=1 , lA 

I  U'  =  I  N+ 1 
IIDD(IN)=I ID(IA) 

CONTINUE 
NUE 
INN=1 , IN 

I  D=  HDD  (INN) 

IF( ID.EQ. lEL)  CO  TO  30 

IJ=IFIX(AIMAG(RID(ID) ) ) 

CURJ=AIMAG(CUREL(ID) ) 

CALL    REFALTdJ,  ID,  FALT,  CURJ) 

FALT2=ALOG10(FALT) 

AI1)  =  REAL(  AIDEV(  ID,  2)  ) 

IA1D=1FIX(AID) 

IF( IAID.GT.0)GO    TO     16 

IF(IAID.LT.0)FALT0=FALT/(10**(REAL(AIDEV(ID,3 


NDEV=-IFIX(REAL(AIDEV(iD,2))) 

KFAL=1FIX( 10*FALT0) 

DO     15     K=15,KFAL,10 
RK=K*. 1 

IF( ALOGIO(RK) .GT. 1 . 2)G0    TO     16 
CALL 
RELAY (NDEV, RK,TVAL1,REAL(AIDEV(ID,6)) 
1,DEV,REAL(AIDEV(ID,5))) 

CALL 
RELAY(MDEV,RK,TVAL2,TMSET, DEV.O.O) 

TIMEO=(10*"TVALl)+.4 

TIMER=TVAL2 

TIMES=ALOG10(TIMEO) 

TMSET1=TIMES-TIMER 

IF(TMSET1 .LE.0)GO  TO  15 

TMSET=TMSET1+TMSET 

15  CONTINUE 

16  IF( I AID.LT.O)GO  TO  25 
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OFSET  =  AIMAG(  AIDEVdD,  3)  ) 
FNSET=FALT2-0FSET 

IF(FNSET.GT.1.5)FALT2=OFSET+1.5 
KVAL=IFIX( 10*0FSET) 
KFAL=1FIX(10*FALT2) 
DO     20    K=KVAL,KFAL 
RK=K* . 1 
IF(  AID.  LT.  35)  CALL    FIJS  E  (  2  ,  AI  D  ,  RK  ,  OF  SET  ,  TVAL  1  ) 

IF( AID.GE. 35)CALL 
CKTBKR(2,IFIX(AID),RK,0FSET,TVAL1,REAL( 
lAIDEVdEL,  5)  )  ) 

CUR=(  10**RK)/(  10**(REAL(  AIDEVdEL,  3)  ) 
)) 

CALL 
RELAY(MDEV, CUR , TV AL2 , TMSET, DEV , 0 . 0) 

TIME0=1  0**TVM\  +  .^ 
T1MEK=TVAL2 
TiriES  =  AL0G10(TIME0) 
TMSET1=TIMES-TIMER 
IF(TMSET1 .LE.0)GO  TO 
TMSET=TMSET+TMSET1 
ONTINUE 
0  TO  30 


20 


20 
25 
30 


CONTI 
AIDEV 
TMSET 
GO  TO 
32  AID=A 
FALT2 
I F  (  A I 
FALTl 
IF(  AI 
NDEV  = 
KFAL  = 
DIF  =  K 
ID=IF 
I  F  (  I  D 
DO  33 

RK2  =  R 
UREA 


RELAY (NDEV, RK 
lAIMA 


C 

G 
NUE 
(I  EL 
=  0.0 

A5 

I  MAG 

=  ALO 

D.GT 

=  FAL 

D.LT 

-IFI 

IFIX 

FAL- 

IX(D 

.  LE. 

K=l 

R 

K*(E 

L(AI 

I 

C 

,TVA 

G(  AI 

C 

T 

T 

T 


, 6)=CMPLX(0.0, TMSET) 


(AIDEVdEL,  2)  ) 

G10(FALT1 ) 

.0)GO  TO  3  5 

T1*(EB1/EB2) 

.0)FALT  =  FALTl/(  1  0* '"' (  RE  AL(  AI  DE  V  (  I  EL  ,  A  )  )  )  ) 

X(AIMAG(AIDEV(IEL,2))) 

(FALT*10) 

15 

IF/50) 

0)ID=1 

5,KFAL, ID 

K=K*. 1 

B2/EB1)^(10*'UREAL(AIDEV(IEL,A)))/10'''* 

DEVdEL,  3)))) 

F(RK.GT. 20)G0  TO  34 

ALL 

LI,  AI  MAG  (AIDEVdEL,  6)  ),DEV, 

DEVdEL,  5))) 

ALL  RELAY(MDEV,RF;2,TVAL2,TMSET,DEV,0.0) 

IME0=1  0**TVALH-.  A 

IMES=ALOG10(TIMEO) 

1MEK=TVAL2 
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33  CONTI 

34  AIDEV 

35  TMSET 
IF(  AI 
OFSET 
FNSET 
I  F  (  F  N 
KVAL  = 
KFAL  = 
DO  40 

IF(  AI 

IF(  AI 

CKTBKR(2,IFIX 

lAIDE 


40 
45 


CONTI 
AIDEV 
TMSET 
IF(  JJ 
IF(IL 
FALTl 
IF(FA 
I  F  (  I  L 


TMSET 1 
IF(TMS 
TMSET= 
NUE 

(  I  E  L  ,  6  )  =  A 
=  0.0 

D.LT.0)GO 

=AIMAG(AI 

=FALT2-0F 

SET.GT. 1 . 

I F I  X (I  0  *  0 

IFIX( 10*F 

K=KVAL,K 

RK=K* . 

D.LT.35)C 

D.GE.35)C 

(AID) ,RK, 

V(IEL,5)) 

CUR=RK 

CURTA= 

CUR=(E 

CURTA= 

IF(CUR 

CALL  R 

TIMEO= 

TIMER= 

TIMES= 

TMSET 1 

IF(TMS 

TMSET= 

NUE 

(  I  E  L  ,  6  )  =  A 
=  0.0 

.EQ.O. AND 
.EQ.O. AND 
=FALTl/( 1 
LTl .GT.6) 
.EQ.O .AND 


=TIMES-TIMER 
ET1.LE.0)G0  TO  33 
TMSETl+TMSET 

IDEVdEL,  6)  +  CMPLX(TMSET,  0.0) 

TO    45 
DEV(IEL,4)  ) 
SET 

5)FALT2=OFSET+l . 5 
FSET) 
ALT2) 
FAL 
1 

ALL  FUSE(2 , AID, RK, OFSET, TVALl) 
ALL 

OFSET, TVALl , AIMAG( 
) 

-REAL( AIDEV (lEL, 3) ) 
1  0**CUR 

R2/EB1  )*(  1  O'VARK) 

CUR/( 10**(REAL( AIDEV ( lEL, 3) ) ) ) 
TA.LT. 1 .5)G0  TO  40 

ELAY(MDEV,  CURTA,TVAL2  ,T!ISET,  DEV,  0.0) 
10*"TVALl+.4 
TVAL2 

ALOGIO(TIMEO) 
=TIMES-TIMER 
ET1.LE.O)GO  TO  40 
TMSET+TMSETl 

IDEVdEL,  f))  +  CMPLX(TMSET,0.0) 

. FALTl .GT.6)FALT1=6 

.JJ.EQ.0)AIDEV(IEL,5)=CMPLX(0.0,FALT1) 

0*=^REAL(  AIDEVdEL,  3)  )  ) 

FALT1=6 

.JJ.NE.0)AIDEVdEL,5)  =  AIDEV(IEL,5)  +  CMPL 


X(FAL 
100 


Tl, 


10.0) 

RETURN 

END 


C 

LONGEST 

C 
******** 


SUBROUTINE  SF ALT  1  ; PROV I DE S  THE  DEVICE  WITH  THE 

TIME  AT  FAULT, FIRST  DEV  IN  ELEMENT. 
*************************************** 

SUBROUTINE  SF ALT  1 ( AI DE V , CU REL , R I D ,  1  EL ,  I , J , TV AL , DEV , 
ICVAL) 
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8 
10 


)))) 


FUSE( 2, AID, 
1(1 


25 


DIMENSION  AIDEV(100,6),CUKEL(100),RID(100),IID(3), 
1IIDD(30) ,DEV(20,0/4) 
COMPLEX  AIDEV, CUREL, RID 
DO  10  JI=1 , 30 

CALL  SERCH(1,RID,IEL,J,JI,IA,IID) 

IF( IA.EQ.O)GO  TO  10 

DO  8  JA=1,IA 

IN=IN+1 
IIDD(IN)=IID( lA) 

CONTINUE 
CONTINUE 
DO  30  INN=1,IN 

ID=IIDD(INN) 

IJ=IFIX( AIMAG(RID(ID) ) ) 

CURJ=AIMAG(CUREL(ID)) 

CALL    REFALTdJ,  ID,  FALT,  CURJ) 

FALT2=ALOG10(FALT) 

AID=REAL( AIDEV (ID,  2) ) 

IAID=IFIX(AID) 

IF( lAID.CT. 32)GO    TO     100 

FALT2  =  REAL(  AIDEVdEL,  4  )  )+.  5 

FALT0=10**(FALT2) 

IF(IAID.LT.0)FALT  =  FALT0/(10*'-(REAL(A1DEV(ID,3 

NDEV=-IFIX(REAL(AIDEV(ID,2))) 

IF(IAID.LT.O)CALL    RELAY (NDEV, F ALT , TV AL 1 , RE AL 
1(AIDEV(ID,6)),DEV,REAL(AIDEV(ID,5))) 

FAL  =  REAL(  AIDEVdEL,  4)  )  +  .  5 

IF(FAL.LT.FALT2)FALT2=FAL 

IF( IAID.GT.0)CALL 
FALT2 ,AIMAG( AIDEV 
D, 3) ) ,TVAL1) 

IF(TVAL.GE.TVAL1 )G0  TO  25 

GO  TO  30 

IF( IAID.LT.0)CVAL=-FALT2 

CVAL=FAL 

TVAL=TVAL1 


30 
100 

A  *  *  * 
*  *  *  * 
C 

SIDE 


CON 

RET 

END 
******* 

******* 


TINUE 
URN 

*********************************** 
*********************************** 


SUBROUTINE  SFALT2;AS  IN  1, COMPUTES  TIME  FOR  LOAD 


OF 


C  ELEMENT  AT  FALT. 

******************************************** 

SUBROUTINE  SFALT2(AIDEV,IEL,I,J,TVAL1,DEV,FALT,CVAL) 
DIMENSION  AIDEV(10  0,6),DEV(20,0/4) 
COMPLEX  AIDEV 
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FALT2=AL0G10(FALT) 
AID=A1MAG(AIDEV( lEL, 2) ) 
FAL=REAL(  AIDEVdEL,  3)  )+.  5 
FALT=10**(FAL) 

FALT3=FALT/( 10**(AIMAG(AIDEV(IEL,3)))) 
IF( AID.LT.O)CALL 
RELAY (NDEV, F ALT3 , TV AL 1 , AIMAG ( A IDEV ( I  EL 
1,6))  ,DEV,AIMAG(  AIDEVdEL,  5)  )  ) 
IF(FAL.LT.FALT2)FALT2=FAL 

IF(AID.GT.O)CALL  FUSE(2,AID,FAL,AIMAG(AIDEV(IEL,4)) 
1 ,TVAL1 ) 
CVAL=FALT2 
RETURN 


:k  -A  'k  -k  -k  -k 
*  *  *  *  ^  * 


3  0 
32 

EV( 
DEV( 

n 


AO 
50 


END 

•k  ^  ■k  k  k  k 
k  k  -k  k  k  k 

SUBR 
DIME 
COMP 
COOF 
IX  =  0 
IF(N 
DO  3 


kkkkkkkkkkkkkkkkkkkk-kkkkkkkkkkkkkkkk' 
kkkkkkkkk-k-k:kkk'kk-k'kkkkk-kkkkkkkkkkkkkk 

OUTINE  XFDEV(NN,I,IEL,AIDEV,COOF) 
NSION  IID(30) , AIDEV( 100, 6) 
LEX  AIDEV 
=  0.0 


CONT 
DO  3 
CONT 
DO  4 


N  .  N  E  .  1 

0  IA=1 
AI 
IF 
IX 
II 

INUE 

2  K=l, 

INUE 

0  JA=1 
OF 
IF 


)G0  TO  50 

,  100 

D=REAL( AIDEV( lA, 1 ) ) 

( IFIX(Ain) .NE. I)G0  TO  30 

=  IX+1 

D( IX)=IA 

IX 

,IX 

SET=AIMAG( AIDEV(IID( JA) , 3) ) 

(REAL(AIDEV(IID(JA),2)).LT.0)OFSET=REAL(AID 


1IID( JA) , 3) ) 

IF(REAL(AIDEV(IID(JA),2)).GE. 3 5 ) OF SET=REAL( AI 

lIID(jA),5))+OFSET-.05 

IF(  AIMAG(  AIDEVdEL,  2)  )  .LT.O.  AND.  REAL(  AIDEVdl 

l(JA),2)).GT.O.AND.REAL(AIDLV(IID(JA),2)).LT.35) 
20FSET=0FSET+. 1 5 

IF( COOF . GT. 0FSET)G0  TO  40 

C0OF=OFSET 
CONTINUE 
GO  TO  60 

1F(  AIMAG(  AIDEV(  lEL,  2)  )  .  GT.  ())COOF  =  AIMAG(  AI  DEVCIEL,  4  )  ) 
IF(  AIMAG(  AIDEV  (  ILL,  2)  )  .LT.0)COOF  =  REAL(  AIDEVdEL,  4)  ) 
IF(AIMAG(AIDEVdEL,2)).GE.35)C00F  =  AIMAG(AlDEVdEL,5)) 
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-.075 

.0 
60 


1+COOF 

IF(  REAL(  AlDEVdHL,  2)  )  .LT.O.  AND.  AIMAG(  AIDEVdEL,  2)  )  .GT 


1 .AND. 
COOF=C 
D00F=1 
COOF=D 
RETURN 
END 


AIMAG(AIDEV(IEL,2)).LT.35)COOF=COOF+.15 

OOF+.025 

0**(C0OF) 

OOF 


*  A  A  *  A  A  A 
******* 

c 

MAGNETl 

C 

******* 

******* 


:***************** 
:***************** 


************** 
***************** 


******** 
******** 

SUBROUTINE  S ETCKB ; C ORD I  NATE S  AND  SELECTS  A  THERMAL 
C 


A5 


A  6 


******* 
******* 

SUBROU 
1  ,OFSE 
DIMENS 
COMPLE 
CON( 1 ) 
CON( 2) 
CON( 3) 
C0N(4) 
C0N(5) 
J=IFIX 
IF( J.E 
IF(  J 
IF(  J 
IF(  J 
IF(  J 
IF(  J 
IF(  J 
IF(  J 
I F  (  E  B  1 
IF(C00 
GO  TO 
RE2=(E 
IF(RE2 
1  F  (  R  E  2 
CALL  S 
I  F  (  I  D  E 
IF( IDE 
CI  =  I 
CJ  =  J  J 
AIDEV( 
I  F  (  J  .  E 
I  F  (  J  .  N 


BREAK 
******* 

******* 

TINE  SE 
Tl , AM AG 
ION  AID 
X  AIDEV 
=  2  .6 
5 
0 
0 
0 

(AIMAG( 
Q. 0)CAL 
0)CAL 
0)D00 
0)COO 
0)DOO 
0)COO 
0)CUR 
0)CUR 
.KE.EB2 
F.LT. CU 
A6 

B2/EB1) 
.GE. CUR 
.LT. CUR 
BRKR(CO 
V.LT.ID 
V.GE.  ID 


ER. 
********** 

********** 


TCKB( CUREL 
,EB1 ,EB2) 
EV( 100,  7)  ,C0N(5) 
, CUREL 


**************** 

*************** 

,I,JJ,IEL,AIDEV,IDEV,OFSET 


=  3 

=  5 
=  6 
=  7 


AIDE 
L  XF 
L  XF 
F  =  AL 
F  =  l  0 
F  =  AL 
F  =  10 
R  =  AI 
R=RE 
)G0 
RR)C 


V(I 
DEV 
DEV 
OGl 

**D 

OGl 

**T) 

MAG 

AL( 

TO 

OOF 


KL, 1))) 

( 1 , J, lEL, AIDEV, COOF) 

(2, J, lEL, AIDEV, COOF) 

0(COOF)-AIMAG(  AIDEVdEL,  5)  ) 

OOF 

0  (  COOF  )-AIMAG(  AIDEVdEL,  5)  ) 

OOF 

(CUREL) 

CUREL) 

A5 

=  C  U  R  R 


*COOF 

R)C00F=RE2 

R)COOF=CUR 

OF, IDEV3,0 

EV3)1DEV=I 

EV3)IDEV=I 


I  EL,  1  )  =  CMPLX(CI ,  C 
Q.O)AIDEV( I  EL, 2)  = 
E.O)  AIDEVdEL,  2)  = 


R 

FSET, OFSETl , DMMY) 

DEV  3 

DEV  +  1 


J) 

CMPLX(0.0, FLOAT( I DEV) ) 
CMPLX(FLOAT( IDEV) ,0.0)+AIDEV 


1 
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KIEL,  2 
IF( J. EQ 
I F  (  J  .  N  E 
IF( J. EQ 
I F  (  J  .  N  E 
SET=10* 
DO    5  0    K 


50 
60 


CO 
IF 
IF 
IF 
RE 
EN 


NTINU 
(  AHAG 
(  J.EQ 
(  J.NE 
TURN 
D 

*  A  *  *  * 
•k  -k  -k  -k  i< 


) 

.0 
.0 
.0 
.0 

*( 

=  1 

T 
I 

E 

.  N 
.0 
.0 


)AIDEV(IEL,A)=CMPLX(0FSET,0FSET1) 

)AIDEV(IEL,3)=CMPLX(0FSET,0FSET1) 

)FLT  =  A*  AIM  AG  (CURED 

)FLT=4*REAL(CUREL) 

OFSET-.  12) 

,5 

RP=CON(K)*SET 

F(TRP.GT. FLT)GO    TO    60 

E.0)AHAG=ALOG10(CON(K) ) 

)AIDEV( I  EL, 5)  =  CMPLX(0 .0, AM AG) 

)AIDEV( lEL, 5)  =  CMPLX( AMAG, 0 . 0 )  +  AI DEV ( I  EL ,  5 ) 


■k-k-kk-kkkkkkkkkkkkk-k-/f-kk:kk'k^:k-k-k 
■k-kkkkk'k'kkk-kk'k-kkkkk-k-k^^kkkkkkk'kk-k 

SELECT; SELECTS  A  PARTICULAR  DEVICE 


•k-k'kk-k'k-kk-k'kk'k-kk-kif:k-k 
■k^kitkk-kkk-kkkkkkkk-k-k 

C  SUBROUTINE 

**************  >V******A*A***A:V  Is- *AAA5k5^-*^A***  A  A  ******* 


5A 


53 
55 


56 
58 


60 


70 


80 


SU 

ID 

IF 

IF 

IF 

IF 

IF 

GO 

TY 

FO 

1 

AC 

FO 

GO 

IF 

CA 

AD 

CT 

GO 

CA 

AD 

GO 

IF 

CA 

AD 

GO 

CA 

AD 

AD 

AD 


BROUT 
EVL1  = 
(IDEV 
(IDEV 
(  IDEV 
(  IDEV 
(  IDEV 

TO  5 
PE  55 
RMAT( 
DEVIC 
CEPT 
RMAT( 

TO  5 
(EB.G 
LL  SB 
EVL1  = 
=  ADEV 

TO  1 
LL  SF 
EVL1  = 

TO  1 
(EB.L 
LL  SB 
EVL1  = 

TO  1 
LL  SR 
EVL1  = 
E  V  L  2  = 
EVL3  = 


INE  SELECT(CURL0,EB,ADEVL1,ADEVL2,ADEVL3,CT) 

IFIX(ADEVLl) 

LI .EQ.0)GO  TO  99 

LI .EQ. 1)G0  TO  60 

LI .EQ. 2)G0  TO  70 

LI .EQ.4)G0  TO  80 

LI .GT.4)G0  TO  53 

8 

'  NO  SUCH  DEVICE  OR  DEVICE  ERROR,  RE-ENTER 
E  NR  1-4 .  .  , ' $) 
56,  IDEVLl 

I) 

4 

T. 600.0)  GO  TO  53 

RKR(CURLO, IDEVL1,ADEVL2,ADEVL3,ADEVL4) 

FLOATdDEVLl  ) 

L4 

00 

US E (CURLO, IDEVLl ,ADEVL2,ADEVL3) 

FLOAT( IDEVLl ) 

00 

T.600)GO  TO  53 

ORFU(CURLO, IDEVLl, AD EVL2,ADEVL3) 

FLOAT( IDEVLl ) 

00 

ELA( CURLO, IDEVLl, NTAP,NCT,PKUP,CT) 

IDEVLl 

PKUP 

NCT*100+NTAP 
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GO    TO     100 
9  9  A1)EVL1  =  99 

ADEVL2=CURL0 

ADEVL3=CURL0 
100  CONTINUE 

RETURN 

END 

C         SUBROUTINE  SFUSE;USED  TO  SELECT  A  CURRENT  LIMITING 

FUSE 

C  THIS  MAY  BE  EITHER  HIGH  OR  LOV/  VOLTAGE 

SUBROUTINE  SFUSE(CURL0,NDEV,0FSET,0FSET1) 

IF(NDEV.LT.0)GO  TO  10 

DEV=-1 1 . 152  2  1+(2  3.8  6  5*ALOG10(CURL0) )-( 16. 487*( AL0G10( 


CURLO 

CU 

OG 

10 

2*(DEV 

39*( 


1)**2)  )  +  (2 . 9288* ( ALOGIOC CURLO) **3) )  +  (3 . 1 9 6 7 2  * ( ALOG 1 0 ( 

2RLO)**4))-(l .48  87  7*( ALOG10(CURLO)**5) )+(.1776431*(AL 

310(CURLO)**6) ) 

NDEV=IFIX(DEV)+1 

IF(NDEV.GT.15)NDEV=15 

DEV=ABS(FLOAT(NDEV)) 

OFSET=.9532-(.04  6  433*DEV)+(.117  8  85*(DEV**2))+(-.0  2935 

1**3 ))+(.00337*(DEV**4 ))+(-. 000 1847*(DEV**5))+(.00a00 

2DEV**6) ) 

0FSET1=0FSET+. 1 

IF(NDEV.LT.-15)NDEV=-15 

IF(NDEV.GT.15)NDEV=15 

RETURN 

END 

C         SUBROUTINE  SBORFU;USED  TO  SELECT  A  BORIC  ACID  HIGH 

VOLTAGE 

C  FUSE 

SUBROUTINE  SBOKFU ( CURLO , NDE V , OF  SET , OF  SET  1 ) 

ALOGC=ALOG10( CURLO) 

IF( NDEV. LT. 0)G0  TO  10 

DEV  =  88.1f)4  3  +  (-2  15.2235*(ALOGC))  +  (188.225*(ALOGC**2))  + 


(-69. 
**5))+ 


13463*(ALOGC**3))  +  (7.4695*(ALOGC**4))  +  (1  .  6 5 2 5 7  * ( ALOGC 

2(-  .  347786*(ALOGC**6) ) 

NDEV=IFIX(DEV)+1 

IF( NDEV .GT. 1 9)NDEV=1 9 . 


H 
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10        DEV  =  ABS( FLOAT(NDEV)  ) 

IF(NDEV.LT.0)DEV=DEV-15 

OFSET=1.3  6  6111+.113A579*DEV+(-.0012  6  54*(DEV**2))+(.O0 
0021  3* 

1  (DEV**3) ) 

0FSET1=0FSET+. 1 

IF(NDEV.LT.-19)NDEV=-19 

IF(NDEV.GT. 19)NDEV=19 

IF(NDEV.LT.0)NDEV=-NDEV+15 

IF(NDEV.GT.0)NDEV=NDEV+1 5 

RETURN  I 

END  ' 

*********  *A**A**>V  **********************  vk- ******* 

C         SUBROUTINE  SBRKR;USED  TO  SELECT  A  THERMAL  MAGNETIC 

BREAKER 

*************************************************** 

SUBROUTINE  SBRKR( CURL0,NDEV,0FSET,0FSET1, AM AG) 

DIMENSION  C0N(5) 

CON( 1 )=2 .6 

CON(2)=3.5 

CON(3)=4.6 

CON(A)=5.6 

CON(5)=7 .0 

TYPE  10 
10        FORMATC'  ENTER  MANUFACTURER   (N.A.  AT  PRESENT  USE 
FOR  FUTURE 

lEXPANSION) ' ) 

TYPE  25 
25        FORMATC'  INPUT  TEMPERATURE  FOR  AMBIENT  SHIFT, 1  FOR 
LESS  THAN 

110  DEG  C,2  FOR  60    DEC  C , ' / / ' <CR> F OR  NO  SHIFT..') 

ACCEPT  28,IOFF 
28        FORMATC I) 

IFCIOF.EQ. 1 )OFS=-. 1 

IFC lOF.EQ. 2)OFS=.05 

IFCMAN.NE.0)GO  TO  80 

1FCCURLO.LE.600. )G0  TO  30 

IFCCURL0.LE.800. )G0  TO  50 

GO  TO  70 
30        IFCCURLO.lt. 300. )MDEV=1 

IFCCURLO.LT.3O0.)OFSET=2.477  12 

IFC  CURLO.GT. 300.0. AND.CURLO.LE. 3  5  0. )MDEV  =  2 

IFCCURLO.GT.300.0.AND.CURLO.LE.35  0.)OFSET=2.54A07 

IFCCURLO.GT.350.0.AND.CURLO.LE.400.)MDEV=3 

IFC  CURLO.GT. 350.0.AND.CURLO.LE.A5O.)OFSET  =  2.6  0206 

IFCCURLO.GT. 4  0  0.0. AND. CURLO. LE.450.)MDEV=4 

IFCCURLO.GT.40  0.0.A]]D.CURLO.LE.450.)OFSET=2.65321 

IFC CURLO.GT. 4 50.0. AND. CURLO.LE. 500. )MDEV= 5 

IFCCURLO.GT.45  0.0.AND.CURLO.LE.5  0  0.)OFSET=2.6  9  897 
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IF(  CUllLO.G 

IF( CURLO.G 

GO  TO  8  0 

50 

IF(CUHL0.L 

IF(CURL0.L 

IF( CURLO.G 

IF( CURLO.G 

GO  TO  80 

70 

IF( CURLO.L 

1F( CURLO.L 

1F( CURLO.G 

IF( CURLO.G 

80 

NDEV=MDEV+ 

0F5ET=0FSE 

0FSET1=0FS 

TYPE  90 

90 

F0RMAT( '  D 

ELEMENT 

?(0)YES.' ) 

ACCEPT  95, 

95 

FORMATC I) 

FCUR=4*CUR 

SET=10**(0 

DO  100  1=1 

TD 

IF 

100 

COl.^TINUE 

105 

AMAG=AL0G1 

IF(NOEL.NE 

RETURN 

END 

*  ■k  -k  -k  -k  -k  * 

kkkkkkkkkkk 

c 

SUBROUTINE 

k  k  k  k  k  k  k 

kkkkkkkkkkk 

SUBROUTINE 

TYPE  10 

10 

FORMATC'  S 

TYPE  20 

20 

FORMATC '  S 

TABLE' , 

$) 

ACCEPT  25, 

25 

FORMATC I) 

NDEV=-NUM 

IFCCURLO.L 

IFCCURLO.G 

IFC  CURLO.G 

I F  C  C  T  .  E  0  .  1 

IFC  CT.EQ. 3 

IFC  CT. EQ.  1 

TAP=CURLO/ 

T.  500. )MDEV  =  6 

T. 500. )0FSET=2 .77815 

E. 7  0  0. )MDEV=7 

E. 700. )0FSET=2 .8451 

T.700. )MUEV=8 

T. 700. )0FSET=2 .90309 

E.  1000. )MDEV  =  9 

E. 1000. )0FSET=3 .0 

T. 1000. )MDEV=10 

T. 1000. )0FSET=3 . 3 

3A 

T+0FS+. 1 122 

ET+. 12 

OES  THE  DEVICE  HAVE  A  MAGNETIC 

NOEL 

LO 

FSET-. 1 2) 

,5 

=SET*CONCl) 

CTD.GT. FCUR)GO  TO  105 

OCCONCD) 

.0) AMAG=0.0 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

SRELA;USED  TO  SELECT  A  RELAY 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

SRELACCURLO,NDEV,NTAP,NCT,PKUP,CT) 
ELECT  MANUFACTUER   CFUTURE  EXPANSION)') 
ELECT  THE  WESTINGHOUSE  CO  RELAY   , REFER  TO 

NUM 


T. 30)CT=10 

E.30.ANn.CURLO.LT.120)CT=30 

E. 120)CT=120 

0)NCT=1 

0)NCT=2 

20)NCT=3 

CT 


lAO 


71 


72 


73 


7A 


75 


76 


77 


78 


79 


8  0 


81 


8  2 


I  F  (  T  A  P  .  L  E 

.  1 

.  0  )  G  0 

TO 

71 

IF(TAP. LE. 

1 

.  2  )  G  0 

TO 

72 

I  F  (  T  A  P  .  L  E  . 

1 

.5)G0 

TO 

73 

IF(TAP.LE. 

2 

.0)  GO  TO  7A 

IF(TAP.L5. 

2 

.5)  GO 

TO 

75 

IF(TAP.LE. 

3 

.0)GO 

TO 

7  6 

IF(TAP. LE. 

3 

.5)G0 

TO 

77 

IF(TAP.LE. 

A 

.0)GO 

TO 

78 

IF(TAP. LE, 

5 

.  0  )  G  0 

TO 

79 

IF(TAP.LE. 

6 

0)GO 

TO 

80 

IF(TAP.LE. 

7  . 

.0)G0 

TO 

81 

IF(TAP.LE, 

8. 

0)GO 

TO 

82 

IF(TAP.LE. 

10.0)GC 

TC 

)  83 

GO  TO  84 

NTAP=1 

TAP=1  .0 

GO  TO  9  0 

NTAP=2 

TAP=1 .2 

GO  TO  90 

NTAP=3 

TAP=1 .5 

GO  TO  90 

NTAP=A 

TAP  =  2  .0 

GO  TO  90 

NTAP=5 

TAP=2 . 5 

GO  TO  90 

NTAP=6 

TAP=3.0 

GO  TO  90 

NTAP=7 

TAP  =  3  .5 

GO  TO  90 

NTAP=8 

TAP  =  4  .0 

GO  TO  90 

NTAP=9 

TAP  =  5  .0 

GO  TO  90 

NTAP=10 

TAP=6 .0 

GO  TO  9  0 

NTAP=1  1 

TAP  =  7  .0 

GO  TO  90 

NTAP=1 2 

TAP=8 .0 

141 


83 

8A 
90 


GO  TO  90 

NTAP=1 3  . 

TAP=10.0 

GO  TO  90 

NTAP=1 A 

TAP=1 2.0 

PKUP=AL0G10(TAP*CT) 

RETURN 

END 
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C 

A  *  A  A 


3 

5 
10 

15 
20 


25 

30 

40 
45 
50 
55 
60 


70 

75 

80 

400 

410 


420 
440 
445 
4  50 


TIT  LI'  : 

A  A  A  A  A  *  *  A  A  A 

0  P  E  N  (  U 
OPEN( U 

0  P  E  N  (  U 
OPEN( U 
OPEN(U 
EXTERN 
DIMENS 
COMPLE 
READ(3 
FORMAT 
READ(3 
FORMAT 
I F  (  I  .  N 
READ(3 
FORMAT 
IF( I .E 
EB(I)= 
GO  TO 
DO  40 
READ(3 
FORMAT 

1  F  (  I  .  E 
CONTIK 
READ(3 
FORMAT 
READ(3 
FORMAT 
I  F  (  N  N  . 
NE  =  NN 
CI  =  I 
CJ  =  J 

R I  D  (  N  N 
GO  TO 
DO  80 

FORMAT 
I  F  (  K  .  E 
CONTIN 
TYPE  4 
FORMAT 
DO  440 
TYPE  4 
PRINT 
FORMAT 
CONTIN 
TYPE  4 
FORMAT 


PL 

A  A 

NI 
NI 
NI 
NI 
NI 
AL 
10 
X 

2, 
(F 

2, 
(I 

E. 

2, 
(I 

Q. 

El 
15 
K  = 
0. 
(I 
Q. 

UE 
A, 
(I 
4. 
(4 
EQ 


OTD. 
A  A  A  A 

T  =  30 
T  =  3  1 
T  =  32 
T  =  33 
T  =  34 
T25 
N  EB 
AIDE 
3)R 
) 
10)1 

) 

0)GO 
20)1 
,F) 
0)G0 


F4 

A  A 

,F 
,F 
,F 
,F 
.F 
HA 
(3 
V, 


;PL 

AAA 


GTS  DEVICES  FROM  COORD. F4 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

=' ADEV. DAT' ) 
=' DEV. DAT' ) 
='RELA.DAT' ) 
=' FCUR. DAT' ) 
= 'TEST. DAT' ) 


ILE^ 
ILE: 

ILE: 

ILE^ 

ILE  = 

N 

0),DEV(20,0/4),AIDEV(100,7),RID(100) 

RID 


TO  5 
,E1 

TO  25 


1,10  0 

30)I,(AIDEV(I,J),J=1,7) 

,  lAF) 

0)G0  TO  4  5 

50)IB,E1 

,F) 

6  0  )  N  N  ,  I  ,  J  ,  K  ,  R  ,  X 

I,2F) 

.0)  GO  TO  70 


)=CMPLX(CI, CJ) 

55 

1=1 , 20 

READ(31,75)K,(DEV(I,J),J=0,4) 
(I,5F) 

Q.0)GO  TO  400 
UE 
10 
('  DEVICE  LISTING') 

1=1 ,NE 
20,1,  ( AIDEV( I, J)  ,  J  =  l  , 7) 
420,1,  (  AIDEVd,  J)  ,  J  =  l  ,7) 
(I,  14F) 
UE 
50 
('  DEVICES  HAVE  BEEN'  S1:T  IF  YOU  WANT  TO  PLOT 
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GRAPH  I'APEK 

1  ENTER(O) 
ACCEPT  A56 

456       FORMAT(I) 
I  F  (  I  P  .  N  E  .  0 
TYPE  460 
FORMAT('  T 


460 
VOLT 


AGE 


461 
GRAP 

462 


480 
485 

490 


495 
550 


*  *  *  -k 
C 


1  IN  THE  S 
TYPE  4  61 
FORMAT('  I 
H, INTEGER. . ' ) 
ACCEPT  462 
FORMATC I) 
CALL  INITT 
CALL  L0G(I 
CALL  FLUSH 
TYPE  485 
FORMATC'  I 
ACCEPT  490 
FORMATC I ) 
I=IFIXCREA 
J=IFIXCAIM 
EB1  =  EBC  1) 
DO  49  5  K=l 
EB 
CONTINUE 
FAC=ALOG10 
FAC1=AL0G1 
CUR=REALC A 
CUR1=REALC 
CUR2=AIHAG 
CUR3=AIMAG 

■k-k-k-ki<-k-k-ki(i<-k-k-k'k 


, IF  NOT  TYPEC 1 ) ' ) 
.IP 

)G0  TO  480 

HIS  VJILL  PLOT  CURVES  BASED  TO  THE  LOWEST 

YSTEM. ' ) 

NPUT  THE  CURRENT  OFFSET  FOR  THE 

, lAOFF 

AOFF) 


NPUT  ELEMENT   . . ' ) 
,IEL 

LCRIDCIEL))) 
AGCRIDC lEL) ) ) 

.IB 

1=AMIN1 CEBl ,ERCK) ) 

CEBC I)/EB1 )-IAOFF 
OCEBCJ)/EB1)-IAOFF 
IDEVCIEL, 3)  ) 
A1DEVCIEL,4) ) 
CAIDEVCIEL, 3) ) 
CAIDEVCIEL,4) ) 


N0\>/  TO  PLOT  LOAD  SIDE  OF  THE  ELEMENT 


5!r****A**A>!cA**A*****:feA**5!c*A*A***A***A****:(:A:>c***5^**AA 


575 


CALL  R 
IFC AIM 
STL=RE 
STP=FA 
I F  C  A I M 
ST1=AL 
YM=CST 
XM=0.0 
CALL  M 
STP=10 
1TP=IF 
I  F  C  I  T  P 
FORMAT 


EFALTC J, lEL, FALT, 9  0  0.0) 

AGC AIDEVCIEL, 2) ) .GT.O)  GO  TO  585 

ALC AIDEVC lEL, 4) ) 

LT/ 10**STL 

AGCAIDEVC lEL, 5) ) .NE.O )STP= AIM AGC AIDEVCIEL, 5) ) 

OGIOCFALT) 

L+FACl )*190+175.0 


OVEAC  XM, YM) 

*STP 

IXC  S 

.  LT. 

C  '  T 


TP) 

15)TYPE  575 

HE  LOAD  SIDE  RELAY  WILE  NOT  BE  VERY 
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SENSITIVE 

1  TO  A  FAULT'  ) 

I1)EV  =  IFIX(ABS(  AIMAG(  AIDEV(  lEL,  2)  )  )  ) 
TYPE  5  7  0,IDEV,ITP 
570       FORMAT( '  IDEV, ' , 21) 
DIF=ITP-1 5. 
ID=IFIX(DIF/200.0) 
IF(ID.LE.O)ID=I 
TVAL3=10.0 

TDSET  =  AIMAG(  AIDEVdEL,  6)  ) 
DO  580  K=l 5, ITP, ID 

RK=. 1*(FL0AT(K) ) 

CALL  RELAY(IDEV,RK,TVAL,TDSET,DEV) 

IF(TVAL.GT.TVAL3)TVAL=TVAL3 

TVAL3=TVAL 

T=900-( (TVAL+2 .0)^180) 

R=(ALOG10(RK)+FACl+STL)* 190+1 75.0 

CALL  DRAWA(T,R) 
580      CONTINUE 
PAUSE  4 

IF(AIMAG(AIDEV(IEL,5)).NE.0)CALL  DRAWA( (9  00.0) ,R) 
GO  TO  59  1 
585       0FSET=CUR1 

.M=IFIX(80*0FSET) 
ADEV=AIMAG( AIDEV( lEL, 2) ) 
KIM=IFIX( 100*(AL0G10(FALT) ) ) 
DO  587  K=M,KIM 

CVAL  =  K*  .01 
IF(ADEV.LT.35)CALL  FUSE(1,ADEV,CVAL,0FSET,TVAL) 
IF( ADEV. GE. 35)CALL 
CKTBKR( 1 , IFIX(ADEV) ,CVAL,OFSET,TVAL, 
1AIMAG(AIDEV(IEL,5))) 

TVAL=TVAL+2 .0 

IF(TVAL2 .LT.TVAL. AND.ND.NE.0)GO  TO  539 

TVAL2=TVAL 

ND=1 

CALL 
DRAWA(900-( (TVAL)*180),( CVAL+FACl )>M90+1 75.0) 
587       CONTINUE 
589       0FSET=CUR3 

CALL  NOVEA(0 .0,0.0) 
ND  =  0 

M=IFIX(80*0FSET) 
DO  588  K=M,KIM 

CVAL=K*.01 
IF( ADEV. LT.35)C ALL  FUSE(2,ADEV,CVAL,0FSET,TVAL) 
IF( ADEV.GE. 35)CALL 
CKTBKR(2,IFIX(ADEV),CVAL,0FSET,TVAL, 
1AIMAG( AinEV( lEL, 5) ) ) 


1A5 


TVAL=TVAL+2 .0 

IF( TVAL2 .LT. TVAL. AND. ND. NE. 0)GO  TO 
TVAL2=TVAL 
KD=1 
CALL 
DRAWA(900-(  (TVAL)'U80),(  CVAL  +  FACl  )'M90+1  75  .0) 


59  1 


•k  -Jc  -k  i<  ic  -k  :k  it  'k  it  ■/<  -k  ■/<  -k  :!(  ■><  'J<  ■!<  'J:  i<  -k  'J:  it  iK  :k  :k  it  ■)<  i<  ■i<  i< 


588       CONTINUE 

itititititititititititititititititititit 

C         NOW  TO  DRAW  THE  LINE  SIDE 

itititititititititititititititititititititititititititititititit-kititkitititit-kitititkititit-k 


591  CALL  MOVEA(0 
AID=(REAL( AI 
IF( AID.GT.O) 
IDEV=IFIX(-A 
STL=REAL( AID 
STP=FALT/ 10* 
IF( REAL(AIDE 
YM=(STL+FAC) 
XM=0.0 

CALL  MOVEA(X 
STP=STP* 10 
ITP=IFIX(STP 
DIF=ITP-15 
ID=IFIX(DIF/ 
I  F  (  I  D  .  L  E  .  0  )  I 
TVAL3=10.0 
TDSET=REAL( A 
DO  5  90  K=15, 
RK  =  K 
CALL 
1F(T 
TVAL 
CALL 

nRAWA(900-( (TVAL+2 .0 
1+FAC)*190+1 

590       CONTINUE 

IF(REAL( AIDE 

DRAWA(  (900  .0)  , ( ( ALOG 
1STL+FAC)>M9 

595  IF(REAL(AIDE 

0FSET=CUR 
ADEV=REAL( AI 
ND  =  0 

M=IFIX(8  0*0F 
KIM=IFIX( 100 


.0,0.0) 

DEV(IEL, 2) ) ) 

GO    TO     595 

ID) 

EV(IEL, 3)) 

*STL 

V(IEL,5)).NE.0)STP=REAL(AIDEV(IEL,5)) 

'•190+175.0 

M,  YM) 

) 

200.0) 
D=l 

IDEV( lEL, 6) ) 
I  T  P  ,  1  D 
*.  1 

RELAY(IDEV,RK,TVAL,TDSET,DEV) 
VAL.GT.TVAL3)TVAL=TVAL3 
3=TVAL 

)*18  0),((ALOG10(RK)+STL 
75.0)) 

V(IEL, 5) ) .NE.O)CALL 

10(RK)+ 

0+175)) 

V( lEL, 2) ) .LT.0)G0    TO    598 

!)EV(  lEL,  2)  ) 


SET) 

*AL0G10(FALT) ) 
DO     597     K=M,KIM 

=  K*  .01 

4  )  C  A  L  L     F  U  S  E  (  1  ,  A  D  E  V  ,  C  V  A  L  ,  0  F  S  E  T  ,  T  V  A  L  ) 

5)  CALL 


CVAL  = 
IF( ADEV. LE. 3' 
I  F  (  A  D  E  V  .  G  E  .  3 
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CKTBKR(  I  ,  IFIX( ADRV)  , CVAL, OFSET, TVAL, 
1REAL(  AIDEVdEL,  5)  )  ) 
TVAL=TVAL+2 .0 

IF(TVAL2 .LT.TVAL. AND. ND. NE. n)GO  TO  601 
TVAL2=TVAL 
ND=1 
CALL 
DRAWA(9  0  0-(TVAL*180)  , ( CV AL+F AC) *  1 9 0+ 1 7 5 . 0 ) 
597       CONTINUE 
601       0FSET=CUR2 

CALL  MOVEA(0 .0,0.0) 
ND  =  0 

M=IFIX(S0*0FSET) 
DO  599  K=M,KIM 

CVAL=K* .01 
IF(ADEV.LE.34)CALL  FUSE(2,ADEV,CVAL,0FSET,TVAL) 
IF( ADEV.GE. 35)CALL 
CKTBKRC  2, IFIX( ADEV) , CV AL , OF  SET , TVAL , 
1REAL(A1DEV(IEL, 5) ) ) 
TVAL=TVAL+2 .0 

IF(TVAL2 .LT.TVAL. AND. NU.NE.0)CO  TO  598 
TVAL2=TVAL 
ND=1 
CALL 
0  0-(TVAL*180) , ( CVAL+FAC) A190+1 75 .0) 
CONTINUE 


DRAWA(9 
599 
598 
600 


TYPE  600 

F0RMAT(  '  FINISHED?  ? '  ) 

ND  =  0 

CALL  MOVEA(0 .0,0.0) 

CALL  FLUSH 

ACCEPT  610,1 
610       FORMAT(I) 

IF(I.EQ.O)GO  TO  2000 

GO  TO  A45 
2000      CL0SE(UNIT=30) 

CL0SE(UNIT=3 1 ) 

CLOSE(UNIT=32) 

CLOSE(UNIT=33) 

CLOSE(UNIT=34) 

CALL  FINITT(0,0) 

STOP 

END 


SUBROUTINE  REFALT;READS  FOUR. DAT  TO  OUTPUT  FAULT 


c 

current: 
c 

A  A  *  *  *  *  * 


AT  THE  BUSSES 

************  Vc*i!c5!c*A*v't  ******  >V**A>V**** 


****** 

SUBROUTINE  REFALT( IFALT, NEL, FALT, CUREL) 
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10        READ(33,15)IEF,NN,TFAL,SFAL 
15        F0RMAT(2I, 2F) 

IF( IFF.EQ.0)GO  TO  20 

IF(  IFF. NE. IFALT)GO  TO  10 

IF(NN.NE.NEL)GO  TO  10 

FALT1=TFAL 

FALT2=SFAL 

FALT=AMIN1(FALT1 , FALT2) 

FL0CUR=3 .0*CUKEL 

IF(FALT.LE.FL0CUR)FALT  =  AMAX1 ( F ALT  1 , F ALT2 ) 

FALT=AMAX1(FALT1 , FALT2) 
2  0        REWIND  3  3 

RETURN 

ErJD 

C         SUBROUTINE  FU S E ; COMPUTES  TIME  VALUE  FOR  A  GIVEN  LOAD 
CURRENT 

AA*A*  *****************  5'cA**AA*AAAA***AA**AA***AAAAAA*A 

SUBROUTINE  FUSE(ID,ADEV,CVAL,OFSET,TVAL) 

TVAL=0.0 

IDEV=IFIX(ADEV) 

IF(IDEV.GT. 15)G0  TO  50 

X=CVAL-OFSET 

IF( ID.NE. 1 )G0  TO  40 

TVAL=4.048  0+(-36.423A*X)+(22  3.7526A(X**2))+(-751.279* 


(X**3)) 

40 
*3))  + 


l+(1347.01*(X**4))+(-1222.67*(X*"5))+(440.«13*(X**6)) 

TVAL=TVAL-2 .0 

GO  TO  100 

X=CVAL-OFSET+. 1 

TVAL=6  .426+(-37  .922*X)+(1  2  1  .90'J'UX*"-2)  )  +  (-2  1  3.420'- (X* 


l(201.771*(X**4))  +  (-96.9836*(X**5))  +  (18.6  0  46*(X*''^6)) 

TVAL=TVAL-2.0 

IF(X.GT. 1 .4)TVAL=-1 .95 

GO  TO  100 
50        X=CVAL-OFSET 

IF( ID.NE. 1) GO  TO  7  0 

TVAL=3 .876  9  7  9+(-4.26276  8*X)+(1.21940*(X**2)) 

TVAL=TVAL-2 .0 

GO  TO  100 
70        X=CVAL-OFSET+. 1 

TVAL  =  4  .614841  +  (-6.173224'''X)+(3.2318  06*(X**2))  +  (-.6  057 
1 36*(X**3)) 

TVAL=TVAL-2 .0 
100       RETURN 

END 
*************>'{************>■{********************* 

C         SUBROUTINE  RELAY ; COM PUTK S  TIME  VALUE  FOR  A  GIVEN 
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CURRENT 
C 

SUHROU 
DIMENS 
TVAL=0 
X=ALOG 
I  F  (  X  .  G 
DO  10 
TVAL=T 
10  CONTIN 
IF(NDE 
TVAL=T 
RETURN 
END 

c 

curve; 

■k  -k  i<  •'<  *  -k  -k  -k  - 

c 

■k-k:kkkkkk' 


VALUE  FOR  A  SPECIFIC  RELAY 

kkkkkkkkkkkkkkkkkkkkkkk-k-kkkkkkk 

TINE  RELAY ( NDE V , CURR , TV AL , TDS ET , DE V ) 

ION  DEV(20,0/4) 

.0 

10( CURR) 

T. 1 .2)X=1  .  2 

1=0,  A 

VAL+DEV(NDEV,  I)*(X*'^I) 

UE 

V.EQ.7)TVAL=TVAL-2.0 

VAL+TDSET 


TITLE : LOG. LOG, MAKES  LOG  LOG  FORM  FOR  GRAPH  OF  RELAY 


C 

?V  k  k  k  k  k  k  -k 


100 
200 
225 


k  k  k  k 

OUTL 

k  k  k  k 

MUST 

k  k  k  k 

SUBR 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
DO  2 
DO  1 
DN=1 
XM  =  n 
IF(X 
XX=( 
XTT  = 
CALL 
CALL 
CO  NT 
CO  NT 
CALL 
CALL 
DO  A 
DO  3 
D  r^'  =  1 
XM  =  M 


k  k 

IN 

u 

kk 

ou 

E 

M 

D 

D 

D 

D 

M 

D 

00 

00 

0=^ 

*D 

M. 

AL 

90 

M 

D 

IN 

ir: 

M 

D 

0  0 

0  0 

Qk 


k  k  k 

E  B 

k  k  k 

SE 
*  ** 

TIN 
RAS 
OVE 
RAW 
R  AVJ 
RA\J 
RAW 
OVE 
RAW 
N  = 
M  = 
*N 
U 

GT. 
OGl 
0-X 
OVE 
RAW 
UE 
UE 
OVE 
RAW 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

ORDER 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

ALOG.LOG,(?SYS:GRALIB 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

E  LOG(IAOFF) 

E 

A(  0.0, 0.0) 

A(  0.0, 7 5 0.0) 

A(  1000.0, 750.0) 

A(  10  00.0,0.0) 

A  (  0  .  0  ,  0  .  0  ) 

A(900. 0,5.0) 

A(  9 0  0.0, 7  5  0.0) 

0,4 

1  ,  10 


35  0  00)  GO  TO  22  5 

0  (  X  M  )  )  =^  1  8  0 

X 

A(XTT,  175.0) 

R  (  0  .  0  ,  5  7  5  .  0  ) 


VC  X  X  X  X 


k  k  k  k  k 


k  k  k  k  k 


n  = 


A(9 0  0. 0,175.0) 
A(0 .0, 175.0) 
0,2 
1,10 


1A9 


XX=(AL0G10(XM) 
XT=^X+1 75 
CALL  MOVEA(95. 
CALL  DKAWRC805 
300       COMTINUE 
AOO       CONTINUE 

DO  4  60  M=0,3 
XN=175+190*M 
CALL  MOVEA(900 
CALL  NUMBERC 10 
CALL  MOVER(0.0 
J=IA0FF+M 
CALL  NUMBERC J, 
A  60  CONTINUE 

CALL  CHRR0T(90 
CALL  M0VEA(9  60 
CALL  ANCHOSC'C 
CALL  ANCUOSCN 
CALL  ANC1I0S(  H 
CALL  ANCH0S('S 
CALL  M0VEA(9  00 
CALL  DRAWA(9  00 
DO  500  M=0,5 
XN=900-180*M 
CALL  MOVEA(XN, 
IF(M.EQ.0)CALL 
IF(M.EQ. 1)CALL 
IF(M.EQ. 2)CALL 
IF(M.EQ. 3)CALL 
IF(M.EQ.4)CALL 
IF(M.EQ.5)G0  T 
500  CONTINUE 
510  CALL  M0VEA(20. 
NUMBERC 10 
CHRROTCO) 
MOVEAC  200 
ANCHOSC 'T 
ANCHOSC ' I 
ANCHOSC ' C 
MOVEACO .0 


)*19() 

0,XT) 
.0,0.0) 


.0,XN) 
,'I2') 
,10.0) 

'II'  ) 

) 

.0,400.0) 
URRE' , 5) 
T     I'  .4) 
AMP' ,5) 

M) 

.0,0. 0 ) 
.0,750.0) 


15  0.0) 
NUMBERC  .01  ,  ' F3.  2'  ) 
NUMBERC . 1 , ' F2. 1 ' ) 
NUMBERC 1 ,  ' I  1 '  ) 
NUMBERC 10, ' 12' ) 
NUMBERC 100, ' 13' ) 

0  510 

0, 140.0) 
00, ' 14'  ) 

.0,755.0) 
IME     ',5) 
N    S  E  '  ,  5  ) 
ONDS' , 5) 
,0.0) 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

RETURN 

END 

A  A  *  *  *  A  *  A  *  *  *  *  *  v';  5k  *  *  *  *  A  :fe  *  *  A  *  A  *  5k  V:  A  Vc  -,'«  A  *  -A-  rt  «  A  *  *  :^-  A  :^  >'c  ii  >■;  >'c  A 

C         SUBROUTINE  CKTBKR ; COMPUTE S  TIME  VALUE  FOR  A  GIVEN 

CURRENT 

C  VALUE 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBP.OUTINE  CKTHKRCID,  MDEV,  DVAL,  OFSET,  TVAL,  AMVAL) 
CVAL=DVAL-OFSET 


150 


IF( ID. NE. 1 )G0    TO     50 

CVAL=CVAL+. 1 122 

IF(CVAL.GT.  AMVAL.  AND.  AI1VAL.  NE.0)GO    TO    90 

IF( CVAL.LT.0)GO    TO    80 

TVAL  =  5.76908  +  (- 10.88261 *(CVAL) )  + ( 2 3 . 6 54 1 9 'U  CVAL* *  2 )  )  + 
(-34.7922 

l*(CVAL**3))  +  (33.77472*(CVAL**4))  +  (-2  1. 1 4 5 3 6* ( CVAL* *  5 
))+(6.4 

231778*(CVAL**f))  ) 

TVAL=TVAL-2 .0 

GO  TO  100 
50        CVAL=CVAL+. 232 

IF( CVAL.LT.0)GO  TO  80 

AMVAL=AMVAL+.07  5 

IF( CVAL. GT. AMVAL. AND. AMVAL. NE. 0)G0  TO  95 

TVAL=5 .85  9  672  2+(-4.39  433  5*CVAL)  +  (-. 506 9 0 28 4 'K  CVAL* *2 ) 


)+( 10.40 
VAL**5)  ) 

8  0 
90 


9  5 

100 

150 


12739*(CVAL*-'3))  +  (-17.06  4  956*(CVAL'^*4))  +  (13.0  11543*(C 

2+(-4 .010724* (CVAL**6) ) 

TVAL=TVAL-2.0 

GO  TO  100 

TVAL=6 .0 

GO  TO  150 

TVAL=-2 .0 

GO  TO  100 

TVAL  =  ALOG10('.025) 

CONTINUE 


RETURN 
END 

>'c  A  A  *  *  *  *  A  A  :%  *  *  ik  A  >'c  A  >^  *  :!;  *  :!c  v'<  *  A  *  :!c  A  1%  *  *  A  vV  A  *  A  *  *  «  *  *  *  5'c  :1c  *  ;'c  :% 
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APPENDIX  C 


DISK -FILE  FORMATS 
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Disk  file  RES. DAT 


a      Double  precision  wire  resistance  coefficients 
o 


a.     Double  precision  wire  resistance  coefficients 
o 


a      Double  precision  wire  reactance  coefficients 
o 


a..  Double  precision  wire  reactance  coefficients 

b 

Disk  file  TEST. DAT 

element  number/line  bus/load  bus/positive  seq  p.u.  impedance/zero 

seq  p.u.  impedance/transformer  neutral  imp  line  bus/transformer  neutral 

imp  load  bus 

00000 

Load  data 

bus  number/p.u.  VA/per  unit  subtransient  reactance 

0000 

Disk  file  XFER.DAT 

line  bus/load  bus/inrush/low  withstand/high  withstand/transformer  size 

LCUR.DAT 

line  bus/load  bus/current  at  line  bus/current  at  load  bus 

FCUR.DAT 

fault  bus/element  number/3  phase  fault/1  phase  fault 

DEV.DAT 

device  number/a  /a^/a-/a„/a, 
o   1   2   3   A 
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ADEV.DAT 

element  number/element  bus  numbers/device  numbers /line  bus  element 
low  pickup,  high  pickup/load  bus  element  low  pickup,  high  pickup/ 
magnetic  settings  line,  load/time  dial  setting  line,  load/CT  ratios 
line  jload/ 
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APPENDIX  D 


SAMPLE  PROGRAM  DATA  INPUT 
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Before  beginning  the  program  sequence,  RES. DAT  and  DEV.DAT 
contain  the  proper  coefficients  listed  in  Table  8,5,   All  source 
programs  have  been  loaded  and  saved.   PL0TD.F4  uses  subroutines  con- 
tained in  the  PDPIO  system  files  so  it  must  be  loaded  using: 

LOAD  PL0TD.F4,  (3  SYSrGRALIB. 

The  first  program  is  started  with  run  INPUT,   Underlined  words  and 
numbers  are  data  input  by  the  user.   Words  not  underlined  are  prompts 
from  the  computer  in  the  input  sequence  that  follows. 


Do  you  want  to  read  from  a  KYBRD? 

<CR>  for  input  from  TTY,  1  for  read.  .  . 

<CR> 

Enter  the  problem  base  in  MVA 
1 

Enter  the  problem  base  in  KV 
15 

Enter  element  data  as  follows:   Bus  (  )  to  (  ) 
(Length) (Size) (Transformer  impedance  (P .U. )) (Transformer  size 
in  MVA) (Transformer  High  Voltage  in  KV) (Transformer  low 
voltage  in  KV) 

1  2  3000  6  -(element  with  no  transformer) 

2  3  356  .001  .005  1  15  4.1  -(element  with  a  transformer) 
Input  transformer  ZN  in  P.U.  refer  to  BUS 

Side,  (BUS  2)  (BUS  3) .  .  . 
.001  .005  .001  ,005 

(This  is  continued  for  each  element  until  all  elements  of 
the  type  in  RES , DAT  have  been  entered) 

0^  -(tells  computer  no  more  elements) 

At  this  point  if  another  wire  type  is  to  be  used  the 
following  procedure  applies  if  not  skip  this  section 
and  proceed  to  (l) 
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Input  Bus  1  voltage.  .  .0^ 

Data  is  accepted  and  a  list  follows 

(Here  the  computer  types  all  computed  information  on  the 
elements  entered  so  that  they  may  be  checked  for  error.) 

Now  enter  the  load  data  in  the  following  manner;.  .  . 
(Bus) (Load  in  KVA) (Pwr  Factor).  .  .or.  .  . (Bus) (Load  in  HP) 
(Pwr  Factor).  .  .or.  .  . (Bus) (Load  in  Amps) (Pwr  Factor) 
What  method  will  you  use?.  .  .VA,  HP,  or  Amps? 
S^  -(writes  element  data  computed  to  disk) 
End  of  execution 

(Load  RES. DAT  with  the  new  wire  type  coefficients) 
RUN  INPUT 

Do  you  want  to  read  from  a  KYBRD? 
<CR>  for  input  from  a  TTY,  1  for  read.  .  . 
I 

Old  data  is  entered  ready  for  new  line  data 
Enter  element  data  as  follows.  .  .(Program  runs  as 
in  the  first  section) 

(ij)  Picking  up  the  program  from  completion  of  element  data, 

the  data  entry  follows. 

Input  Bus  1  Voltage.  .  .15000 
Input  Bus  2  Voltage.  .  .15000 
Input  Bus  3  Voltage.  .  .4100 

Input 'Bus  8  Voltage.  .  .440 

Data  is  accepted  and  a  list  follows 

(Computer  types  computed  element  data) . 

Number  of  buses  is  8 

Now  enter  the  load  data  in  the  following  manner;.  .  . 

(Bus) (Load  in  KVA) (Pwr  Factor).  .  .or.  ,  . (Bus) (Load  in  HP) 

(Pwr  Factor).  .  .or.  .  . (Bus) (Load  in  Amps) (Pwr  factor) 

What  method  will  you  use?   KVA,  HP,  or  Amps? 

HP 

Enter.  .  . 

4  800  .80 

What  method  will  you  use?.  .  .KVA,  HP,  or  Amps? 

HP 

8  500  .90 

(This  is  continued  until  all  loads  are  entered? 

What   method  will  you  use?.  .  .KVA,  HP,  or  Amps? 
S^ 

End  of  Execution 
RUN  LDFLO 
End  of  Execution 
RUN  ZBUS 
Enter  the  zero  sequence  impedance  for  the  loads  at  each  bus.  .  . 
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8  .001  .005  -for  each  load  where  the  zero  sequence  Impedance 

will  not  be  assumed  to  equal  the  positive  sequence 
'.  impedance 

0^  -tells  the  computer  there  are  no  more  entries 

Enter  the  mutually  coupled  elements  in  per  unit  in  the  following 
format,  TEST. DAT  has  the  element  numbers.  .  ., (Element) (Element) 
(Zero  Sequence  Value) 

2  3  0  .001  (The  zero  must  be  entered  because  the  coupling  is  a 
reactive  value) 
0^  -(no  more  or  no  coupling) 
End  of  execution 
RUN  FALT 

Enter  the  fault  impedance  (Real) (Imaginary) 
0  0 

End  of  execution 
RUN  COORD 

All  load  currents  are  ready,  use  device  number  list 
Select  device  for  the  load  at  bus  4 
4 


Select  the  device  for  the  load  at  Bus  8 

3^ 

Input  temperature  for  ambient  s.hift,  1  for  less  than 

110  Deg  C,  2  for  60  Deg  C 

<CR>  for  no  shift.  .  . 

<CR> 

Does  the  device  have  a  magnetic  element ? (0) Yes 

0 

Select  device  for  Bus  6  protecting  element  feeding  load  at  Bus 

4 

Select  the  Westinghouse  Co.  relay,  refer  to  Table  1^ 

Does  the  relay  have  an  inst.  element  if  yes  0 

0 


(Until  all  backup  devices  have  been  set) . 

Enter  the  device  for  the  Bus  6  of  element  5 
4 

Select  the  Westinghouse  Co.  relay,  refer  to  Table  1^ 
Does  the  relay  have  an  inst.  element  if  yes  0 
0 


(Until  all  devices  in  the  system  have  been  set)  . 
Device  listing  -(Here  the  computer  lists  the  matrix  ADEV) 
End  of  execution 
Run  PLOTD 

Devices  have  been  set  if  you  want  to  plot  graph  paper  0  if  not 
type  1 
0 
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This  will  plot  curves  based  on  the  lowest  voltage  in  the 

system 

Input  the  current  offset  for  the  graph,  integer  1 

Input  element  number 

]_ 

Finished? 

1^   -(no,  0  yes) 

Devices  have  been  set  if  you  want  to  plot  graph  paper 

0  if  not  type  1^ 

1 

Input  element  number 

6 

Finished? 

0^ 

End  of  execution. 

If  a  device  did  not  plot  as  desired  its  settings  can  be  changed 
in  ADEV.DAT. 
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